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NOTATION 


The  reading  of  section  2.1  is  recommended  for  a 
complete  understanding  of  the  notation  below. 

0 = origin  of  U-V  global  coordinate  system 

1 = imaginary  number  (i  = -1)  associated  with  the  V-axis 

when  shaking  force  is  considered 

i,j,k  = orthogonal  unit  vectors 


; = cos  d> 

pq  9pq 


s = sind> 
pq  9pq 

uk  = angular  velocity  of  the  input  link 


m = mass  of  link  pq 

.H  T, 


XD  >y  = local  coordinates  of  the  center  of  mass  of  link 
44  1 4 pq 

kpq  = radius  of  gyration  of  link  pq  about  the  center  of 


mass 


x°  »y°  >k°  = lumped  mass  parameters  of  link  pq  (section 

pq  pq  pq  2.10) 

a = fixed  length  of  link  pq 

W T. 

xrpq’yrpq  “ Vapq  and  yr*apq’  respectively 

I = mass  moment  of  inertia  of  link  pq  about  the  center 
44  of  mass 

Ipq  = mass  moment  of  inertia  of  link  pq  about  point  "p" 


Vll 


i i ^ 


pq 


’m 


pq 


M = 


T 


pq 


t = 


= shaking  force  acting  on  a linkage  due  to  the 
inertia  of  link  pq 

shaking  force  acting  on  a linkage  due  to  all  the 
inertias 


= shaking  moment  about  "0"  acting  on  a linkage  due 
to  the  inertia  of  link  pq 

shaking  moment  about  "0"  acting  on  a linkage  due  to 
all  the  inertias 


- input  torque  of  a linkage  due  to  the  inertia  of 
link  pq 


input  torque  of  a linkage  due  to  all  the  inertias 
reaction  force  acting  on  link  pq  at  point  "p" 
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Major  Department:  Mechanical  Engineering 

This  dissertation  describes  an  approach  to  the 
optimization  of  the  dynamic  behavior  of  single  degree  of 
freedom  planar  linkages,  with  known  geometry  and  motion, 
via  the  digital  computer.  The  first  phase  of  this  approach 
consists  of  the  dynamic  synthesis  of  one  or  two  of  the 
dynamic  properties  considered  in  this  work,  i.e.,  shaking 
force,  shaking  moment,  and  input  torque.  The  result  of 
synthesis  is  a fruitful  reduction  in  the  number  of  free 
design  variables.  The  second  phase  is  essentially  an 
optimization  procedure  where  the  values  of  the  free 
parameters  remaining  from  synthesis  are  determined  such 
that  the  bearing  reactions  and  the  dynamic  properties  not 
considered  in  the  synthesis  phase  are  optimized. 
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An  interactive  CAD  package  for  the  optimization  of  the 
dynamics  of  four-bar  linkages  is  presented,  and  two 
examples  based  on  the  same  mechanism  are  considered.  The 
computer  programs  written  in  APL  language  are  contained  in 
two  appendices. 
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CHAPTER  I 
INTRODUCTION 

1 . 1 Purpose 

Dynamics  of  mechanisms  and  machinery  has  become  a 
field  crucially  important  due  to  the  increasing  demand  for 
machines  with  higher  speeds  and  higher  precision  of  opera- 
tion. Thus,  there  is  a growing  need  for  efficient 
computer  aided  design  (CAD)  techniques  to  optimize  the 
dynamic  behavior  of  such  machines.  As  an  attempt  to 
contribute  to  the  fulfillment  of  this  need,  the  objective 
of  this  work  is  to  meliorate  the  dynamic  synthesis  method 
presented  by  Elliott  in  Ref.  [1],  and  then  use  it  as  a 
basis  for  an  interactive  computer  package  for  the 
optimization  of  the  shaking  force,  shaking  moment,  input 
torque,  and  reaction  forces  of  planar  four-bar  linkages. 

1 . 2 History  and  Literature  Review 
Mechanisms  are  used  everywhere  in  the  modern  world; 
from  simple  everyday  mechanical  items  to  complex  machines: 

toys,  photographic  cameras  and  cassette  recorders 
household  items 

agriculture,  printing  and  textile  machinery 
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packaging  and  assembly  machinery 
machine  tools  and  presses 
vehicles  and  spacecrafts 

prostheses,  robots  and  industrial  manipulators, 

etc . • • • 

Although  this  list  is  far  from  being  complete,  it  is 
sufficient  to  illustrate  the  great  technological  impor- 
tance that  the  engineering  field  of  mechanism  and  machine 
design  has  in  modern  life. 

The  high  level  of  today's  mechanisms  and  machines  is 
a result  of  man's  old  desire  to  understand  and  control  the 
laws  of  nature  and  his  intrinsic  spirit  of  optimization 
which  has  made  him  search  for  better  ways  of  doing  things 
since  his  first  appearance  on  earth. 

The  evolution  of  mechanisms  and  machines  practically 
parallels  the  historic  development  of  mankind  [2,3,4]. 

The  earliest  records  of  the  invention  and  use  of  simple 
mechanical  devices  are  the  prehistoric  drawings,  tools, 
weapons,  traps,  and  other  implements  discovered  by  archae- 
ologists. These  simple  devices  were  created  by  the  trial 
and  error  method  of  experimentation.  In  fact,  this  method 
remained  for  many  centuries  as  the  only  approach  adopted 
by  the  ancient  inventors  until  science  sporadically 
started  to  change  this  situation. 

When  and  where  science  actually  began  is  difficult  to 
determine  [2,5].  The  ancient  Egyptians  probably  were,  to 
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a certain  extent,  familiar  with  the  basic  laws  of  mechani- 
cal equilibrium  without  which  the  pyramids  and  other  grand 
structures  could  not  have  been  built.  The  "science  which 
is  concerned  with  the  motions  and  equilibrium  of  masses" 
was  christened  "mechanics"^  by  the  great  Greek  philosopher 
Aristotle  (384-322  B.C.)  who  summed  up  in  his  Physics  the 
facts  of  mechanics  known  to  the  ancients;  however,  his 
basic  law  of  motion  of  bodies  was  incorrect  1.2,3],  which 
was  found  only  19  centuries  later  by  the  Italian  Galileo 
Galilei  (1564-1642)  who  was  the  first  to  formulate  the  law 
of  falling  bodies.  Proceeding  from  Galileo's  work  and 
from  the  achievements  of  his  contemporaries,  the  great 
English  scientist  Sir  Isaac  Newton  (1642-1727)  discovered 
the  fundamental  laws  of  mechanical  motion  and  the  univer- 
sal law  of  gravitation  and  stated  them  in  a very  clear  and 
concise  form  in  his  Principia  (1687)  [6], 

Galileo,  Newton  and  many  other  scientists  of  the 
sixteenth  and  seventeenth  centuries  laid  down  the  scienti- 
fic foundation  on  which  a new  engineering  was  built. 

Until  then,  engineering  had  been  an  art  based  on  empirical 
rules  passed  from  one  generation  to  another  and  derived 
very  much  from  trial  and  error  experiments  L 4 , 7 ] - 
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The  quoted  definition  of  the  science  of  mechanics  was 
given  by  the  Austrian  physicist  Ernst  Mach  (1838-1916) 
[2,  p.  3]. 
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During  the  eighteenth  and  nineteenth  centuries,  the 
era  marked  by  the  industrial  revolution  and  the  triumph  of 
the  machine,  science  and  engineering  both  grew  rapidly 
[4,8].  Several  new  branches  of  the  science  of  mechanics 
were  created.  In  1875,  the  German  Franz  Reuleaux  (1829- 
1905)  set  the  basis  for  mechanism  science  with  his  pioneer 
book  Theoretisch  Kinematik  which  was  translated  into 
English  in  1876  [9]  and,  in  1888,  L.  Burraester  (1840-1927) 
published  his  classical  book  Lehrbuch  der  Kinematik 
[10].  The  Russian  school  of  mechanism  science  was  also 
created  in  this  period  by  P.L.  Chebyshev  (1821-1894). 

These  two  schools,  the  German  and  the  Russian,  were 
greatly  developed  during  the  first  half  of  the  twentieth 
century.  Although  this  short  period  was  marked  mainly  by 
the  growth  of  these  two  schools,  there  were  also  many 
significant  contributions  to  this  field  from  other  parts 
of  the  world,  e.g.,  Allievi  (Italy),  Rosenauer  and  A.  H. 
Willis  (Australia),  Poeschl  and  Federhofer  (Austria), 
Nicaise,  Koenings,  Bricard,  and  Manheim  (France). 

In  the  United  States,  it  was  not  until  1940  that 
Svoboda  [11]  developed  numerical  methods  to  design  a 
four-bar  linkage  to  generate  a desired  function.  In  1951, 
the  publication  by  Hrones  and  Nelson  [12]  of  an  "Atlas" 
containing  approximately  10,000  coupler  curves  offered  a 
very  practical  approach  for  some  design  engineers.  Un- 
doubtedly, the  publications  from  Prof.  F.  Freudenstein 
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in  the  fifties  played  a crucial  role  in  the  development  of 
mechanism  science  in  the  U.S.  The  research  activity 
related  to  mechanisms  has  accelerated  a great  deal  during 
the  last  three  decades,  particularly  in  the  United  States. 

To  date,  the  scientific  and  engineering  field  of 
mechanisms  has  accumulated  a voluminous  amount  of  knowl- 
edge. As  implied  in  the  above  historical  overview,  this 
is  due  to  the  work  of  creative  minds  from  many  different 
parts  of  the  world.  A great  deal  of  this  knowledge  is 
concentrated  in  the  subfield  of  "Kinematics"  the 
importance  of  which  cannot  be  overemphasized;  a mechanism 
with  poor  kinematic  characteristics  will  not  perform 
satisfactory.  However,  as  stated  in  the  beginning  of  this 
chapter,  high  speed  and  high  precision  of  operation  cannot 
be  achieved  by  a machine  with  poor  dynamic  behavior. 

The  number  of  investigations  into  the  balancing  of 
mechanisms  has  increased  considerably  during  the  last  two 
decades.  One  of  the  earliest  works  in  this  period  was 
done  in  1966  by  Sherwood  [13]  who  suggested  a method  for 
minimizing  the  fluctuation  of  kinetic  energy  of  a four-bar 
linkage  by  using  the  principle  of  dynamically  equivalent 
systems.  Chi-Yeh  Han  [14]  presented,  in  1967,  a method 
for  balancing  the  shaking  force  and  shaking  moment  of  a 
punch-and-reader  unit  by  means  of  a single  counterweight 
attached  to  the  input  crank  with  certain  phase  angle. 


6 


In  1968,  Lowen  and  Berkof  [15]  performed  a survey  of 
the  literature  about  force  and  moment  balancing  for 
mechanisms  with  rigid  body  links.  In  that  same  year, 
Kamenskii  [16]  classified  the  various  known  balancing 
methods  and  established  quality  criteria  for  static,  dyna- 
mic, and  combined  balancing  of  certain  mechanisms.  Still 
in  1968,  Gheronimus  [17]  applied  Chebyshev's  methods  of 
best  approximation  of  functions  to  the  balancing  of 
mechanisms . 

In  1969,  Berkof  and  Lowen  [18]  introduced  the  method 
of  "linearly  independent  vectors"  by  forcing  the  center  of 
mass  of  four-  and  six-bar  linkages  to  stay  stationary. 

This  method  became  the  inspiration  for  several  following 
works.  Also  in  1969,  Sherwood  and  Hockey  [19]  extended 
the  optimization  method  suggested  by  Sherwood  in  1966 
[13],  by  making  it  more  general. 

Skreiner  [20],  in  1970,  discussed  how  the  process  of 
mechanism  design  brings  together  in  one  work  the  various 
considerations  usually  treated  one  at  a time  in  the 
literature . 

In  1971,  Benedict  and  Tesar  [21]  showed  how  to  calcu- 
late the  dynamic  response  of  mechanical  systems  by  means 
of  "kinematic  influence  coefficients";  and  Berkof  and 
Lowen  [22]  applied  the  least-square  theory  to  the  optimi- 
zation of  the  shaking  moment  of  a fully  force-balanced 
four-bar  linkage. 
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In  1972,  Hockey  [23]  proposed  a technique,  based  on 
previous  works  [13,19,24],  for  the  minimization  of  the 
fluctuation  of  input  torque  in  planar  mechanisms  under  the 
influence  of  external  loads.  Tepper  and  Lowen  [25]  stated 
certain  general  theorems  concerning  full  force  balancing 
of  planar  mechanisms  by  redistribution  of  internal  mass. 

In  1973,  Tepper  and  Lowen  [26]  introduced  the  "theory 
of  isomomental  ellipses"  and  showed  that  the  RMS  shaking 
moment  of  unbalanced  planar  mechanisms  is  constant  with 
respect  to  all  points  along  certain  concentric  ellipses. 
Berkof  [27]  demonstrated  how  the  use  of  counterweights  and 
a physical  pendulum  link  can  provide  complete  force  and 
moment  balance  of  an  inline  four-bar  linkage.  Stevensen 
[28]  pointed  out  that  "very  little  of  the  developments  for 
the  balancing  of  mechanisms  may  be  used  directly  in  the 
balancing  of  machines,"  and  presented  a general  method  for 
the  complete  balancing  of  a harmonic  unbalance  of  any 
machine . 

In  1975,  Paul  [29]  presented  a computer  oriented 
overview  of  the  various  techniques  available  for  genera- 
tion of  the  equations  of  motion  of  linkages,  for  integrat- 
ing them  numerically,  and  for  calculating  the  bearing 
reactions . 

Matthew  and  Tesar  [30,31],  in  1977,  considered  the 
synthesis  of  spring  parameters  to  satisfy  specified  energy 


8 


levels  and  to  balance  general  forcing  functions  in  planar 
mechanisms . 

Benedict  and  Tesar  [32],  in  1978,  made  a significant 
contribution  to  mechanism  science.  They  presented  a 
formulation  based  on  velocity  and  acceleration  influence 
coefficients  which  has  proven  to  be  very  efficient  for 
kinematic  and  dynamic  analysis  of  multi-degree  of  freedom 
mechanisms . 

Walker  and  Oldham,  in  1978  [33]  and  1979  [34],  con- 
sidered the  balancing  of  frame  forces  in  planar  linkages 
possibly  containing  prismatic  joints  and  presented  a 
method  for  checking  whether  a linkage  can  be  fully  force- 
balanced  using  counterweights  alone. 

In  1981,  Elliott  and  Tesar  [35]  introduced  a general 
theory,  extracted  from  [1],  that  allows  complex  planar 
mechanisms  to  attain  pre-specif ied  values  of  any  conceiv- 
able dynamic  property.  As  declared  in  section  1.1,  that 
theory  will  be  enhanced  in  this  work  in  order  to  make  it 
more  easily  applicable  to  the  synthesis  of  shaking  force, 
shaking  moment,  and  input  torque  of  such  mechanisms. 

Still  in  1981,  Tricamo  and  Lowen  [36,37]  presented  a con- 
cept for  reducing  the  residual  shaking  forces  of  theoreti- 
cally fully  force  balanced  mechanisms,  and  developed  an 
experimental  balancing  machine  that  reduced  the  actual 
shaking  force  components  of  a four-bar  linkage  by  more 
than  fifty  percent. 
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In  1982,  Tricamo  and  Lowen  [98]  showed  how  the  maxi- 
mum shaking  force  of  a four-bar  linkage  can  be  lowered  to 
a pre-specif ied  value  by  shrinking  the  force  hodograph  of 
the  mechanism  with  the  help  of  two  counterweights.  In 
that  same  year,  they  [39]  extended  the  previous  concept  by 
using  a third  counterweight  to  simultaneously  minimize  the 
maximum  values  of  input  torque,  shaking  moment  and  bearing 
reactions  while  obtaining  the  prescribed  maximum  value  of 
the  shaking  force. 

Freeman  and  Tesar  [40]  presented  in  1982  a technique 
that  allows  the  selection  of  any  generalized  coordinate ( s ) 
of  a complex  mechanical  system  to  be  the  reference  para- 
meter(s)  for  the  development  of  the  system's  dynamic 
model.  Once  this  model  is  obtained,  it  can  be  transferred 
to  any  other  reference  parameter(s)  by  simple  transforma- 
tion equations. 

The  works  [41-57]  represent  some  of  the  other  contri- 
butions to  the  field  of  dynamics  of  planar  mechanisms 
during  the  past  decade. 

In  the  nineteenth  century,  the  industrial  revolution 
considerably  enhanced  man's  physical  powers.  In  the  pre- 
sent century,  a second  industrial  revolution  is  taking 
place,  with  computers  offering  an  enhancement  of  man's 
mental  capabilities.  The  use  of  computers  in  engineering, 
particularly  in  the  field  of  mechanisms,  has  been  rising 
exponentially.  During  the  last  decade,  several 
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general-purpose  computer  programs  for  analysis  and  design 
of  mechanisms  have  been  developed  and  described  in  the 
literature.  Notable  among  these  are  IMP  [58],  DAMN  [59], 
DRAM  [60],  ADAMS  [61],  DADS  [62],  DYMAC  [63],  KINSYN  [64], 
MECSYN  [65],  LINCAGES  [66],  RECSYN  [67],  ANIMEC  [68], 
MEDUSA  [69],  ISD-FORSS-II  [70],  VECNET  [71],  PLANET  [72], 
and  SPACEBAR  [73]. 


CHAPTER  II 
DYNAMIC  SYNTHESIS 


The  objective  of  this  chapter  is  to  provide  a solid 
foundation  for  the  effective  dynamic  synthesis  of  complex 
planar  linkages,  and  then  derive  the  formulations  perti- 
nent to  the  dynamic  synthesis  of  the  four-bar,  the  slider- 
crank,  the  Watt  I and  II,  and  the  Stephenson  I,  II  and  III 
mechanisms . 


2.1  Preliminary  Remarks 

This  work  is  concerned  only  with  mechanisms  having 
rigid  links.  Hence,  the  distance  between  any  two  given 
points  of  a link  will  always  be  considered  invariant  in 
the  derivations  to  follow. 

An  x-y  Cartesian  coordinate  system  will  be  conven- 
iently fixed  to  every  moving  link  of  a mechanism,  as  shown 
by  the  examples  given  in  Fig.  2.1.  Accordingly,  the 
following  comments  are  appropriate  at  this  point: 

1.  An  ordered  pair  of  natural  numbers  will  be  used  to 
designate  a moving  link  and  additionally  indicate 
the  location  of  the  x-y  coordinate  system  of  that 
particular  link.  For  example,  the  origin  of  the 
x-y  coordinate  system  of  representative  link  62  of 
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12 


(b) 


Figure  2.1  Examples  of  Coordinate  System  Placement 
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a certain  mechanism  is  attached  to  point  b of  that 
link,  and  the  positive  x-axis  passes  through  point 
2 of  the  same  link. 

2.  A general  moving  link  will  be  designated  by  "pq," 
where  "p"  and  "q"  have  the  above  meaning. 

3.  The  x-y  coordinate  system  of  link  pq  will 
generally  be  referred  to  as  its  "local"  coordinate 
system. 

4.  The  U-V  coordinate  system  fixed  to  the  ground  of 
the  mechanism  will  analogously  be  called  the 
"global"  coordinate  system. 

2 • 2 Procedure  to  Locate  the  x-y  Coordinate  Systems 

When  choosing  the  locations  for  the  x-y  coordinate 
systems  of  the  moving  links  of  a mechanism,  some  care 
should  be  taken  so  that  the  method  contained  in  this  work 
can  be  adequately  applied.  The  justifications  for  this 
will  be  given  either  explicitly  or  implicitly  in  this 
chapter,  as  they  become  opportune.  For  now,  it  is 
desirable  to  have  some  guidelines  indicating  how  to 
properly  perform  the  above  task. 

Any  moving  link  (binary,  ternary,  quaternary,  etc.) 
and  its  local  coordinate  system  should  correspond  to  one 
of  the  three  cases  shown  in  Fig.  2.2.  As  a complement  to 
that  figure  consider  the  following  remarks: 
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Figure  2.2  The  Three  Types  of  Local  Coordinate  System 
Placement 


Figure  2.3  Equivalence  of  Representations  for  a Link  pq 
With  Prismatic  Joint 
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1.  The  origin  of  the  local  coordinate  system  should 
be  attached  to  the  center  of  a revolute  joint,  if 
there  is  one. 

2.  The  zero  angle  between  the  x-axes  of  links  pq  and 
rq  in  cases  (b)  and  (c)  yields  simpler  analytical 
expressions . 

3.  Since  it  is  unimportant  to  this  work  to  know  which 
link  is  sliding  internally  in  a prismatic  pair, 
the  two  representations  in  Fig.  2.3  are 
equivalent. 

Given  a linkage  mechanism,  the  placement  of  the  x-y 
coordinate  systems  must  be  made  considering  the  whole 
linkage,  and  the  following  rules-of-thumb  should  be 
satisfied: 

Rule  1 : If  a link  has  a grounded  pin  joint,  the  origin  of 

its  local  coordinate  system  should  be  attached  to 
the  center  of  that  joint.  However,  if  the 
grounded  joint  is  prismatic,  the  x-axis  should  be 
parallel  (or  coincident)  to  the  axis  of  that 
sliding  joint  (Fig.  2.1b). 

Ru_le_2:  Each  revolute  joint  of  a linkage  should 

contribute  to  the  placement  of  the  x-y  coordinate 
system  of  at  least  one  of  the  links  jointed  by 
it. 

It  should  be  noticed  that  for  all  the  mechanisms 
shown  in  Fig.  2.4,  Rule  2 is  automatically  fulfilled  if 
Rule  1 is  satisfied.  Analogously,  only  Rule  1 has  to  be 
considered  for  mechanisms  composed  only  of  binary  links 
(Fig.  2.1).  On  the  other  hand,  the  consideration  of 
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Figure  2.4  Six-Bar  Mechanisms  Satisfying  Rules  1 and  2 
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Rule  2 becomes  important  for  some  mechanisms,  as,  for 
example,  the  six-  and  eight-bar  linkages  illustrated  in 
Fig.  2.5. 


2.3  Mass  Parameters  and  Motion  Specification 

of  Moving  Links 

Figure  2.6  shows  the  general  representation  of  a link 
with  planar  motion  which  will  be  used  throughout  this 
dissertation.  Four  "mass  parameters"  will  be  associated 
with  every  moving  link,  i.e.,  the  mass  mnn,  the  radius  of 
gyration  Kpq  about  the  center  of  mass,  and  the  coordinates 
(Xpq,  Ypq)  of  the  center  of  mass. 

To  apply  the  formulations  of  this  chapter  in  the 
dynamic  synthesis  (or  analysis)  of  a linkage,  it  is 
necessary  first  to  know  its  motion  up  to  the  second 
derivative  at  every  considered  instant  of  time.  This 
condition  can  be  satisfied  if,  for  every  moving  link,  the 
position,  velocity,  and  acceleration  of  points  "p"  and  "q" 
are  known  or,  alternatively,  if  <}>pq,  tf’pq  and  $pq  are  known 
in  addition  to  the  motion  of  "p"  (Fig.  2.6). 

The  shaking  force,  shaking  moment,  and  input  torque 
°f  3 linkage  due  to  its  inertia  can  be  calculated  by 
adding  the  individual  contributions  from  each  moving 
link.  The  next  sections  will  deal  with  these  dynamic 
properties  for  the  two  cases  indicated  in  the  above 
paragraph.  First,  expressions  for  links  with  "p-<J>  motion 
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Note  : Joints  indicated  bg  arrans  have  not  contributed  to  the 
placement  at  ang  at  the  local  coordinate  sgstems  . 


Figure  2.5  Watt  I and  an  Eight-Bar  Mechanism  Illustrating 
Rule  2 
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Figure  2.6  A Generic  Link  With  Planar  Motion 
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specification"  will  be  derived,  then  links  with  "p-q 
motion  specification"  will  be  considered. 

These  expressions  will  be  written  in  a form  suitable 
for  dynamic  synthesis.  They  are  essentially  the  same 
expressions  derived  algebraically  by  Elliott  and  Tesar 
[1 >35] . The  objective  here  is  to  develop  a better  insight 
into  their  formulations  to  apply  them  effectively  in  the 
rest  of  this  work. 

2.4  Inertial  Shaking  Force  (Links  With 
p—  <t>  Motion  Specification) 

The  shaking  force  acting  on  a linkage  due  to  the 
inertia  of  link  pq  can  be  expressed  as  (Figs.  2.6  and 


2.7). 


pqypq*pq(cpq+13pq 


where  cpq  = cos*pq  and  spq  = sintpq 


Re-grouping  the  terms  results  in 


(2.4.1) 


where , 
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Figure  2.7  Diagram  of  a Moving  Link  for  Shaking  Force 
Derivation 


Figure  2.8 


Diagram  of  a Moving  Link  for  Shaking  Moment 
and  Input  Torque  Derivations 
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Ipq 

= U +i  V 

P P 

1 

2pq 

mpqxpq 

l 

2pq 

= W'S 

l 

3pq 

= m y 

pqJ  pq 

l 

.. 

3pq 

= -<J>  (c 
?pq 

Equation  (2.4.1)  was  considered  free  of  any  linear 
dependency  in  Refs.  [1,35].  However,  that  is  not  so  since 
^3pq  = ^2pq’  anc^  Eq.  (2.4.1)  must  be  rewritten  as 


s 


l l ,1  1,1 

Y.  D.  +(Y„  +iY,  )D„ 

1 Pq  Ipq  2pq  3pq  2pq 


(2.4.2) 


because  no  linear  dependency  between  the  D-coef f icients 
can  be  allowed  during  synthesis. 


2.5  Inertial  Shaking  Moment  (Links  With 
p— d>  Motion  Specification) 

The  Shaking  Moment  about  the  origin  0 acting  on  a 

linkage  due  to  the  inertia  of  link  pq  can  be  expressed  as 

(Figs.  2.6  and  2.8) 


3M  = (r  +rr,  / )xm  (r  +r~  , ).k  + I <t> 
pq  -p  -G/p'  pq  -p  -G/ p ' - pq9pq 


= m (r  xr  .k)+m  (r^  < xr  .k)+m  (r  xr„,  .k)+I  <j> 

PQ  -P  -P  - Pq  -G/p  -p  - pq  -p  -G/p  pq9pq 
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= m (U  V -V  U ) + 
PQ  P P P p' 


+ m x (V  c -Us  )-m  v (Vs  +U  c ) + 
pq  pq'  p pq  p pq  ; pqypq  P pq  P pqJ 

+ VWpq<Vpq+Vpq,+Wpq*pq('Vpq+Vpq)  * 

+ mpq:pq*pq(-UpSpq+Vpq)-Wpqipq(UpCpqtVpq>  + 

_ 2 _ 2 _2 

+ m (x  +y  +k  )d> 
pq  pq  ^pq  pq'^pq 

Re-grouping  the  terms  results  in 


3M  = Y1  D3  +Y3  D3  +Y 1 D3  +Y1  D3 
Pq  Ipq  Ipq  2pq  2pq  3pq  3pq  4pq  4pq 


(2.5.1 ) 


where 

%q  - 

U V -V  U 
P P P P 

°2pq  ■ 

••  ••  •• 

(V  c — U s ) +d>  (U  c +V  s 
P pq  P pq'  9pq  P pq  P pq 

)+W-V 

°3pq  - 

••  ••  •• 

-(Vs  +U  c ) + d>  (-U  s +V  c 
P Pq  P pq  9pqv  p pq  p 

, • 2 
) -(b  (U 
pq  vpq  p 

°4pq  - 

<|> 

pq 

y! 

4pq 

,-2  _ 2 _ 2 . 
m (x  +y  +k  ) 

pq  pq  ypq  pq' 

Note : 

l 

Yjpq,  j=1, 2, 3 were  defined  in 

section  2 

2.6  Inertial  Input  Torque  (Links  With 
P — 4>  Motion  Specification) 

The  torque  acting  on  the  input  shaft  of  a linkage  due 
to  the  inertia  of  link  pq  can  be  expressed  as  (Fig.  2.6 
and  2.8) 
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1 • • ••••.•• 

Tnn  = [(r  +r„/  ).m  (r  +vn  , ) + 4>  I <t>  ] -s-o> 

pq  -p  -G/p  pq  -p  -G/p  pq  pq  pq J 1 

where  = angular  velocity  of  the  input  link.  Thus, 


Tnn  = ^rn,mnnrn^(rp/nimn,rJ  + (r  -m  ?n/  ) + * I <*>  ] + u 

PQ  “P  Pq-P  “G/p  pq-p  -p  pq-G/p  no  nnTr-,nJ 


pq  pq  pq  i 


PV  * 

••  •• 

• •• 

(-U  S_+V  c ) 

-m  y <f> 

(U  c 

p pq  p pq 

pq'pq  pq 

p pq 

(— U s +V  c ) 

-m  y <P 

(U  c 

P pq  P pq 

pq  pq  pq 

p pq 

+ m x 


" mpqXpqVq(%Cpq+^pSpq)+inpqypq*pq(^pSpq'^pCpq)  + 

_ 2 _2  _2 

+ m (x  +y  +k  )<f>  <f>  1-5-u 

pq  pq  apq  pq  pq  pqJ  i 

Re-grouping  the  terms  results  in 

1 1 5 1 5 1 5 1 5 

pq  '1PqD'1Pq+^2pqD2pq+Y3pqD3pq+Y4pqD4pq 

where 

5 


D 


D 


i pq 

5 

2pq 


(U  U +V  V )tio. 

P P P p'  1 

• ••  ••  ••  # 

[<i>  (-U  s +V  c ) + <j)  (-U  s +V  c ) 

pq  p pq  p pq'  p pq  VpCpq; 

• 2 

- <J>  (U  C +V  S )]-S-u). 

pq  p pq  p pq  J i 


(2.6.1 ) 


D 


3pq 


• " • « 

[-<)>  (U  c +V  s )-<(»  (U  c +V  s ) + 
pq  p pq  p pq'  pq^  pcpq  pSpq' 


.2  . 

+ ( U s -V  c ) ]*u> . 

pq  P pq  P pq  J i 
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5 • 

°4pq  = *pq*Pq*“i 


2.7  Inertial  Shaking  Force  (Links  With 
p-q  Motion  Specification) 

As  was  shown  in  [1,35j,  the  shaking  force  acting  on  a 
linkage  due  to  the  inertia  of  link  pq  can  be  expressed  as 


’F  = Z y2  D2  ; i = 1 to  4 

ipq  ipq 


(2.7.1 ) 


where 


Ipq 

2 

2pq 

2 

3pq 

2 

4pq 


m (1-  x fa  ) 

pq  pq  pq 


m (-y  fa  ) 

pq  °pq  pq 


m (y  fa  ) 

pq  w pq  pq 


m (x  fa  ) 

pq  pq  pq 


(2. 7. 2. a) 
( 2 . 7 . 2 . b ) 
(2.7.2.C) 
(2.7.2.d) 


(a  = length  of  link  pq  = distance  between  points 
"p"  and  "q") 


and 


D 


D 


D 


D 


ipq 

2 

2pq 

2 

3pq 

2 

4pq 


U +iV 
P P 


-V  +iU 
P P 


-V  +iU 

q q 


U +iV 

q q 


(2. 7. 3. a) 
(2.7.3.b) 
(2.7.3.C) 
(2 . 7 . 3 . d ) 
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Equation  (2.7.1)  was  also  regarded  in  [1,35j  as  containing 

, . 2 2 
no  linear  dependency.  But,  since  D0  = iD.  , and 

2 2 2pq  1pq 
°3pq  = iD4pq»  it  must  be  modified  to 


(Y1pq+lY2pq)D1pq+(Y4pq+iY3pq)D4pq 


(2.7.4) 


which  is  now  in  a form  suitable  for  synthesis. 

Clearly,  all  the  Y's  above  have  units  of  mass,  and 
2 2 

D1pq  anc*  D4pq  are  ^he  equivalent  accelerations  of  points 
"P"  and  "q,"  respectively.  Therefore,  according  to  Eq. 
(2.7.4),  the  mass  nip^  of  link  pq  can  be  transferred  from 
the  center  of  gravity  to  points  "p"  and  "q"  as  "composite" 
point  masses  expressed  in  complex  form.  As  depicted  in 
Fig.  2.9,  this  new  system  of  two  composite  point  masses  is 
exactly  equivalent  to  the  original  link  model,  for  shaking 
force  considera tions . 

The  concept  of  composite  point  masses  allows  one  to 
write  the  expression  of  the  total  shaking  force  acting  on 
a linkage,  in  a single  step,  without  linear  dependencies , 
by  just  looking  at  the  diagram  of  the  mechanism  (Fig. 

2.10) . 


The  mass  mpq  at  the  center  of  gravity  of  a link  pq 
was  previously  transferred  to  points  "p"  and  "q"  by  means 
of  the  four  dimensionless  "transfer  coefficients"  inside 
the  parentheses  in  Eqs.  (2.7.2).  Analogously,  any  point 
mass,  be  it  composite  or  not,  can  be  transferred  by  means 
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Figure  2.9  Two  Shaking  Force  Equivalent  Link  Models 


(Y 


^+1  ^Y313  + Y323^ 


^Y113  +'Y123^+1^y2 


313  4 Y323  ^ D113 


N0tB  ’ D1I2  3 Df23  * and  D113  " D123 


Figure  2.10  Shaking  Force  Derivation  for  a Four-Bar 
Linkage 
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of  the  same  coefficients  if  Xp^  and  yp^  are  replaced  with 
the  coordinates  of  the  point  where  the  mass  to  be 
transferred  is  located.  This  is  exemplified  for  a 
composite  point  mass  (A+iB)  in  Fig.  2.11. 

Now  it  can  be  understood  why  rule  2 given  in  section 
2.2  must  be  satisfied  for  the  six-  and  eight-bar 
mechanisms  in  Fig.  2.5.  If  that  rule  is  overlooked,  the 
minimum  number  of  composite  point  masses  with  linearly 
independent  D-coeff icients  may  not  be  obtained  (Fig. 

2.12).  Although  Fig.  2.12  is  based  on  shaking  force  only, 
the  same  justification  holds  for  shaking  moment  and  input 
torque,  as  will  become  apparent  in  the  following  sections. 

2‘8  inertial  Shaking  Moment  and  Inertial  Input  Torque 
(Links  With  p-q  Motion  Specification) 

As  was  shown  in  [1,35],  the  shaking  moment  about  the 
origin  0 acting  on  a linkage  due  to  the  inertia  of  link  pq 
can  be  expressed  as 


s 


(2.8.1 ) 


where 


(2. 8. 2. a) 


3 


2 9 2 9 


(2.8.2. b) 


3 


(2.8.2.C) 
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Figure  2.11  Transformation  of  a Composite  Point  Mass  at 
"r"  into  Two  Composite  Point  Masses  at  "p" 
and  "q" 
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Figure  2.12  Examples  of  the  Influence  of  Local  Coordinate 
Systems  on  the  Number  of  Composite  Point 
Masses 
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v 3 

/ — 2 -2  -2  . 2 
m (x  +y  +k  Ua 

pq  pq  ^ pq  pq  pq 

Y4pq" 

(2.8.2.d) 

and 

^4 

Ipq  = 

U V -V  U 
P P P P 

(2. 8. 3. a) 

°2pq  = 

(UV-VUJ+(U  v -v  5 ) 
pqpq  q P q P 

(2.8.3.b) 

o4, 

3pq 

(U  U +V  V )-(U  U +V  V ) 
P q P q V q P q p; 

(2.8.3.C) 

4 

D, 

4pq 

U V -V  u 

q q q q 

(2.8.3.d) 

It  was  also  shown  in  [1,35]  that  the  torque 

acting  on  the 

input  shaft  of  a linkage  due  to  the  inertia 

of  link  pq  can 

be  expressed  as  (Appendix  A demonstrates  by 

means  of 

kinematic 

influence  coefficients): 

:t 

pq 

= y.  D?  ; i = 1 to  4 

L ipq  ipq  ’ H 

(2.8.4) 

where 

°1pq  = 

(U  U +V  V )i(1). 

p p p p • “l 

(2. 8. 5. a) 

°2pq  = 

t(UnU  +VnV  )+(U  U +V  V )]*u. 
pq  pqqp  qpJ,“i 

( 2 . 8 . 5 . b ) 

3pq 

[-(UnV  -V  U ) + (U  V -V  U )Um. 

p q pqqp  q p J ,aji 

(2.8.5.C) 

D® 

4pq 

(2. 8 . 5 . d ) 

A glance  at  Eqs.  (2.8.2)  shows  that  all  the  Y's  have 
units  of  mass.  Thus,  examining  Eq.  (2.8.1)  and  Eqs. 
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(2.8.3  3),  it  can  be  concluded  that  for  shaking  moment  and 
input  torque  considerations,  the  mass  at  the  center  of 
gravity  of  a link  pq  and  the  moment  of  inertia  about  that 
point  can  be  replaced  with  four  point  masses  located  at 
"p,"  "q,"  and  two  other  points  whose  actual  positions  are 
immaterial  to  this  work.  It  should  be  noticed  though, 
that  point  masses  Y|pq  and  Yjpq  are,  in  general, 
continuously  changing  their  positions  on  link  pq,  as 
opposed  to  point  masses  Y^pq  and  Y^pq  which  are  always 
located  at  "p"  and  "q,"  respectively.  This  fact  is 
represented  by  the  two  broken  lines  in  Fig.  2.13. 

It  is  important  to  observe  that  the  dimensionless 
coefficients  that  multiply  mpq  in  Eqs.  (2.8.2)  can  be  used 
to  transfer  any  point  mass  of  a link  pq  to  the  same 
locations  of  point  masses  Y?pq  (i  = 1 to  4)  just  by 
removing  Kpq  from  those  coefficients  and  replacing  x and 
Ypq  the  coordinates  of  the  point  mass,  as  shown  in 

Fig.  2.14. 

If  a link  pq  is  grounded  by  a pin  joint,  the  point 
mass  Y1pq  located  at  "p"  is  motionless  (o!jpq  = Dipq  = °) 
and  should  not  be  considered;  in  addition,  for  input 
torque  considerations,  the  point  masses  Y|pq  and  Y|pq  are 
also  located  at  "p"  (D2pq  = D^pq  = 0),  and  should  be 
disregarded,  too.  If  a ground  link  has  constant  angular 
velocity,  further  reduction  of  the  number  of  linear 
independent  D— coefficients  will  occur.  Thus,  for  shaking 
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Figure  2.14  Redistribution  of  a Point  Mass 
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4 4 

moment,  D2pq  = D4pq»  as  demonstrated  in  Fig.  2.15;  and  for 
input  torque,  D4pq  = 0.  Finally,  D2pq  = = 0 if  the 

origin  of  the  global  coordinate  system  (the  shaking  moment 
reference  point)  coincides  with  the  center  of  the  grounded 
revolute  joint.1  Figure  2.16  shows  the  three  types  of 
grounded  link  conf igurations  and  how  they  decrease  the 
number  of  linearly  independent  D-coeff icients  in  the 
expressions  for  the  dynamic  properties  of  a linkage 
mechanism. 

Figures  2.17  and  2.18  illustrate  shaking  moment  and 
input  torque  derivations,  respectively,  for  a four-bar 
linkage;  and  Figure  2.19  is  a collection  of  all  the  basic 
formulations  considered  thus  far. 

Except  for  the  slider-crank,  all  the  mechanisms  that 
are  going  to  be  considered  in  this  work  have  only  revolute 
joints.  Hence,  the  contents  of  the  last  three  sections 
will  be  very  useful  for  the  rest  of  this  chapter.  The 
expressions  derived  in  sections  2.4,  2.5,  and  2.6  for 
links  with  p-4>  motion  specification  were  expressed  in  a 
generic  form.  Therefore,  with  the  help  of  Appendix  C one 
can  easily  derive  shaking  force,  shaking  moment,  and  input 
torque  expressions,  free  of  linear  dependencies , for 
linkages  containing  prismatic  joints. 


1 


Appendix  B explores  this  fact  further. 
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Figure  2.15  Grounded  Link  With  Constant  Angular  Velocity 
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Figure  2.17  Shaking  Moment  Derivation  for  a Four-Bar 
Linkage 


V 

0 * 


Figure  2.18 


Input  Torque  Derivation  for  a Four-Bar 
Linkage 


Figure  2.19  Collection  of  Basic  Formulations 
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2.9  The  Method  of  Dynamic  Synthesis 
The  synthesis  method  to  be  used  in  this  work  is  the 
one  presented  in  Refs.  [1,35J,  whose  description  is  given 
below. 

According  to  the  previous  sections,  the  three  dynamic 
properties  that  were  considered  can  be  expressed  in  the 
following  form, 

dP  = [fl(Y's)](D)1+[f2(Y's)](D)2+  ...  + [fn(Y's)](D)n 

(2.9.1 ) 

where  the  contents  of  the  brackets  multiplying  the  n 
linearly  independent  D-coeff icients  are  real-valued 
functions  of  the  Y's,  when  either  shaking  moment  or  input 
torque  is  considered,  and  complex-valued  functions  of  the 
Y's,  if  shaking  force  is  dealt  with. 

For  any  linkage  mechanism  with  known  geometry  and 
motion,  the  D-coeff icients  can  be  calculated  for  any 
position  of  the  input  link.  Thus,  if  dP  is  pre-specif ied 
at  n distinct  values  of  the  input  angle  along  one  cycle  of 
operation,  the  following  matrix  equation  can  be  formed 

{dP}  - = [D]  { f ( Y ' s ) } , (2.9.2) 

1 Jnx1  Jnxnl  J nxl 

Since  the  square  matrix  above  is  generally  non-singular , 


the  result  is 
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{ f ( Y ' s ) } = [D]  1{dP]  (2.9.3) 

When  the  n numerical  values  obtained  from  the  right- 
hand-side  of  the  above  equation  are  substitued  into  Eq. 
(2.9.1),  the  synthesis  of  dP  is  accomplished.  If  the 
result  is  not  satisfactory,  the  original  set  of  synthesis 
specifications  can  be  modified  and  the  whole  process  re- 
peated again. 

Matrix  equation  (2.9.3)  gives  n linearly  independent 
algebraic  equations  (2n  for  shaking  force)  containing  the 
mass  parameters  of  the  linkage.  They  can  be  used  as 
equality  constraints  to  reduce  the  number  of  design  vari- 
ables when  optimizing  the  behavior  of  the  dynamic  proper- 
ties not  considered  in  the  synthesis  and/or  the  magnitude 
of  the  reaction  forces  at  the  joints. 

In  Refs.  [1,35],  it  was  stated  that  it  is  always 
possible  to  concurrently  satisfy  both  shaking  force  and 
input  torque  specifications.  However,  this  assertion  is 
not  true  since  for  some  mechanisms  the  set  of  equations 
obtained  from  shaking  force  synthesis  is  not  linearly 
independent  of  the  set  of  equations  obtained  from  input 
torque  synthesis,  as  will  be  confirmed  in  the  sections  to 
follow. 
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2.10  Creating  a Suitable  Notation 
In  order  to  simplify  the  equality  constraints  repre- 
sented by  matrix  equation  (2.9.3),  the  lumped  mass  para- 
meters below  will  be  used  for  links  with  fixed  link 
length: 


m x ia 
pq  pq  pq 


(2.10.1 ) 


m y *a 

pqJpq 


pq 


m 


-2  _2  _2 


u +y 


+K 


pq  pq  pq  pq 


)*a 


pq 


(2.10.2) 

(2.10.3) 


So,  from  here  on,  the  mass  parameters  of  a link  with  p-q 
motion  specification  will  be  mpq,  x°q,  y°q,  and  k°q.  The 
original  mass  parameters  mpq,  xpq,  ypq,  and  Epq  will  be 
referred  to  as  "basic  mass  parameters"  whenever  this 
distinction  becomes  necessary. 

If  the  mass  parameters  defined  above  are  substituted 
into  Eqs.  (2.7.2),  the  result  is 


Y1pq  mpq_Xpq 

Y2pq  = ^pq 
xl  - y° 

3pq  ^pq 

Y2  - x° 

4pq  pq 


(2. 10. 4. a) 
(2.10.4-b) 
(2.10.4.C) 
(2. 10. 4. d) 


and  from  Eqs.  (2.8.2) 
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i pq 

= m -2xu  +k° 

pq  pq  pq 

(2. 10. 5. a) 

3 

2pq 

0 , 0 
= X -K 

pq  pq 

(2. 10. 5. b) 

3 

3pq 

0 

= y 

Jpq 

(2.10.5.C) 

3 

4pq 

= k° 

pq 

(2. 10.5. d) 

Another  simplication  will  be  related  to  Eq.  (2.9.1) 
which  has  been  written  in  a general  form.  Thus,  for  a 
specific  equation,  (D)1  might  represent  D*  , in  which 
case  f^CY's)  will  turn  into  f*  . For  example,  the 
equation  in  Fig.  2.10  may  be  rewritten  as 


s 


F 


2 2 
f D 
412  412 


2 

+f 

443 


2 

D 

44  3 


A third  and  final  simplification  will  be  achieved  by 
means  of  the  following  equalities: 


rpq 

= x *a 
r pq 

= y la 

rpq 

r pq 

where  "r"  is  a point  on  link  pq. 

In  the  subsequent  sections,  the  equality  constraints 
which  result  from  shaking  force,  shaking  moment,  and  input 
torque  syntheses  will  be  derived  for  seven  kinds  of 
linkage  mechanisms. 
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2.11  Equality  Constraints  for  the  Four-Bar  Linkage 
The  functions  multiplying  the  D-coef f icients  in  Fig. 
2.20  take  the  form  below  when  Eqs.  (2.10.4)  and  Eqs. 
(2.10.5)  are  substituted  into  them.  They  become  equality 
constraints  when  their  values  are  fixed  by  the  synthesis 
equation  (2.9.3). 

Shaking  Force: 
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Shaking  Moment  and  Input  Torque: 
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It  is  important  to  notice  that 
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III 

+ f4 

4 1 2 
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22  3 
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323 
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+f 

24  3 


Hence,  for  a four-bar  linkage,  one  cannot  perform  shaking 
moment  synthesis  together  with  the  synthesis  of  either  of 
the  other  two  dynamic  properties.  Fortunately,  there  is 
no  impediment  to  the  double  syntheses  of  shaking  force  and 
input  torque  for  this  type  of  mechanism. 


2.12  Equality  Constraints  for 
the  Slider-Crank  Mechanism 

The  following  functions,  which  can  also  be  considered 
as  equality  constraints,  result  from  Figure  2.21. 
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ShaKing  Moment  and  Input  Torque: 
.4  n 
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Therefore,  in  contrast  to  the  four-bar  linkage,  no  double 
syntheses  can  be  carried  out  for  the  slider-crank 
mechanism. 


2*13  Equality  Constraints  for  the  Watt  I Mechanism 
From  Figure  2.22,  the  following  functions  (or 
equality  constraints)  can  be  written. 
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If  link  12  is  the  input  link  (D4  = D4  and  D6  =0), 
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The  following  relationships  exist  between  the  pre- 


vious functions: 
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Thus,  for  the  Watt  I linkage,  no  double  syntheses  can  be 
performed,  either. 


2.14  Equality  Constraints  for  the  Watt  II  Mechanism 
From  Figure  2.23,  results 
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The  following  relationships  can  be  written  from 
above : 
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Therefore,  the  double  syntheses  of  shaking  force  and  input 
torque  can  be  performed  for  the  Watt  II  mechanism. 


2.15  Equality  Constraints  for  the 
Stephenson  I Mechanism 


The  following  equalities  result  from  Figure  2.24 
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The  above  functions  are  related  in  the  following  way 
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Hence,  no  double  syntheses  can  be  executed  for  the 
Stephenson  I mechanism. 


2.16  Equality  Constraints  for  the 
Stephenson  II  Mechanism 

The  following  functions  can  be  written  from  Figure 
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Thus,  no  double  syntheses  can  be  performed  for  the 
Stephenson  II  mechanism,  either. 
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2.17  Equality  Constraints  for  the 
Stephenson  III  Mechanism 

Figure  2.26  yields  the  following  equalities 
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4|g  = 0 if  iulz  = const.,  and  D476  = 0 if  uj76=  const. 
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From  above,  the  following  relationships  may  be 
deduced. 
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Hence,  it  is  feasible  to  carry  out  the  double  syntheses  of 
shaking  force  and  input  torque  for  the  Stephenson  III 
mechanism. 


2.18  Comments 

It  is  important  to  notice  from  the  previous  sections 
that  the  masses  of  the  links  grounded  by  a revolute  joint 
never  appeared  explicitly  in  the  equality  constraints. 

This  was  achieved  because  the  origin  of  the  local 
coordinate  system  of  each  one  of  those  links  was  attached 
to  the  center  of  the  grounded  pin  joint,  in  accordance 
with  "Rule  1"  of  section  2.2.  As  a reinforcment  to  the 
importance  of  this  fact,  the  masses  of  links  grounded  by 
pin  joints  do  not  appear  explicitly  in  the  derivations  for 
the  reaction  forces  either,  as  demonstrated  in  Appendix  D. 

All  the  sets  of  equality  constraints  obtained 
previously  contain  more  variables  than  equations,  hence 
some  of  the  variables  need  to  be  pre-specif ied  in  order  to 
solve  for  the  others.  Fortunately,  the  remaining 
variables  can  always  be  calculated  by  judiciously  solving 
the  equations  sequentially,  even  for  the  cases  of  double 
syntheses . 
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Table  2.1  presents  some  key  information,  from  the 
last  seven  sections,  relevant  to  the  cases  of  dynamic 
synthesis  that  convey  fewer  free  mass  parameters  to  the 
optimization  process  that  will  be  described  in  the  next 
chapter . 
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If  link  43  is  the  input  link;  If  link  76  is  the  input  link. 


CHAPTER  III 
DYNAMIC  OPTIMIZATION 


In  this  chapter,  the  main  aspects  of  a proposed 
interactive  software  package  for  the  dynamic  optimization 
of  linkage  mechanisms  will  be  considered.  A package  based 
on  the  contents  of  this  chapter  has  actually  been  created, 
and  will  be  described  in  the  next  chapter. 

3 • 1 Optimization  in  Mechanism  Design 
The  development  of  the  digital  computer  brought  a new 
impetus  to  the  general  field  of  optimization.  Previously 
infeasible  approaches  to  optimization  became  practical 
giving  new  directions  to  innumerous  areas  of  human  knowl- 
edge , in  particular  the  area  of  mechanism  design. 

Some  of  the  earliest  works  that  applied  the  digital 
computer  on  linkage  optimization  problems  were  performed  by 
Freudenstein  and  Sandor  [74,73]  in  1959.  They  synthesized 
path  generating  mechanisms  by  means  of  complex  number 
theory  and  an  IBM  650  computer.  Since  then,  the  number  of 
contributions  to  mechanism  design  optimization  has  in- 
creased exponentially.  Undoubtedly,  most  of  these  con- 
tributions have  emphasized  the  iterative  methods  [76-83], 
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i.s.,  methods  that  involve  "mathematical  programming" 
techniques . 

As  opposed  to  iterative  methods,  random  search 
methods  look  for  the  best  solution  by  exhaustively  con- 
sidering the  points  of  a grid  which  are  representative  of 
the  whole  design  space.  Although  the  random  search"' 
methods  are  in  general  computationally  more  expensive  than 
the  iterative  ones,  there  are  some  techniques  that  can 
increase  their  efficiency  considerably. 

An  early  use  of  the  random  search  concept  in  linkage 
optimization  occurred  in  1962  when  Roth,  Sandor  and 
Freudenstein  [ 89 J worked  on  the  synthesis  of  four-bar  path 
generating  mechanisms  with  optimum  transmissions  charact- 
eristics . 

In  1969,  Eschenbach  and  Tesar  [90]  introduced  the 
technique  of  "sequential  filtering"  by  means  of  a random 
search  method  for  the  optimal  design  of  crank  and  rocker 
four-bars.  Essentially,  a large  set  of  linkages  satisfy- 
ing four  prescribed  multiply  separated  positions  and  cer- 
tain necessary  conditions  was  generated  computationally 
and  ranked  according  to  numerous  desirable  conditions 
acting  as  sequential  filters.  Later,  Spitznagel  [91 J,  and 
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In  this  work,  the  term  "random  search"  has  a broad 
sense,  that  is,  a random  search  method  does  not 
necessarily  select  the  design  points  randomly. 
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Spitznagel  and  Tesar  j.92]  improved  this  approach  by  intro- 
ducing the  idea  of  "grid  expansion."  Through  this  techni- 
que, a region  on  the  grid  about  the  linkages  of  higher 
quality  is  selected,  then  a more  refined  grid  is  generated 
about  the  selected  region. 

The  philosophies  of  "sequential  filtering"  and  "grid 
expansion"  have  been  implemented  in  the  interactive  soft- 
ware package  that  will  be  described  in  Chapter  IV. 

3 • 2 Definition  of  the  Problem 
Given  a set  of  N linkages  with  unknown  mass  para- 
meters, possibly  coming  from  a previous  kinematic  optimi- 
zation process,  it  is  desirable  to  find  by  means  of  an 
interactive  software  package  the  best  distribution  of 
internal  masses  for  each  of  these  linkages  in  order  to 
optimize  their  dynamic  behavior.  Then,  from  the  N 
"optimum"  solutions,  select  the  linkage  that  has  the  best 
performance  with  respect  to  all  the  criteria  of  interest 
to  the  designer.  It  is  also  desirable  to  be  able  to 
balance  a pre-existing  linkage  with  known  mass  parameters. 

To  solve  such  a problem,  the  combination  of  synthesis 
and  optimization  will  be  considered  next. 

3-3  Integration  of  Synthesis  and  Optimization 
The  flowchart  in  Fig.  3.1  shows  the  primary  elements 
necessary  to  the  integration  of  dynamic  synthesis  and 
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N Promising  Linkages 
/\ 

\ 


Linkage  1,  Linkage  2 Linkage  i Linkage  N 


Figure  3.1 


Integration  of  Dynamic  Synthesis  and 
Optimization  in  Computer  Software 
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optimization  in  an  interactive  software  package  applicable 
to  the  problem  described  in  the  preceding  section. 

The  dynamic  synthesis  process  considered  in  Chapter 
II  represents  the  first  stage  of  the  whole  optimization 
method.  As  explained  in  section  2.9,  the  equality 
constraints  obtained  from  synthesis  decrease  the  number  of 
free  design  variables  which  are  carried  over  to  the  main 
optimization  stage.  This  fact  was  summarized  in  Table  2.1 
for  ten  cases  of  dynamic  synthesis.  However,  as  indicated 
in  Fig.  3.1,  synthesis  can  be  considered  successful  only 
if  the  necessary  conditions  can  be  fulfilled  and  the  over- 
all behavior  of  the  dynamic  property ( ies ) being  synthe- 
sized is  considered  satisfactory. 

In  the  remaining  part  of  this  chapter,  the  parts  of 
the  chart  in  Fig.  3.1  needing  elucidation  will  be  con- 
sidered. 


3.4  The  Necessary  Conditions 
A mechanism  undergoing  dynamic  optimization  must  have 
its  mass  parameters  subjected  to  side  constraints.  In 
this  work,  these  constraints  are  represented  by  the  in- 
equalities , 


pqmin 

< 

m 

pq 

< 

m 

pqmax 

(3.4.1 .a) 

pqmin 

< 

X 

pq 

< 

xpqmax 

(3.4.1 .b) 

pqmin 

< 

y 

pq 

< 

y 

J pqmax 

(3.4.1 .c) 
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kpqrain  - kpq  - Kpqmax  (3.4.1.d) 

where  mpq,  xpq,  ypq,  and  kpq  are  the  basic  mass  parameters 
of  a generic  link  pq  representing  all  the  moving  links  of 
such  a mechanism.  Thus,  using  the  definition  of  the 
(lumped)  mass  parameters  given  in  section  2.10,  it  is 
possible  to  transform  the  constraints  above  into  the 
following  inequalities: 


where 


m . < m < m 

pqmin  - pq  - pqmax 

x°  < x°  < 

pqmin  - pq  - Apqmax 

y°  . < y°  < y° 

k°  . < k°  < k° 

pqmin  - pq  - pqmax 


(3. 4. 2. a) 
(3 • 4 . 2. b ) 
(3.4.2.c) 
(3.4.2.d) 


x . = m x — ^ 

pqmin  pqmax  pqmin"  pq 


lf  xpqmin  1 0 


or  = m . .x  . ia 

pqmin  pqmin  pq 


o 

pqmax  ~ mpqmax *xpqmaxvapq 


or  = m ,x  ia 

pqmin  pqmax'  pq 


if  x . > 0 

pqmin  - 


if  x > 0 

pqmax  - 

if  x < 0 

pqmax  - 


(ypqmin  and  ypqmax  are  calculated  similarly) 

kpqmin  = mpqmin ^MIN ^ xpqmiJ ’ I xpqmax* ^ + 

+ [MIN ( | y . | , | y I )]2+k2  Ua2 

kpqmin1 ’ 1 kpqmax1 ' J pqmin  J pq 


82 


and 


k 


o 

pqmax  " pqrnax 


{ [MAX( | x . |,  |x 
1 pqmin ' ’ 1 


pqmax 


pqmax 


pqmax 


pq 


2 


The  inequality  constraints  (3*4.2)  for  each  moving  link  of 
a mechanism  experiencing  optimization,  together  with  the 
equality  constraints  obtained  from  dynamic  synthesis  (sec- 
tion 2.9)  represent  "necessary  conditions"  that  must  be 
satisfied.  However,  in  some  cases,  these  necessary  condi- 
tions cannot  be  satisfied  simultaneously  by  any  of  the 
design  points,  as  illustrated  in  Fig.  3.2  for  the  design 
space  of  a generic  objective  function  F = F(u1fu2). 
Therefore,  it  is  wise  to  check,  right  after  synthesis,  if 
there  are  any  points  in  the  design  space  that  fulfill  all 
the  necessary  conditions  concurrently . How  to  accomplish 
this  can  probably  be  better  explained  by  means  of  a speci- 
fic illustration.  Thus,  taking  as  an  example  the  shaking 
moment  synthesis  of  the  four-bar  in  section  2.11,  the 
following  constraints  are  obtained: 
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u2 

^Zmax  - 

Equality  Constraint 
From  Synthesis 

uZmin  - 

fl 

' • " V*  ' 

' ■ ■ ■ ' ■ 

' :■  '*  «•  - . • ,*2:  ' ' ‘ 

l 

Imin  U 

> 

(max  Uj 

a)  Only  points  on  segment  HB  satisfg  the  inequality  and 
equality  constraints  simultaneously. 


b)  There  is  no  point  in  the  design  space  satisfying  the 
equality  and  inequality  constraints  simultaneously. 


figure  3*2  Two  Examples  of  Possible  Necessary  Conditions 
for  a Generic  Objective  Function  F=F(u1 ,u2) 
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It  should  be  noticed  that  all  the  f's  at  the  left-hand- 
side  of  the  equations  above  are  known  numbers.  These 
equations  and  the  inequality  constraints  for  the  three 
moving  links  can  be  satisfied  simultaneously  only  if  the 
following  inequalities  are  true: 

o o 

^i2min  - ^312  - ^i2max 

O a o 4 

x . +m  . -2x  +ku  < f < 

12mm  23mm  23max  23min  - 412  - 
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o 

2 3max 


< f 


223 


< X 


o 

2 3max 


< f 


343 


4 3max 


4 


o 

4 3max 


< x 


2 4 3 


Essentially,  the  same  checking  procedure  applies  to  other 
linkages  and/or  dynamic  synthesis  cases. 
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5 • 5 The  Desirable  Conditions 
When  designing  a mechanism,  the  designer  always 
"desires"  to  optimize  certain  important  characteristics  of 
the  mechanism.  To  accomplish  this,  he  must  have  some 
criteria  by  which  the  quality  of  such  attributes  can  be 
measured.  These  criteria,  which  will  be  referred  to  here 
as  "desirable  conditions,"  are  generally  in  conflict. 
Hence,  the  designer  must  look  for  a solution  that  reflects 
the  best  compromise  for  his  dilemma.  This  can  be  achieved 
by  means  of  an  objective  function  composed  of  weighted 
desirable  conditions,  as  described  in  the  next  section. 

There  are  various  well-known  desirable  conditions 
that  can  be  used  by  a designer,  but  ideally,  he  should 
have  the  freedom  to  customize  his  own.  In  any  case,  a 
desirable  condition  should  always  fall  within  the 
following  broad  classification: 

Glass  I - conditions  dependent  only  on  the  geometry 
of  the  mechanism. 

Subclass  la  - conditions  of  Class  I with  no 

dynamic  implications,  e.g.,  the 
location  of  a fixed  pivot. 

Subclass  lb  - conditions  of  Class  I with  dynamic 
implications,  e.g.,  the  max.  value 
of  the  g-function  of  an  output 
link. 

Class  II  - conditions  dependent  on  the  mass  parameters 
of  the  mechanism,  e.g.,  maximum  value  of 
the  shaking  force. 
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A complete  package  for  the  computer  aided  design 
(CAD)  of  mechanisms  must  include  a module  for  kinematic 
optimization,  followed  by  another  for  dynamic  optimiza- 
tion. The  former  should  have  its  objective  function  based 
on  conditions  from  Class  I,  and  the  latter  on  conditions 
from  Class  II.  As  indicated  in  the  rectangle  at  the 
bottom  of  Fig.  3.1,  there  must  be  a final  comparative 
stage  at  the  very  end  of  the  whole  optimization  process. 

In  that  stage,  only  conditions  from  Subclass  la  and  Class 
II  need  to  be  considered. 

It  is  opportune  to  emphasize  that  the  conditions  that 
belong  to  Subclass  lb  represent  dynamic  "predictors,"  that 
is,  they  forecast  the  dynamics  of  a mechanism  when  only 
the  geometry  is  known.  Such  predictors  have  been  success- 
fully applied  to  cam  systems  by  Tesar  and  Matthew  [51], 
and  to  Geneva  mechanisms  by  Taat  and  Tesar  [50] . Although 
they  can  act  as  powerful  desirable  conditions  in  a pre- 
vious kinematic  linkage  optimization  package,  like  the  one 
described  in  [93],  they  do  not  possess  the  same  value  for 
the  dynamic  optimization  package  considered  in  this 
work.  Here,  it  is  preferable  to  address  the  dynamic  pro- 
perties of  a linkage  directly  by  means  of  the  following 
desirable  conditions  from  Class  II: 
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a) 

b) 

c) 

d) 

e) 

f) 


'l  s t 

Maximum  magnitude  of  total  shaking  force;  aF 

° ’ max 


Max.  magnitude  of  total  bearing  reactions;  F 


Imax  ’ 


2 max 


, etc . 


s,„t 


Peak-to-peak  value  of  total  shaking  moment;  M 


ptp 


Root-mean-square  value  of  total  shaking  moment; 

sMt 

rms 


Maximum  absolute  value  of  total  shaking  moment; 


V 

max 


Lt 


Peak-to-peak  value  of  total  input  torque;  T 


ptp 


g)  Root-mean-square  value  of  total  input  torque; 

iTt 

rms 

h)  Maximum  absolute  value  of  total  input  torque; 

V 

max 

The  objective  function  in  the  software  package  to  be 
described  in  the  next  chapter  must  be  established  from  the 
desirable  conditions  above. 


3 • 6 The  Objective  Function 
Given  a set  of  k conflicting  desirable  conditions 
(from  the  list  presented  in  the  previous  section),  it 


The  adjective  "total"  in  this  list  denotes  summation  of 
the  contributions  from  internal  masses  and  external 
loads . 
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becomes  necessary  to  create  the  following  objective 
function 


S = z uiSi  ; i = 1 , 2,  3,  ...  , k (3.6.1 ) 


where 

S = cumulative  score 

i = desirable  condition  counter 

si  = score  signifying  performance  of  the  linkage 
relative  to  the  i™  condition 

ui  = weighting  factor  signifying  importance  of  i^ 
condition  relative  to  the  other  i-1  desirable 
conditions. 

The  k weighting  factors  are  defined  by  the  designer 
according  to  his  experience  and  intuition. 

The  score  is  calculated  as  follows: 


where 


(3.6.2) 


wi  = worst  value  for  the  ith  condition 

= best  value  for  the  i^  condition 

Ai  = actual  value  for  the  ith  condition  (Bi<Ai<Wi) 

ni  = exponent  allowing  scores  to  be  assigned  in  a 
non-linear  manner  (nj_>0) 


Since 


y 


in  this  work,  = 0 always,  Eq.  (3.6.2)  becomes 
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; 0 < A.  < W.  (3.6.3) 

According  to  Eq.  (3.6.3)  and  its  graphical  representation 
in  Fig.  3.3.,  the  lower  the  actual  value  Ai , the  more  the 
score  Si  approaches  1,  its  highest  possible  value.  The 
exponent  n^  must  also  be  set  by  the  designer  in  conformity 
with  his  experience  and  intuition. 

It  is  desirable  to  normalize  the  cumulative  score  S, 
so  that  its  highest  possible  value  is  always  the  same, 
under  the  unrealistic  but  appropriate  assumption  that  all 
the  S^'s  can  simultaneously  attain  the  value  1.  In  this 
work,  such  referential  maximum  value  has  been  chosen  to  be 
100.  Thus,  from  Eq.  (3.6.1),  it  is  possible  to  write 

+ u2  + ...  + pk  = 100  (3.6.4) 

If  u ' i > w 1 £ > • • • » u'k  are  values  chosen  by  the 

designer  for  the  weighting  factors  of  all  the  k desirable 
conditions,  such  that 

0 < u?  < 10  ; i = 1,  2,  ...  , k 

then 

ui  = if  y;  t o ; i = 2,  3,  ...  , k (3.6.5) 
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Figure  3.3  Some  Possible  Scoring  Curves  for  Desirable 
Condition  i 
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Substituting  Eqs.  (3.6.5)  into  Eq.  (3.6.4)  yields 

u(u'+u'+  ...  + u')*y  ' = 100 

112  K 1 


or 


u ^ = IOOu^tZu?  ; i = 1 to  k 


(3.6.6) 


Extending  this  result  to  all  the  other  weighting  factors, 
gives 


V-  = IOOuItEu'  ; i = 1 to  k 
J J 1 

j = 1 » 2,  ...  , k 


(3.6.7) 


CHAPTER  IV 

AN  INTERACTIVE  CAD  PACKAGE 


The  theory  presented  in  Chapters  II  and  III  has  been 
implemented  in  an  interactive  software  package,  which  is 
currently  applicable  only  to  four-bar  mechanisms. 

However,  it  can  be  expanded  to  encompass  other  linkages, 
with  relative  ease.  This  chapter  is  not  intended  to  be  a 
user's  manual.  In  fact,  the  package  is  intuitive  and 
self-explanatory,  but  it  is  recommended  that  the  user  be 
aware  of  the  flowchart  presented  in  the  preceding  chapter. 

4.1  General  Description  of  the  CAD  Package 
The  software  package  addressed  by  this  chapter  has 
been  written  in  APL  which  is  a computer  language  very 
suitable  for  manipulation  of  arrays.  The  APL  functions 
that  compose  the  package  can  be  classified  into  two 
groups,  namely,  the  "Graphics  System"  and  the  "Application 
Program."  The  Graphics  System  consists  of  the  primitive 
graphic  functions  used  to  create  graphics  elements  such  as 
segments  of  lines,  polygons,  circles,  and  character 
strings  on  the  screen  of  a display  terminal.  These 
i unctions , which  are  listed  in  Appendix  E,  are  very  much 
dependent  on  the  low-level  architecture  of  the  terminal 
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and  tha  X-Y  coordinate  system  of  the  physical  screen. 

They  are  compatible  with  the  Tektronix  4010/4012/4013/4014 
and  4015  graphics  terminals  (Fig.  4.1),  or  any  other 
terminal  capable  of  simulating  a Tektronix  from  the  4010- 
series  mentioned  above.  The  Application  Program  consists 
of  the  functions  which  directly  implement  the  theory 
discussed  in  Chapters  II  and  III.  These  functions  are 
listed  in  Appendix  F. 


4 . 2 The  Application  Program 

The  Application  Program  is  based  on  the  flowchart  in 
Fig.  3.1  which  has  been  reproduced  for  convenience  in  this 
chapter  as  Fig.  4.2.  This  chart  is  not  complete  in 
details.  It  has  been  created  to  illustrate  the  main 
logical  decision  making  stages  involved  in  the  package. 

Figure  4.3  is  a photograph  of  part  of  the  screen  of  a 
Tektronix  4015  terminal.  It  shows  the  initial  stage  of 
the  computer  aided  design  process  where  information  about 
the  geometry  and  motion  of  the  four-bar  is  entered  into 
the  computer.  This  photograph  does  not  show  the  image  of 
a referential  four-bar  linkage  drawn  at  the  lower  right 
corner  of  the  screen.  Such  image  was  similar  to  the  one 
shown  in  Fig.  4.4.  An  important  feature  of  the 
Application  Program  is  that  it  allows  the  user  to  choose 
the  locations  for  the  x-y  coordinate  systems  of  the  moving 
links.  The  user  may,  for  example,  choose  the  local 
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Figure  4.1 


Tektronix  4015  Computer  Display  Terminal 


95 


N Promising  Linkages 
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' \ 


Linkage  1.  Linkage  2 Linkage  I Linkage  N 


Figure  4.2  Integration  of  Dynamic  Synthesis  and 
Optimization  in  the  CAD  Package 
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Figure  4-3  Initial  Stage  of  the  CAD  Package 
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Figure  4-4  A Four-Bar  Linkage  and  a Possible  Set  of  Local  Coordinate  Systems 
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systems  shown  in  Fig.  4.4.  For  all  subsequent  input  and 
output  operations,  the  Application  Program  performs  all 
the  necessary  transformations  between  the  user's 
coordinate  systems  and  the  set  of  x-y  coordinate  systems 
which  is  internal  to  the  program  and  completely 
transparent  to  the  user. 

The  rest  of  this  chapter  will  address  some 
peculiarities  of  the  application  program  with  respect  to 
external  loads,  synthesis,  optimization,  and  the  basic 
mass  parameters. 


4 • 3 External  Loads 

The  Application  Program  offers  the  user  two  ways  for 
entering  information  about  external  loads.  In  one  of 
them,  the  numerical  data  for  the  external  loads  must  be 
typed  in  one  at  a time  for  each  position  of  the  input 
crank.  However,  if  the  data  to  be  entered  assume  a 
constant  value  for  certain  interval  of  motion,  this  value 
needs  to  be  entered  only  once.  The  other  way  to  input  the 
external  loads  into  the  computer  is  by  means  of  a 
"calculator"  portrayed  on  the  screen,  and  the  crosshair 
cursor,  as  shown  in  Fig.  4.5.  The  cursor  can  be 
positioned  by  two  thumbwheels  at  the  far  right  of  the 
keyboard  (Fig.  4.1),  or  optional  joystick.  This 
calculator  uses  the  RPN  logic  and  has  a stack  with  four 
storage  registers  (X,  Y,  Z and  T).  Every  time  a number  is 
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Figure  4.5  A Computer  Generated  Calculator 
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entered  into  tne  X register,  or  an  operation  is  performed, 
a new  representa tion  of  the  stack  appears  near  the 
previous  ones.  Figure  4.5  shows  six  such  representations , 
where  IA,  N1 , and  R1  signify  input  angle,  first  number 
entered  through  the  keyboard  of  the  terminal,  and  the 
result  of  first  operation  performed  by  the  calculator, 
respectively.  The  calculator  is  useful  when  there  is  a 
known  function  of  the  input  angle  from  which  the  data  to 
be  entered  can  be  calculated  for  a specified  interval  of 
motion.  For  example,  Fig.  4.6  shows  the  plot  of  a 
fictitious  load  versus  input  angle  created  by  means  of  the 
calculator.  It  is  important  to  understand  that  the 
computer  does  not  retain  the  analytical  expressions  for 
each  of  the  "function  segments,"  but  instead  it  keeps  only 
the  function  value  for  each  of  the  equally  spaced 
positions  of  the  input  crank. 

Either  the  global  coordinate  system  or  the 
appropriate  local  coordinate  system  can  be  used  as  a 
reference  for  an  external  force.  In  addition,  the  point 
of  application  of  a force  on  a moving  link  can  be 
continuously  changing. 

When  the  static  analysis  is  finished,  the  application 
program  asks  the  user  if  he  wants  to  see  the  plot  of  the 
shaking  force,  shaking  moment,  input  torque  or  any  of  the 
bearing  reactions,  due  to  the  external  loads. 


3. OOFOO 
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Figure  4.6  A Fictitious  Load  Entered  by  Means  of  the  Calculator 
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4 • 4 Dynamic  Synthesis 

Even  though  it  is  not  explicitly  indicated  by  the 
flowchart  in  Fig.  4.2,  there  is  a stage  in  the  Application 
Program,  before  dynamic  synthesis,  where  the  user  is 
allowed  to  enter  mass  parameters  for  the  four  bar,  in  case 
he  wishes  to  do  so.  This  is  useful  if  the  user  wants  to 
balance  an  existent  linkage  or  if  he  has  a good  feeling 
about  what  the  optimum  mass  parameter  might  be. 

As  demonstrated  in  Chapter  II,  either  shaking  moment 
synthesis  or  the  double  synthesis  of  shaking  force  and 
input  torque  is  possible  for  a four-bar  linkage.  When  the 
synthesis  stage  is  reached,  the  Application  Program 
displays  the  curve  of  the  dynamic  property  to  be 
synthesized  that  results  from  the  external  loads  and/or 
the  internal  masses  of  the  four-bar. 

Figure  4.7  shows  a shaking  moment  plot  and  seven 
precision  points  selected  with  the  crosshair  cursor,  and 
Fig.  4.8  shows  the  synthesized  shaking  moment  curve 
passing  through  all  the  seven  precision  points.  Synthesis 
can  be  performed  as  many  times  as  the  user  wants  until  a 
satisfactory  curve  is  obtained.  If  more  than  seven 
precision  points  are  selected  for  shaking  moment 
synthesis,  and  more  than  three  for  input  torque  synthesis, 
they  are  approximated  in  the  least  squares  sense.  Figure 
4.9  shows  a shaking  force  hodograph  and  two  precision 
points,  and  Fig.  4.10  shows  the  synthesized  shaking  force. 
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Figure  4.7  A Shaking  Moment  Plot  and  Seven  Precision  Points 
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Figure  4.8  A Synthesized  Shaking  Moment 
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Figure  4-9  A Shaking  Force  Hodograph  and  Two  Precision  Points 
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Figure  4.10  A Synthesized  Shaking  Force 


107 


4 . 5 Optimization 

Before  the  optimization  stage  can  be  initiated,  the 
user  has  to  define  the  objective  function.  Therefore,  for 
each  of  the  desirable  conditions  considered  (from  the  list 
in  section  3*5),  the  user  must  enter  the  worst  value, 
weighting  factor,  and  the  exponent  of  Eq.  (3.6.3).  To 
select  this  exponent,  the  user  indicates  a trial  value, 
and  the  Application  Program  plots  the  corresponding  score 
curve  on  the  screen  (Fig.  4.11).  This  process  can  be 
repeated  as  many  times  as  necessary  until  a value  that 
meets  the  user's  intuition  is  found. 

As  mentioned  in  the  last  chapter,  optimization  in  the 
package  is  performed  by  means  of  a random  search  method 
improved  by  the  techniques  of  "sequential  filtering"  and 
"grid  expansion."  Thus,  for  every  design  point  (m23,  x°3) 
belonging  to  a specified  grid  in  the  region  of  the  design 
space  defined  by  the  constraints 


[ 

2 3 nun 

< m 

2 3 

< m 

2 3max 

(4.5.1a) 

o 

2 3niin 

x o 
< X 

2 3 

/ 0 

< x , 

2 3 max  ’ 

(4.5.1b) 

the  equality  constraints  obtained  from  synthesis  are 
solved  sequentially,  and  each  calculated  mass  parameter  is 
in  turn  compared  to  its  specified  limits.  If  the  design 
point  is  not  rejected  during  this  process,  the  Application 
Program  proceeds  to  calculate  the  desirable  conditions 
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Figure  4.11  Two  Score  Curves  (ni  = 0.6  and  2) 
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also  sequentially,  however,  before  the  calculation  of 
condition  i+1 , the  condition  i is  checked  against  its 
permissible  worst  value.  If  the  constraint  imposed  by 
this  extreme  value  is  not  satisfied,  the  design  point  is 
rejected  and  a new  one  selected.  The  "grid  expansion" 
technique  (section  3.1)  allows  the  user  to  specify  a two- 
dimensional  grid  in  the  region  of  the  design  space  defined 
by  the  inequalities  (4.5.1),  as  indicated  in  Fig.  4.12.  A 
second  and  more  refined  grid  can  be  generated  about  the 
best  point  from  the  first  grid.  This  process  can  be 
repeated  as  many  times  as  desirable.  At  the  end  of  the 
optimization  stage,  the  optimum  mass  parameters  found  are 
displayed  on  the  screen,  as  shown  in  Fig.  4.13. 

4.6  The  Transformation  from  Lumped  to 
Basic  Mass  Parameters 

The  lumped  mass  parameters  defined  in  section  2.10 
greatly  simplify  the  formulation  and  solution  of  the 
synthesis  and  optimization  problems.  In  fact,  for 
mechanisms  more  complex  than  the  four-bar,  these  problems 
would  become  extremely  complicated  without  the 
introduction  of  such  parameters.  However,  since  the  basic 
mass  parameters  are  the  ones  that  ultimately  interest  the 
designer,  the  inverse  transf ormation  from  the  "lumped 
space"  to  the  "basic  space"  is  necessary.  In  the  package, 
this  transf orma tion  is  performed  during  the  optimization 
phase  for  every  point  of  the  lumped  design  space  that 
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Figure  4.12  The  Specification  of  a Grid  in  the  Design  Space 
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Figure  4.13  The  End  of  the  Optimization  Stage 
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satisfies  the  necessary  conditions  imposed  on  the  lumped 
mass  parameters.  The  reason  why  the  transformation  is  not 
performed  at  the  end  of  optimization,  for  the  optimum 
solution  only,  is  because  the  transformed  point  may  not 
satisfy  all  the  side  constraints  on  the  basic  mass 
parameters,  as  exemplified  in  Fig.  4.14.  The  following 
algorithm  has  been  included  in  the  Application  Program  to 


perform 

the  transformation  above: 

Link  23 

( the  coupler ) : ^ 

i) 

Calculate  x , y , and  k from  known  m , 

2323  23  23 

0 0 , , 0 
x , y , and  k 
2 3 2 3 2 3 

ii) 

Proceed  if  side  constraints  for  link  23  are 
satisfied,  otherwise  reject  point. 

Link  12 

( the  input  link) : 

i) 

Find  the  interval  Am  from  known  x°  , as  illu- 

1 2*  12 
strated  in  Fig.  4.15. 

Note:  Am  = [m  . ,m  ] if  x°  = 0. 

12X  palin’  i2max  12 

ii) 

Similarly,  find  Am  from  known  y°  . 

i2y  J 12 

1 


The  reference  for  this  algorithm  is  the  four-bar  drawn 
on  page  47* 
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m23  x23  “ *23  x 23m  in  = n>23min  x23min 


nOTE : xjjitnou,,,  satisfies  the  side  constraint 

x23min  i x23  i X23max 

Bul  x23iound  does  not  satisfy  the  constraint 

x23min  i x23  S X23max 


Figure  4.14 


An  Illustrative  Example  of  the  Calculation  of 
x23  from  known  m23  and  x°3 
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X|2  > 0 


Figure  4.15  Two  Examples  of  the  Calculation  of  Am 


12  X 
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iii ) Find  Am  = Am  f|  A 

12  1 2X  i2y 

If  Am  = c+>  (empty  set),  reject  point. 


Link  43: 

i)  Find  Am  from  known  x°  by  the  same  method 

43X  43 

used  for  link  12. 


ii)  Similarly,  find  Am  from  known  y°  . 

4 3Y  4 3 


iii)  Find  Am'  = Am  A Am 

4 3 4 3X  4 3y 

If  Am'  = <j>,  reject  point. 


iv)  Find  Am  , from  known  x°  , y°  , and  k°  , COn- 
43*  43  43  43 

considering  the  following 


o /_2  _2  -2  2 

k = m (x  +y  +k  ) ^a 
43  43  43  43  43  43 


(4.6.1 ) 


Substituting  x2  = [(a  x°  )*m  ]2  and 

43  43  43  43 


2 2 

y = [(a  y°  )im  ] into  Eq.  (4.6.1),  gives 
43  4343  43  & 


2 _ ~ 1 


- 2 


a k = a [(x  ) +(y  ) Jm  +m  k (4.6.2) 

43  43  43  43  J 4 3 43  43  43 


“ 1 - 2 

C,  = C m +m  k 
K xy  43  4343 


(4.6.3) 


where 


C 


k 


2 q 

a k >0  and  C 
4343  xy 


2-(y°3)2]  > 0 


If  C = 0,  find  Am  , by  the  same  method  used  to  find 
xy  4 3k 

Am  and  Am  . Otherwise  refer  to  Fig.  4.16. 

4 3 x 4 3 y 
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and 
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c*2 

Ul43p«afc  c r 

R43p®ak  3 ^ 

= ^ ^4}min 

Figure  4.16  The  Calculation  of  Am  . 
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Note : 


- 2 - 2 
If  k , < k . , 

43peak  43min’ 


reject  point  (Fig. 


4.15). 


v)  Find  Am  = Am ' n Am  , . 

43  43  43k 

If  Am  = a,  reject  point;  otherwise  it  can  be 
43 

accepted. 


The  algorithm  just  presented  can  be  applied  to  any 
linkage  since  each  moving  link  of  a given  mechanism  will 
always  correspond  to  one  of  the  three  types  of  links 
considered  above. 


CHAPTER  V 
EXAMPLES 


In  this  chapter,  the  redistribution  of  the  internal 
masses  of  a fully  force  balanced  four-bar  linkage  acted  by 
an  external  torque  will  be  considered  to  illustrate  the 
two  basic  ways  in  which  the  CAD  package  described  in 
Chapter  IV  can  be  applied.  This  will  be  performed  by 
means  of  one  example  involving  the  double  syntheses  of  the 
input  torque  and  shaking  force  of  the  fully  force  balanced 
four-bar,  and  another  example  addressing  the  shaking 
moment  synthesis  of  the  same  mechanism. 

5 • 1 The  Known  Linkage  and  Load 
The  mass  parameters  of  the  mechanism  in  Fig.  5.1  make 
the  first  four  functions  given  on  page  45  vanish,  that  is, 

Re (f2  ) = Im( f 2 ) = Re  (f2  ) = Im(f2  ) = 0 

4 12  412  443  443 

Thus,  the  inertial  shaking  force  acting  on  the  foundation 

of  the  mechanism  is  zero,  and  so  is  the  shaking  force  due 

to  the  external  torque  T . This  torque  has  the  following 

4 3 

specification : 
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Parameter 

Link 

12 

23 

13 

11 

length  [mm] 

25 

102 

75 

76 

mass  [Kg] 

0.150 

0.103 

0.300 

— 

x [mm] 

-8.9 

51 

-13 

- — 

g [mm] 

0 

0 

0 

— 

k [mm] 

20 

51 

35 

— 

Ilote  : The  shaking  moment  reference  point  for  this 
chapter  is  (U.V)  = (0.D3S  m , D) 


Figure  5.1  A Fully  Force  Balanced  Four-Bar  Linkage 
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t43  = 0 ; *12  = 0°  to  35° 

T = -70cos[(<t>  -0.6109)  *0.5694 1+70  Nra;  <t>  = 55°  to  240° 

j 12  12 

T = 0 ; <(>  = 240°  to  360° 

4 3 12 

where  <t>12,  in  the  argument  of  the  cosine  function  above, 
must  be  in  radians.  The  graphical  representation  of  T43 
is  presented  in  Fig.  5.2. 

Figures  5.3  and  5.4  show  the  total  shaking  moment  and 
the  total  input  torque  of  the  mechanism  described  in  Fig. 
5.1.  The  maximum  values  of  the  bearing  reactions  are 


F _ f 

imax  i+max 


5390  N 


F 

2max 


5720  N 


F 

3 max 


4870  N 


5 . 2 Example  1 : Input  Torque  and  Shaking  Force  Synthesis 

The  following  assumptions  apply  to  this  example: 

1.  The  foundation  of  the  mechanism  is  massive,  and 
the  shaking  force  is  not  a very  important 
property. 

2.  It  is  desirable  to  reduce  the  root-mean-square 
value,  the  peak-to-peak  value,  and  the  maximum 
absolute  value  of  the  input  torque. 
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Figure  5.2  External  Torque  T43  in  Nra 


PRESS  ‘RETURN'  TO  CONTINUE 
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Figure  5.3  Total  Shaking  Moment  in  Nm  of  the  Fully  Force  Balanced  Four-Bar 


PRESS  'RETURN'  TD  CONTINUE 


123 


in 


Figure  5.4  Total  Input  Torque  (in  Nm)  of  the  Fully  Force  Balanced  Four-Bar 
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Figures  5.5  and  5.6  show  the  input  torque  and  shaking 

force  synthesized  by  means  of  the  CAD  package  described  in 

the  last  chapter.  In  accordance  with  Chapter  II,  the 

number  of  free  design  variables  was  reduced  from  9 to  2, 

by  the  synthesis  process.  The  mass  parameters  that  were 

carried  over  to  the  optimization  phase  were  m and  x°  . 

2 3 2 3 

The  objective  function  given  by  Eqs.  (5.6.1)  and  (5.6.3) 
was  specified  as  indicated  in  Table  5.1. 


Table  5.1 

Specification  of  Objective  Function 
for  Example  1 


* 

Desirable  Condition 

f 

yi 

W. 

l 

n . 

l 

Fjmax  ; d=1>*-->4 

10 

5000  N 

1 

V, 

ptp 

6 

350  Nm 

1 

V 

rms 

6 

100  Nm 

1 

sMt 

max 

6 

190  Nm 

1 

The  optimum  solution  obtained  is  presented  in  Table 
5.2.  The  corresponding  shaking  moment  plot  is  shown  in 
Fig.  5.7,  and  the  maximum  values  of  the  bearing  reactions 
are 


F 

imax 


= 2990  N 


F 


2 max 


2840  N 
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Figure  5.5  Total  Synthesized  Input  Torque  (in  Nra)  Obtained  in  Example 
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Figure  5.6  Synthesized  Shaking  Force  (in  N)  Obtained  in  Example 


PRESS  'ftETURN'  TO  CONTINUE 
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Figure  5.7  Total  Shaking  Moment  (in  Nm)  for  the  Optimum  Solution  in  Example 
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F = F 
3inax  4max 


2210  N 


Table  5.2 

Optimum  Solution  for  Example  1 
(Score  = 40) 


Parameter 

Link  12 

23 

43 

mass  [Kg] 

0.045 

0.0787 

0 . 066 

x [mm] 

12.7 

50.8 

38.1 

y [mm] 

2.1 

0.5 

1.6 

k [mm J 

— 

34.9 

27.9 

5.3 

Example  2:  Shaking  Moment  Synthesis 

If  the  foundation  of  the  four-bar  in  Fig.  5.1  is  not 
massive,  the  shaking  moment  (Fig.  5.5)  may  cause  too  much 
vibration  on  the  surroundings  of  the  mechanism.  The 
synthesized  shaking  moment  shown  in  Fig.  5.8  represents  an 
attempt  to  reduce  these  undesirable  vibrations.  As  in 
example  1 , the  number  of  free  design  variables  was  also 
reduced  to  2 by  this  synthesis  process.  The  objective 
function  for  the  optimization  stage  was  specified  as 
indicated  in  Table  5.3. 

The  optimum  solution  obtained  is  presented  in  Table 
5.4.  The  corresponding  input  torque  and  shaking  force 
plots  are  shown  in  Figs.  5*9  and  5.10,  and  the  maximum 
values  of  the  bearing  reactions  are 


PRESS  RETURN 
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Figure  5.8  Synthesized  Shaking  Moment  (in  Nm)  Obtained  in  Example 
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Table  5.3 

Specification  of  Objective  Function 
for  Example  2 


* 

Desirable  Condi tion 

1 

wi 

. 

l 

n . 

l 

Fimax’ 

Ft 

i+max 

6 

5000  N 

1 

Ft 

2 max  ’ 

Ft 

3max 

5 

5000  N 

1 

sFt 

max 

10 

4000  N 

2 

V. 

p tp 

8 

130  Nm 

1 

ITt 

rms 

8 

35  Nm 

1 

I t 
max 

8 

80  Nm 

1 

* Section  3.5. 

Table 

Optimum  Solution 
(Score  = 

5.4 

for  Example  2 
49.3) 

Parameter 

Link  12 

23 

43 

mass  [Kg] 

0.062 

0.0647 

0.051 

x [mm] 


11 .2 


1.0  17.2 


y [mm] 
k [ram] 


3.3 


2.1  10.9 


36.2 


32.7 


PRESS  ‘RETURN’  TO  CONTINUE 
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Figure  5.9  Total  Input  Torque  (in  Nm)  for  the  Optimum  Solution  in  Example 
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Figure  5.10  Shaking  Force  (in  N)  for  the  Optimum  Solution  in  Example 
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F 


2380  N 


imax 


F 


2230  N 


2max 


F 


2080  N 


3 max 


F 


2040  N 


4max 


5.4  Comments 


The  previous  examples  have  been  included  in  this  work 
to  illustrate  the  applicability  of  the  CAD  package 
described  in  the  last  chapter.  Clearly,  the  solution 
obtained  in  example  2 is  better  than  the  one  achieved  in 
example  1.  However,  improved  solutions  could  very 
probably  be  obtained,  in  both  examples,  if  more  time  had 
been  dedicated  to  the  synthesis  and  optimization 
processes.  This  was  not  possible  due  to  high  instability 
of  the  digital  computer  used. 


CHAPTER  VI 

CONCLUSIONS  AND  RECOMMENDATIONS 
6.1  Conclusions 

As  indicated  in  the  introductory  chapter  of  this 
work,  mechanisms  and  machines  play  a very  important  role 
in  modern  life.  Despite  the  existence  of  innumerable 
sophisticated  mechanical  systems  nowadays,  there  is  still 
a need  for  powerful  and  efficient  means  to  design  faster 
and  more  precise  machines.  It  is  the  belief  of  the  author 
that  this  dissertation  represents  a substantial  contribu- 
tion to  the  fulfillment  of  such  need.  The  theory 
contained  in  Chapter  II  is  a profound  rearrangement  of  the 
work  done  by  Elliott  in  Ref.  [1],  The  concepts  of  point 
masses  and  composite  point  masses  greatly  facilitate  the 
derivation  of  the  expressions  for  the  shaking  force, 
shaking  moment,  and  input  torque  of  any  planar  linkage. 
Thus,  by  just  looking  at  the  diagram  of  a mechanism,  one 
can  immediately  write  down  these  expressions  in  a form 
free  of  linear  dependencies  and  hence  suitable  for  dynamic 
synthesis.  The  lumped  mass  parameters  defined  in  section 
2.10  considerably  simplify  the  formulation  and  solution  of 
the  synthesis  and  optimization  problems.  The  optimization 
technique  presented  in  Chapter  III  allows  the  designer  to 
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impose  inequality  constraints  on  the  mass  parameters  so 
that  the  chance  of  obtention  of  an  impractical  solution  is 
highly  reduced.  The  CAD  package  described  in  Chapter  IV 
confirmed  the  validity  and  usefulness  of  the  theory 
contained  in  Chapters  II  and  III. 

6 . 2 Restrictions  and  Recommenda tions  for  Future  Works 
The  following  items  are  limitations  of  this  work  and 
recommendations  for  future  works: 

1.  A future  research  to  extend  the  method  contained 
here  to  systems  composed  of  several  inter- 
connected planar  linkages  is  highly  desirable. 

An  extension  to  encompass  spatial  linkages  would 
also  be  valuable. 

2.  A sound  combination  of  dynamic  synthesis  and 
spring  synthesis  would  be  very  useful,  specially 
for  balancing  mechanisms  that  operate  under  large 
loads  that  cannot  be  compensated  by  the  internal 
masses  only. 

3.  The  CAD  package  presented  here  is  applicable  only 
to  four-bar  mechanisms.  Future  software  packages 
should  consider  other  linkages  besides  the  four- 
bar. 

4.  The  "calculator"  implemented  in  the  CAD  package 
described  here  cannot  be  easily  used  to  enter  a 
load  which  is  a function  of,  for  example,  an 
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output  link  angle.  In  addition,  the  package  does 
not  permit  the  user  to  enter  an  external  load  by- 
means  of  a graphics  tablet,  if  the  plot  of  such 
load  is  available. 

5.  A profound  investigation  to  improve  the  effi- 
ciency of  the  optimization  phase  is  recommended. 

6.  The  use  of  a rapid  raster  display  terminal 
instead  of  a storage  tube  one  is  also  recom- 
mended, since  the  former  is  generally  more 
suitable  for  dynamic  picture  manipulation  and 
selective  erasure  of  just  some  parts  of  the 
picture  on  the  screen. 

7.  An  important  part  of  the  design  of  a mechanism 
consists  of  the  definition  of  the  complete  shape 
of  its  moving  links,  once  the  optimum  mass 
parameters  are  known,  and  the  verification  for 
strength  of  these  links.  A complete  CAD  system 
for  mechanism  design  should  offer  powerful 
interaction  capabilities  during  this  crucial  part 
of  the  design  process.  Additionally,  the  real- 
time animation  of  the  solution  mechanism  would 
certainly  be  appreciated  by  the  designer. 

8.  Finally,  the  designer  should  have  some 
flexibility  to  customize  his  own  desirable 
conditions  to  be  used  by  the  optimization 
procedure . 


APPENDIX  A 

DERIVATION  OF  INERTIAL  INPUT  TORQUE 
BY  MEANS  OF  KINEMATIC  INFLUENCE  COEFFICIENTS 


The  concept  of  velocity  and  acceleration  influence 
coefficients  can  be  used  profitably  for  kinematic  and 
dynamic  analysis  of  multi-degree  of  freedom  mechanisms 
[21,32,40,45].  Here,  this  concept  will  be  used  for  the 
derivation  of  the  equivalent  inertia  (I*q)  at  the  input 
shaft  of  one-degree  of  freedom  linkages  with  constant 
input  angular  velocity.  Once  Ip^  is  obtained  in  a proper 
form,  the  expression  for  the  inertial  input  torque  can 
easily  be  deduced. 

If  the  formulation  derived  in  Fig.  A.1  for  the 
generic  point  "r"  is  applied  to  the  center  of  gravity,  the 
kinetic  energy  of  link  pq  can  be  written  as 


KE 


pq  * 0-5(mpq[(8pq)2  + («pq)L«pq^q(g*q)2}4  = 0.5  I*  * 


:Pq  = ”pq[(ipq)2+(ino)2]+mn„k2n(g+„)2 


pq ' 


pq  pq  pq 


pq" 

(A.1  ) 


Substituting 


gU 


= gU-g*  ( 


pq 


x s _+yr 


P Pq  Pq  pq  pq  pq 


) 


and 
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figure  A.1  The  Concept  of  Velocity  Influence  Coefficients 
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-v 


°pq  gp+gpq^XpqCpq-ypqSpq^ 
into  Eq.  (A.1 ) , gives 

Ipq  = VL(gp)2-2gpgpq(Wpq+Wpq)  + (Sp)2  + 


+ 2SPSPq^XpqCpq"ypqSpq^+mpq^Xpq+^pq+^pq^Spq^ 


(A. 2) 

Now,  considering  point  "q,"  the  following  can  be  written 


Similarly, 


g„  = gp-gjq(apqspq) 


u 

!q 


< f>  ( u Ux 

gpqSpq  “ (gp-gq)tapq 


V V <b  , x 

Sq  “ Sp+Spq  apqGpq 

„<(>  „ / V Vs 

gpq  pq  " Sq_Sp  :apq 


(a. 3) 


(A. 4) 


Squaring  Eqs.  (A. 3)  and  (A. 4),  and  combining  them,  yields 


&pq  LkSp  °q ' ' °q  &p > J-apq 


U Un  2 , v vN  2 


(A. 5) 


Substituting  Eqs.  (A. 3),  (A. 4)  and  (A. 5)  into  Eq.  (A. 2), 
and  re-grouping  the  terms,  gives 


'pq 


"pq < (#  ‘ + ‘ # ’♦2lpq  Up  ( gqv-g; ) -g“(  g“-ga ) J /apq 

- 2ypq[g^(g“-ga)tga(g^)]/ann}  ♦ 


-2  _ 2 _ 2 


pq- 


+ m (x  +y  +k  )f(g  — e ) +(e  — s ) J/a 
pq  pq  ypq  pq  L^6P  sq'  vsq  sp;  J/  pq 
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Re-grouping  the  terms  once  more,  yields 

2 2222*22 
I _ = [ra  (a  -2a  x +x  +y  +k  )/a  JL(gu)  +(gv)  ] + 
pq  pq  pq  pq  pq  pq  ■'pq  pq;/  pqJL^bp'  Vbp' 

+ [m  (a  x -x  -y  -k  )/a  ]2( gugu+gvgv ) + 

pq  pq  pq  pq  °pq  pq  pq J sp6q' 

+ (mpqi'pq/apq)2(SpSq-«p8q)  + 

+ [mpq ( '‘pq ^pq^pq  > /a pq  H < g,  > 2 U gq ) 2 ] 


The  input  torque  due  to  the  inertia  of  link  pq  can  be 
written  as 


'pq 


dKE  /d<j).  = d(0.5Iw  .)/d<j>.  = 0.5[dl  /d<t>.]<D. 


pq 


l 

3 


pq  i 


v / U,U  V, 

-pq  - Y1pq(Sphp+Sphp)“i  + 


pq 


V / U.U  U,  U V,  V V,  Vn  - 

+ Y„  (g  h +g  h +g  h +g  h )(jj. 

2pq  SP  q 5q  P SP  q °q  py  l 


v t v,  u u,  v u,  v v,  u-.  - 

+ Y3pq(8phq+*qhp'gphq‘gqhp,"i 


+ y!  (gUhU  + gVhV)a,2 
4pq  sq  q sq  q l 


(A. 6) 


where  Y^,  (i  = 1 to  4)  are  the  same  coefficients  given  in 
Eqs.  (2.8.2),  and  h^,  hp,  h^,  and  h^  are  acceleration 
influence  coefficients  with  the  following  property: 


and 


U . = 


V . 
J 


u 

h .cj  . +g  .a 


3 l 
2 


J 1 


h .a) . +e  .a  . 
J 1 S0  1 


; J = P,q 
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wnere 

Therefore 

1 


where 


= angular  acceleration  of  input  link  = 0. 
, Eq.  (A. 6)  can  be  rewritten  as 

36  36  36  36 

1pqD1pq+Y2pqD2pq+Y3pqD3pq+Y4pq°4pq 


D 


D 


D 


Ipq 

6 

J2pq 

6 

3pq 

6 

4pq 


(U  U +V  V )*«. 
P P P P i 


LVq+Vq)  + (UqVVp)j4“i 
[-(flpVVq)  + tVp-Vp)^“i 


(UqVW*“i 


APPENDIX  B 

EFFECTS  OF  A GROUNDED  LINK 
ON  SHAKING  MOMENT  EXPRESSION 


From  Fig.  B.1,  the  following  equation  can  be  written 

SM  = r xF  .k  +1  J 
12  ~1  “ JL  2 “ 12  12 

Substiting  F = m(xex+yey)  into  the  above  equation, 
results  in 

= mx xe  .kj+my[r  Xe„.kj+m(x  +y  +k  )[J  ] (B.1) 

i c.  i -A  ~ -i  -y  - 12 

Thus,  it  can  be  concluded  that  the  three  kinematic 

coefficients  inside  the  brackets  of  Eq.  (B.1)  are,  in 

general,  linearly  independent,  and  so  are  the  three 

coefficients  multiplying  the  brackets.  Clearly,  if  1' 

v 1 2 


1 

For  simplification,  the  subscripts  "12"  will  be  dropped 

out  of  m , x , y , k , and  a , since  link  12  is 
12  12  12  12  12 

the  only  one  addressed  here. 


142 


143 


V 


(a) 


V 

k i 


00 


lL 

k i 


Figure  B.1 


A Grounded  Link 

a)  Diagram  of  an  Input  Link  12 

b)  Free-Body  Diagram  of  Input  Link  12  When 
All  the  Other  Links  are  Considered 
Massless 
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constant,  the  third  term  which  corresponds  to  the  input 

torque  contribution  to  shaking  moment  vanishes,  and  link 

12  will  allow  only  two  specifications  for  the  latter 

dynamic  property.  If  "0"  is  located  at  joint  1 (r  =0), 

then  the  first  two  terms  go  to  zero  (F  does  not 

-12 

contribute  any  moment)  and  only  one  shaking  moment 
specification  can  be  made  for  link  12.  Of  course,  no 
specification  can  be  made  if  the  two  cases  above  occur 
simultaneously.  All  this  is  in  agreement  with  section 
2.8.  To  expose  even  further  the  relationships  between  Eq. 
(B.1)  and  the  formulation  given  in  section  2.8,  the 
following  rearrangement  of  Eq.  (B.1)  is  appropriate: 

S„  -r  -1”  “1-  . , _2  _2  2 

M12  = mxLr^xa  £2 .k]+my [r ^ .a  £2J+m(x  +y  +k  )U  ] 

= (m^/a)[UiV2-ViU2]+(my/a)LUiU2+ViV2]  + 

2 _ 2 _ 2 -2  2- 
+ m(x  +y  +k  )a  [a  <t>  j 

12 

Everything  stated  in  the  previous  paragraph  will  still 
hold  if  the  above  equation  is  modified  as  follows: 


■«  “"1  2 2 2 “2  ••  •• 

5M12  = { mxa  -m( x +y  +k  )a  ^U1V2“V1U2]  + 


••  ••  2 2 2 ~ 2 2 •• 

+ (my/a ) [U  U +V  V ]+m(x  +y  +k  )a  [a  * +U  V -V  U ] 
1212  J 121212 

2 2 2 “ 2 ••  •• 

= m(ax-x  -y  -k  )a  ^U1^2_V1U2 ]+(my/a) [u1u2+ViV2 ] + 
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, _ 2 _ 2 _ 2 , ” 2 r ••  ••  , 

+ m(x  +y  +k  )a  [U  V -V  U J 

2 2 2 2 

= Y3  [D4  ] +Y 3 [D4  ] +Y 3 [D4 

212  212  312  312  412  412 


] 


A similar  analysis  can  be  carried  out  for  a grounded 
link,  that  is  not  the  input  link  and  the  same  conclusions 
achieved . 


APPENDIX  C 

ELIMINATION  OF  LINEAR  DEPENDENCIES 
FOR  CASES  INVOLVING  PRISMATIC  JOINTS 


The  formulations  derived  in  sections  2.4,  2.5,  and 
2.6  for  links  with  p-<j>  motion  specification  were  left  with 
a very  attractive  configuration.  Hence,  one  can  easily 
comprehend  and  identify  the  linear  dependences  indicated 
in  Fig.  C.1  excepting,  probably,  case  4 for  shaking  moment 
and  input  torque  whose  proof  is  given  below. 


n=3  and  m=4  (shaking  moment): 


(D,  -D„  +D . )+a  = L (U  V -V  U )-(U  V -V  U ) - 

Ipq  2pq  4pq ' pq  p p p p'  v p q p q; 


- (U  V -V  U )+(U  V -V  U )]+a 
qp  q p v q q q q'J  pq 


“ [(Ep*Ep>-%*Eq)-<Eq><Ep>*RQ*E„>l-i*a 


q -q 


pq 


where  r^  and  r^  are  the  position  vectors  of  points  "p"  and 
"q"  with  respect  to  the  global  coordinate  system.  Con- 
tinuing , 


.4  4 4 

(d„_-d0__+d 


)*a  = [r  x(r  -r  )-r  x(r  -r  )].k*; 


1pq-“2pqT“4pq"“pq  ‘ i*.p*  'ip' W 'Ep'Eq'  j .«« 


= [(Ep-£q)«(Ep-Eq)]-!i«pq 
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note  2:  n=l  & m=Z  for  SF.  n=3  2.  m=4  for  SM.  n=5  &.  m=S  far  :T . 

Figure  C.1  Linear  Dependencies  Involving  Prismatic  Joints 
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= Hr  -r  )x  [*  k*  (r  -r  )-♦*  (r-r  ) J I .K*a2 


■p  -q  pq-  -p  -q  pq  -p  -q'"  - ”pq 

’pq  = °4rs  Q.E.U. 


The  proof  for  n=5  and  m=6  (input  torque)  is  very 
similar  to  the  previous  one.  One  just  needs  to  recognize 
from  Fig.  2.19  that 


D 


and 


D 


D 


i pq 

6 

2pq 

6 

4pq 


(Ep-Ep>*"i 


[ ( r . r )+(r  .r  ) ]*to . 

-p  -q  -q  -P  1 


(VEq)4“i 


APPENDIX  D 

PROOF  THAT  THE  MASS  OF  A LINK  GROUNDED  BY  A 
PIN  JOINT  NEED  NOT  BE  CONSIDERED  INDIVIDUALLY 
WHEN  CALCULATING  REACTION  FORCES 


The  shaking  force  acting  on  a linkage  due  to  the 
inertia  of  one  of  its  links  corresponds  exactly  to  the 
resultant  force  acting  on  that  link.  Therefore,  Eq. 
(2.4.1)  can  be  applied  to  link  pq  in  Fig.  D.1  in  the 
following  way: 


1111  1 

ra  r~  = Y0  Dn  +Y,  D-,  : (D„  = 0) 

pq-G  2pq  2pq  3pq  3pq  ’ v Ipq  ' 


= in  x [d>  (-s  +ic  )-<k  (c  +is  )]  + 

pq  pq L vpq ' pq  pq'  9pqv  pq  pq ' J 


+ in  y [-d>  (c  +is  ) + 6 (s  -ic  )1 

Pq^Pq L 9pqv  pq  pq'  9pq^  pq  pq'J 


Rearranging  and  considering  separately  the  two  global 
components  of  the  resultant  force  above,  yields 


= summation  of  the  U-components  of  all  external  and 
reaction  forces  acting  on  link  pq 


Xpqapq(  <*’pqSpq-(*’pq"pq^+ypqapq^~<*>pqCpq+<,)pqSpq) 


and 


EF  = x°  a (<f>  c 
Pqv  pq  pq v pq  pq 


•2  . n ” *2 

<t>  s ) +y  a ( -<f>  s — 4>  c ) (D.2) 

pq  pq  ypq  pq  pq  pq  pq  pq'  V ' 
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1130 


Figure  D.1 


A Grounded  Link  That  Fulfills  Rule  1 of 
Section  2.2 
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The  third  equation  of  dynamic  equilibrium  to  be 
considered  in  the  calculation  of  reaction  forces  can  be 
the  following: 


EM  = summation  of  all  the  moments  about  point  "p" 


■ V*pq+V+Spq>* 

Equations  (D.1)  to 
explicitly  because  Rule 
satisfied. 


pq  kpqapq<*’pq 

(D.3)  do  not  contain  mpq 
1 of  section  2.2  has  been 


(D.3) 


APPENDIX  E 

APL  FUNCTIONS  COMPOSING  THE  "GRAPHICS  SYSTEM" 


(The  majority  of  the  functions  in  this  appendix  have  been 
adapted  from  functions  created  by  Matthew  B.  Reischer.) 


152 


V GFACT  N 

[1]  n GFACT  - MAKES  REGION  N THE  ACTIVE  DRAWING  REGION  . 

[2]  n N - NUMBER  OF  THE  REGION  TO  BE  ACTIVATED  . 

[3]  GVARG+- ,GVRDTlGVRDTl ;l]\W; ] 

V 


Cl] 

[2] 

[3] 

[4] 

[5] 

[6] 

[7] 

[8] 

[9] 

[10] 
[11] 
[12] 

[13] 

[14] 

[15] 


V C GFARC  RSAI  \ NSEG\T\CO\S  \ START  ;V',M 


GFARC 

C 

RSAI [ 1 ] 
RSAll 2] 
RSAlL 3 ] 
RSAlL 4 ] 


DRAWS  AN  ARC  IN  THE  ACTIVE  REGION 
POSITION  OF  THE  CENTER  IN  LOG.  COORDS 
RADIUS  OF  THE  ARC  IN  LOGICAL  COORDS 
STARTING  ANGLE  IN  DEGREE  ( POS . OR  NEG.) 
ANGLE  DEFINED  BY  THE  ARC  (POS.  OR  NEG.) 
NO.  OF  SEGMENTS  IN  ARC  ( OPTIONAL ) 


NSEG+-2  0 


-*■(  3=pRSAl)  /JUMP 
N SEG+-RS  Al  [ 4 ] 

JUMP:T+( 2 2)p CO, (-S) , (S+-10T) ,C0*-2OT*-( RSAll 3 ]xotl 80 ) + NSEG 
V*-  (i?5AJ[l]*(  20ST  ART)  , loST  ART+-RS  Al  [ 2 ] *0  * 1 8 0 ) , ( 2,NSEG)pO 


M+-  0 


i5+ ( l^xNSEG+1 ) p'*-T+ . x7[ iM+M+ll 1 
V*-V+( pV) p( (NSEG+1 )p£[l] ) , (NSEG+ 1 ) pC[ 2 ] 

(0 ,NSEGpl)GFVABSt(  (llpA/1  l=G74flff[l0  1 1 ] ) / ' L ( ( p7 ) p . 5 ) + » ) 

,’7' 


V 


[1] 

[2] 

[3] 

[4] 

[5] 

[6] 

[7] 

[8] 

[9] 

[10] 
[11] 
[12] 

[13] 

[14] 

[15] 

[16] 


C GFCIRC  R\T\ S\CO\M\ V ; NSEG 


GFCIRC  - DRAWS  A CIRCLE  IN  THE  ACTIVE  REGION 
C - POSITION  OF  THE  CENTER  IN  LOGICAL  COORDS 

R - RADIUS  OF  THE  CIRCLE  IN  LOGICAL  COORDS ; 

OF  SEGMENTS  IN  CIRCLE 


IF  pi?=2,  Rl  2]=  NO 
ROTATION  MATRIX 
NUMBER  OF  SEGMENTS  IN  CIRCLE 


T 

NSEG 
NSEG*- 120 
-*■(  2*p  ,R)  / START 
N SE  G*-~  li  R 
J?«-  ( i 0 ) p 1 ♦ /? 

START:  T+-(  2, 2 )p  CO,S,(  -S+IOT)  ,C0*-2OT+-(  02  ) + NSEG 
V*-(0,R)  , (2,NSEG)pO 
M+-  0 


4 5+(14x^S£'ff+l)p,-*-2,+  .x7[  ; M+-M+ 1 ] 1 
V*-V+(  pV)  p(  (NSEG+1  )p<7[  1]  ) , (NSEG+1  ) pC[2] 

( 0 ,N SEGpl) GFV ABSt( ( llpA/1  l = GVARGllO  1 1 ] ) / ' L ( ( pV ) p . 5 ) + ’ ) 


V GFCLA ; T 

[1]  n CLEARS  SCREEN  AND  LEAVES  TERMINAL  IN  ALPHA  MODE  . 

[2]  □ ARBOUT  27  12 

[3]  T-+UDL  1 


V 


154 


V GFCLG ; T 

[1]  A CLEARS  SCREEN  AND  LEAVES  TERMINAL  IN  GRAPHICS  MODE  . 

[2]  UARBOUT  27  12  29 

[3]  T+UDL  .5 

V 


[1] 

[2] 

[3] 

[4] 

[5] 
C6] 

[7] 

[8] 

[9] 

[10] 
[11] 
[12] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 

[19] 

[20] 
[21] 
[22] 
[23] 


V PI  GFDALINE  P2 ; A ; B ; P1P2PH ; R12 \MAG ; R12UNIT ; RATIO ; INT ;V1 ;V 
2 ;A iPViIBiliKEEP 

a GFDALINE  - DRAWS  AN  INTERMITENT  LINE  SEGMENT  FROM  PI 
TO  P2U0G. 

fl  COORDS.)  IN  THE  CURRENTLY  ACTIVE  REGION. 

A DASHTYPE  - GLOBAL  VARIABLE  THAT  DEFINES  THE  TYPE  OF  T 
HE  LINE  SEG : 

A DASHTYPE*!  DOTTED  LINE 

A =2  -*■  SHORT-DASHED  LINE 

A =3  -*■  MEDIUM-DASHED  LINE 

A =4  -*■  LONG-DASHED  LINE 

A =5  -*■  DOT-DASHED  LINE 

Z'P1P2PH*-'  , ( (5pv/l  1 0 0*KEEP+GVARGI  10  11  12  13  ])/'GFLTP 

' ) , * PI, [1. 1 ] P2  f 

R12UNIT+-R12iMAG+-(  + / (R12-* / P1P2PH  ) * 2 )*  . 5 

I«-0 

GVARGllO  11  12  13>1  100 
AGAIN:  A*- 1 6 8 10  12  1 [ DASHTYPE+Il 
B+-  6 4 5 5 12  SlDASHTYPE+I] 

INT+lRATIO*-MAG*A  + B 

Vl+l  . 5++\PlP2Ptf[  ;1]  INT,  2 )pRl  2 UNIT*A+B 

V2<-1  . 5 + + \ (P1P2PBC  ; 1 1 + A*  R12UN  IT)  ,A 

PV+( 0 "1+ ( 2 , 2*INT+1 )p71 , [2 . 1]72 ) ,P1P2PB[ ;2 ] 

IB+-(  2 X-INT+  1 ) pO  1 

IB[  1 + 2* INT]-*- (RATIO- INT)  £ . 3 

IB  (7P7i4B5i(  ( 12p  a/1  l = GIM/?ff[10  1 1 ])/'  L ((  pPV ) p . 5 ) + »),’  PF  ’ 
-^/4C^I^x1I^(I=0  ) f^D ASHTYPE-  5 
G7.4BS[10  11  12  13  l+KEEP 


V LL  GFDARECT  UR 

[1]  A GFDARECT  - DRAWS  A RECTANGLE  WITH  INTERMITENT  LINE  SEG 
MENTS. 

[2]  a DASHTYPE  - GLOBAL  VARIABLE  DEFINING  THE  TYPE  OF  THE  LI 
NE  SEG. 

[3]  A ( SEE  COMMENTS  IN  ' GFD A L IN E » AND  ' GFRECT' ) 

[4]  LL  GFDALINE  BP[l],II[2] 

[5]  (Z7tf[l]  ,£I[2]  ) GFD  ALINE  UR 

[6]  UR  GFDALINE  II[l],Bff[2] 

[7]  (LI[1] ,URl2l)GFDALINE  LL 


V 
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V GFDREG  R 

[1]  n GFDREG  - DELETES  REGION (S)  FROM  THE  REGION  DEFINITION 
TABLE  . 

[2]  A R - NUMERIC  SC  ALAR  (.OR  NUM  VECTOR ) SPECIFYING  THE 

REGION (S)  TO 

[3]  r BE  DELETED  FROM  GVRDT. 

[4]  GVRDT<-(~GVRDTl  illeR) /ZllGVRDT 

V 


[1] 

[2] 

[3] 


GFDRW  P\ B 

A GFDRW  - DRAWS  A LINE  SEGMENT  FROM  CURRENT  POSITION 
fl  OF  GRAPHIC  CURSOR  TO  P.  THE  TERMINAL  MUST 

a ALREADY  BE  IN  GRAPHICS  MODE , AND  WILL  BE  LEFT 


[4] 

[5] 

[6] 

[7] 


V 


A 

A P 
A 


IN  GRAPHICS  MODE. 

A POINT  EXPRESSED  IN  LOGICAL  COORDS.  IN  THE 
ACTIVE  REGION  TO  WHICH  THE  LINE  SEGMENT  IS  TO 


BE  DRAWN  . 

1 GFVABS  £ ( (llpA/1  1 0 0 = GVARGll0  11  12  1 3 ] ) / ’ L ( ( pP ) p . 5 ) 
+' ) , 'P» 


[1] 

[2] 

[3] 

; 4 : 

5] 

L6] 


Z+GFGCURitiTXiA 

A GFGCUR  - BRINGS  UP  THE  CROSSHAIR  CURSORS 
A Z - Z C 1 ] X-COORD 

a Z [ 2 3 Y-COORD  ( LOGICAL ) 

A Z[3]  UAV\(WHATEVER  KEY  YOU  HIT) 

V\TX+-  ’ UARB OUT  7 0-*-5  ' 

UARBOUT  29  27  26 


[7]  Z-«-  <t>  3 2 64  4>  3 2 1§2  2p  ( ~ 1 +(2A  V i 1 + A+B ) - 3 2 96  32  64 

[8]  Z-<-  ( GFPTL  Z)  ,UAV\1*A 
V 


V PI  GFLINE  P2\V 

[1]  A GFLINE  - DRAWS  A LINE  SEGMENT  FROM  PI  TO  P 2 . THE  TER 
MINAL  MUST 

[2]  r ALREADY  BE  IN  GRAPHICS  MODE , AND  WILL  BE  LEF 

T IN  GRAPHICS  MODE  . 

[3]  A P1,P2  - POINTS  EXPRESSED  IN  LOGICAL  COORDINATES 

[4]  n IN  THE  CURRENTLY  ACTIVE  REGION  . 

[5]  0 1 GFVABSi.(  ( llpA/l  l 0 0 = GVARGll0  11  12  1 3 ] ) / ' L ( ( pV ) p . 5 
) + ’ ) , 1 V*-P1  ,[1.5  ]P2  ’ 

V 

V PP+GFLTP  LP 

[1]  A GFLTP  - CONVERTS  THE  LOGICAL  COORDINATES  OF  THE  POINT 

(5)  SPECIFY I ED 

[2]  A BY  LP  TO  PHYSICAL  SCREEN  COORDINATES  . 

[3]  A LP  -A  2 *N  ARRAY  OF  LOGICAL  COORDS  OF  N POINTS  IN 
THE  CURRENTLY  ACTIVE  REGION  . 

[4]  a IF  N=1  EITHER  pLP=2  1 OR  pLP=2  WILL  BE  OK  . 

[5]  a PP  - CONVERTED  PHYSIC AL  SCREEN  COORDS.  (pPP=pLP)  . 

[ 6 ] i l = p p LP 


[7]  PP^( L . 5+GVARGl 12]+SPCl ; ]*GVARGllOl ) , [ . 5]L . 5 +GVARGI 1 3 ] + LP 
C 2 ; ] x GV  ARGL 11] 

[8]  ->0 

[9]  L:PP*-l  .5  + (~2  + GVARG)  + LP*GVARGllO  11] 

V 

V ffFA/07  P 

Cl]  A GFMOV  - MOVES  THE  CURSOR  INVISIBLY  TO  THE  POINT  P . 
[2]  A THE  TERMINAL  WILL  BE  LEFT  IN  GRAPHICS  MODE  . 

C3]  A P -A  POINT  (PC  1], P[ 2])  EXPRESSED  IN  LOGICAL 

C4]  A COORDINATES  IN  THE  CURRENTLY  ACTIVE  REGION  . 

C 5 ] 0 GFV  AB  S ±( ( 1 lp  a /l  1 o 0 = GVARGLl0  11  12  13 ] ) /' L ( ( pP) p . 5 ) + 

' ) , 'P' 

V 

7 GFPBUF  RS 

Cl]  A GFPBUF  - PURGES  THE  WHOLE  BUFFER  IF  0 e RS,  OTHERWISE 

C 2 ] A PP4SSS  StfPJ  TPff  COLUMNS  OF  GVBUF  CORRESP. 

C 3 ] A TO  THE  REGIONS  INDICATED  BY  RS  . 

C 4 ] A RS  - NUMERIC  SCALAR  OR  NUMERIC  VECTOR  . 

C 5 ] ->-LBxt~0ePS 
C 6 ] GVBUF+ 4 lpO 
C 7 ] +0 

C 8 ] LB  :GVBUF+(~GVBUFHnleRS)  / GVBUF 

V 

V L+GFPTL  P 

Cl]  A GFPTL  - CONVERTS  PHYSICAL  COORDINATES  IN  P TO  ITS 
C 2 ] a CORRESP.  LOG.  COORDS.  IN  THE  ACTIVE  REGION. 

C 3 ] a P -A  2 xtf  ARRAY  OF  PHYSICAL  COORDS.  OF  N POINTS  0 

F THE  SCREEN. 

f>]  n IF  N= 1,  EITHER  pP=2  1 OR  pP=2  WILL  BE  OK. 

C 5 ] A L - CONVERTED  LOGICAL  COORDINATES  (pL=pP). 

C 6 ] +((ppP)=2) / LINE 

C 7 ] L+-GVARG  C6  7 ]+  ( P- GVARG  C 2 3 ] ) * S7/SPGC  1 0 11] 

C 8 ] -*-0 

C 9 ] LINE: L+GVARG C 6 ]+ (PC  1 i 1-GVARGl 2 ] )iGVARGl 10 ] 

CIO]  L+L.l.llGVARGZ7]+(Pl2i]-GVARGZ3])*GVARGlll] 

V 

V Z+-GFRBUF  RS;  Q\  T\PV ; IB 

Cl]  A GFTBUF  - TRANSMITS  THE  WHOLE  BUFFER  (.GV  BUF)  IF  0 € RS 
C 2 ] A OTHERWISE  TRANSMITS  ONLY  REGIONS  INDICATED 

C 3 ] a BY  RS  AND  FOUND  IN  GVBUF  . 

C 4 ] A PS  - NUMERIC  SCALAR  OR  NUMERIC  VECTOR  . 

C 5 ] -»-£l  x i ~0ePS 

C 6 ] Q+(  l+~ 1 + p GV BUF) pi 

C 7 ] -*-L  1 + 1 

C 8 ] Ll:Q+liGVBUFl>±;]eRS 
C 9 ] IB+Q/liGVBUFLl ;] 

CIO]  PV+Q/0  WGVBUF12  3;] 

C 1 1 ] T+,§( 2 9 x~ip ) , C 1 ]7C ; 1 +P7C  2 ; ] ] , C 1 ]XC ; 1 +P7C 1 ; ] ] 

C 1 2 ] Z*-(T*0)/T 

C 1 3 ] GVCUR+-,  2 ~lJtPV 

V 
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[1] 

[2] 

[3] 

[4] 

[5] 

[6] 

C 7 ] 

[8] 

C9] 

[10] 


V LL  GFRECT  UR; I 

A GFRECT  - DRAWS  A RECTANGLE  . 

A LL  - ( X ,Y ) COORDINATE  OF  THE  LOWER  LEFT  CORNER  . 

A UR  - (X,Y)  COORDINATE  OF  THE  UPPER  RIGHT  CORNER. 

A BOTH  COORD'S  EXPRESSED  IN  LOGICAL  COORD'S 

A IN  THE  CURRENTLY  ACTIVE  REGION  . 

■+LB*  i ( l = ppLL)v3=+/p£L 

(JpO  111  l)GFVABS$((I+-5x~l  + pLL)  ,2)p($LL)  ,LLll;l  ,URl2;l 
,(*tUR)  ,URil;l  ,LLl2;l  ,§LL 
-*■0 

LB : UR-*- , UR 

0 1111  GFVABS  LL,(LLlll,URl2l)  ,UR , (URZll ,LLl2l)  .[1.5] 
LL+,LL 


V 


V N GFREG  C 

[1]  A GFREG  - ASSIGNS  A NUMBER  AND  PHYSICAL  CO'ORDS  TO  A R 

EGION  . 

[2]  A N - REGION  NUMBER  . 

[3]  « C - PHYSICAL  CO'ORDS  (LLX ,LLY ,URX , URY ) . 

[4]  a GVRDT  - REGION  DEFINITION  TABLE ( GLOBAL ) : 

C5]  A (NUMBER, LLXP, LLYP ,URXP,URYP,LLXL,LLYL,URXL,URYL ,S 
CLX ,SCLY ,0  RGXP ,0  RGYP) 

[6]  +l+nLCll]+2*NeGVRDTl ;1] 

[ 7 ] GVRDT+-GVRDT , [ 1 ] ( p GVRDT)  [ 2 1*N  ,C 

[8]  +0 

[9]  GVRDTlGVRDTl  ; 1 ] t ; 1 + 1 4 

V 


7 N GFSCL  T 

[1]  A GFSCL  - MODIFIES  THE  REG.  DEFINITION  TBLE,  GFRDT,  BY 

[2]  A ASSOCIATING  LOG.  COORDS.  WITH  THE  REGION  N 

[3]  A INITIALIZED  BY  GFREG. 

[^3  n T ~ LOGICAL  REGION  COORDS.  (LLX ,LLY ,URX ,URY) . 

[5]  N+-GVRDTI  ;1  ] i 

[6]  GVRDTlN  \5  + \ 6]^r,  ( ( GVRDTiN  ; 4 ] - GVRDTlN  ; 2 ] ) * T\_  3 ] -T[  1 ] ) , (GVR 
DT[N;5]-GVRDT[N;3])tT[4]-T[2] 

[7]  GVRDTlN  ; 12  131+-GV  RDTlN  -,2  3]-T[l  2 ] x GVRDTlN  ; 1 0 11] 

V 


[1] 

[2] 

[3] 

[4] 

[5] 

[6] 

[7] 


V GFSTUP 

A INITIALIZES  REGION  DEFINITION  TABLE(GVRDT)  BY 
a ASSIGNING  TO  IT  THE  BASE  REGION (REGION  1):  FULL 
A SCREEN ; BOTH  LOGICAL  AND  PHI  SIC AL  ORIGINS  AT 
A THE  LL  CORNER ; BOTH  X AND  Y SCALE  FACTORS  = 1 . 

A INITIALIZES  BUFFER(GVBUF)  AND  SETS  THE  TRANSMISSION 
A CODE(TRC)  TO  AMEND  AND  TRANSMIT  . 


GVRDT* 
GVBUF * 


•(1  1 3 ) p 1 , 0 
(4,l)p0 


0,1023,779,0,0,1023,779,1,1,0,0 


TRC+' AT' 


V 
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V GFTBUF  RSiQiTiPViIB 

[1]  r GFTBUF  - TRANSMITS  THE  WHOLE  B UFFER (GV BUF)  IF  0 e RS 

[2]  R OTHERWISE  TRANS.  ONLY  THE  REGIONS  INDICATED 

L3]  R BY  RS  AND  FOUND  IN  GVBUF  . 


[4] 

R RS  - NUMERIC  SCALAR 

OR  NUMERIC  VECTOR  . 

[5] 

-+Ll*i~0eRS 

[6] 

G«-(  1 + ~1  + p GVBUF)  p 1 

[7] 

-*-L  1 + 1 

[8] 

LI iQ+l+GVBUFlUileRS 

[9] 

IB+Q/1\GVBUF\_1-,  ] 

[10] 

PV+-Q/ 0 l\GVBUF\_2  3;] 

[11] 

2’-,6?(2  9x~jfl)  , [ l ]y[  ;1+P7[2;] 

].[1]J[ ;1+P7[ 1 ; ] ] 

[12] 

UARB0UT(T* 0 )/T 

[13] 

V 

GVCUR+- , 2 "1  \PV 

7 

P GFTXT  TEXT 

[1] 

R MOVES  CURSOR  TO  P[ 1 2]  ( 
ACTIVE  REG.) . 

LOG .COORDS . OF  A POINT  IN  THE 

[2] 

fl  IF  pP=  3 , BACKSPACES  THE 
3].  THEN 

NUM.  OF  TIMES  SPECIFIED  BY  P[ 

[3] 

fl  DISPLAYS  THE  CHAR.  VECTOR  TEXT. 

[4] 

GFMOV  2\P 

[5] 

OARBOUT  31 

[6] 

-*■(  2+0LC\.  l])xt2  = pP 

[7] 

OARBOUT  P[ 3 ] p 8 

[8] 

TEXT 

7 


V IB  GFVABS  XY;T;PV 

Cl]  R GFVABS  - TAKES  AN  X ,Y  VECTOR  OF  COORDS  AND  DRAWS  THEM 
AS  ABSOLUTE 

[2]  r COORDS  IN  THE  ACTIVE  REGION. 

[3]  r XY  - ARRAY  OF  ABSOLUTE  COORDS,  pXY=(2,N)  . 

[4]  n IB  - CONNECTIVITY  VECTOR ; IB[I]  = 0:  MOVE  CURSOR 

INVISIBLY  TO  ZY[;J] 

[5]  R =1:  DRAW  FROM  CU 

RRENT  CURSOR  POSITION  TO  XYlill 

[6]  +2  + ( l + DLOx  x 2 = pps.'PV+'  , ( ( 5p  v/l  l 0 0*GVARGll0  11  12  13]) 
/'GFLTP'),'  XY' 

[7]  PV+-  2 IpPV 

[ 8 ] -+ER*  i v/  ( ( T /PV)>  1023,779  ) , ( L/P7)<0 

[9]  GVCUR+,2  ~1 \PV 

[10]  -*-(  1 eT_RC_)  /2  + 1 + ULC 

[11]  GVBUF+GVBUF  ,IB  , [ 1 ~\PV  , [ 1 ] ( p ,IB)  pGVARG\_  1 ] 

[12]  +(~'T' eTRC) / 0 

[13]  ^,fi?(29x~JS)  , [ 1 ]J[  ;1+P7[  2 ; ] ] , [ 1]  J[  ;l+PVl  1 ; ] ] 

[14]  OARBOUT(T*0 )/T 

[15]  +0 

[16]  ERiUARBOUT  13 

[17]  ’ GRAPHICS  ERROR:  ATTEMPT  TO  DRAW  OFF  OF  THE  SCREEN' 

[18]  ) SI 

[19]  - 


V 
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[1] 

[2] 

[3] 

[4] 

[5] 
C6] 

[7] 

[8] 
[9] 
CIO] 
[11] 
[ 12] 

[13] 

[14] 

[15] 

[16] 

[17] 

[18] 
[19] 


V IB  GFVINC  XY;T',QlO;PV 

A GFVINC  - TAKES  AN  X ,Y  VECTOR  OF  COORDS  AND  DRAWS  THEM 
RELATIVE 


A XY 
A IB 

BLY  TO  XY 


TO  THE  CURRENT  PO S . OF  THE  GRAPHICS  CURSOR. 

- ARRAY  OF  RELATIVE  COORDS , pXY=(2,N) 

- CONNECTIVITY  OF  LINES,  0:  MOVE  CURSOR  INVISI 


A 

CURSOR  POSITION  TO  XY 


DRAW  FROM  CURRENT 


+2  + (l  + dLC)xi2  = ppP7«-0  14  + \ ( 2 IpGVCUR)  , i(  ( 5p  v/l  l 0 0 XGVAR 
<?[10  11  12  13])/'  GFLTP1  ) , ' XY ' 

PV-*-  2 lp  PV 

+ER*\ v/( ( £ / PV )>1023,779) , ( L/P7) <0 
GVCUR+,2  ~1 \PV 
+(-' A' eTRC) /2+1+ULC 

GVBUF+GVBUF , IB ,ll]PV ,ll]( p, IB) pGVARGl  1] 

~ 1 21 1 eTRC ) /0 

T+-  ,$(29  x~  jp ) ,[1]J[  ;1+P7[2;]]  ,[1]J[  ;1+P7[l;]] 

UARB0UT(T* 0 )/r 
-*•0 

13 

’ ffP/lPPICS  ERROR:  ATTEMPT  TO  DRAW  OFF  OF  THE  SCREEN ' 

)SI 


V 


V ffP7P  NS',  I ',  II 

[1]  A DRAWS  VIEWPORT  BOUNDARY  ASSOCIATED  TO  THE  REGION (S)  I 
NDICATED  BY  NSiSCALAR 

[2]  a OR  VECTOR  BUT  ALWAYS  NUMERIC). 

[3]  NS+( l,I-*-p,NS)pNS 

[4]  GFACT  Ii>( iO)ptfS[l;I] 

[5]  (GVRDTlGVRDTl ;1 ]ill;6  7 ] ) GFRECT  GVRDTl GVRDTl ; 1 ] i II; 8 9] 

[6]  -►-2  + l + []Z;C’xO*r«-I-l 

V 


APPENDIX  F 

APL  FUNCTIONS  COMPOSING 
THE  "APPLICATION  PROGRAM" 


V PHI+X  AEG  Y ; Z 

[1]  A CALCULATES  ANGLE(S)  BETWEEN  X ,Y  VECTOR(S)  AND  X-AXIS. 

[2]  PHI*-(  (O.  5)*Z**Y)  + (~Z)*~3OYiX+Z<-X=0 

[3]  PH I+-PH  1+02*0  > PH I+-PHI+(X<0  )x(oi)x(  x7)+y=o 

V 


7 Z-+BASICMP;  L\MUV\CM  \CXY 

[1]  A TRANSFORMS  LUMPED  MASS  PARAMETERS 

[2]  A INTO  BASIC  MASS  PARAMETERS 

[3]  NM12X-+M12MIN  tM12MAX'>  Z-0 

[5]  -(LS2 ,LB1 ,QLC[ l]+l )[ 2+*X012] 

[ 6 ] AMI  2X+-  1+2*  (Ml  2MIN*X1  2 MAX  )>L)+M  1 2 MIN  , ( L+-XO 1 2*A12)iX12MA 

X 

[7]  AM12X+-AM12X  , l8  + (8*(M12MAX*X12MIN)  > L)  <$>'  M12MAX  L*X12MIN' 

[8]  -LSI 

[9]  LB2:AM12X+(~1  + 2*(M12MIN*X12MIN)<LUM12MIN,  ( L-+X01  2 * A 1 2 ) 4 XI 

2 MIN 

[10]  AM12X+-AM12X  ,t8i  (8*(M12MAX*X12MAX)<L)<p'  M12MAX  TAX12MAX' 

[11]  LSI  : AM12Y+-M1 2 MIN  ,M12 MAX 

[12]  — (LS4,LS3,QLL’[l]  + l)[2  + xy,912] 

[13]  AMI  27— ( 1 + 2 x ( M 1 2 MIN * 7 1 2 M AX ) > L ) 4 Ml  2 MIN , ( L- YO 1 2 * A 1 2 ) 4- 7 1 2 MA 
X 

[14]  AM12Y+-AM12Y  (8*(M12MAX*Y12  MIN ) > L ) <j> ' Ml  2 MAX  L*Y12MIN  ' 

[15]  — LS3 

[16]  LS4  :AM12Y+(~1+2*(M12MIN*Y12MIN)<L)\M12MIN  , (L+-Y01  2*A12  )*71 

2 MIN 

[17]  AMI  2Y+-AM1  27  , ±84  ' 8 x ( VI  2M4Xx  71 2M.4X ) < L ) 4> ' Ml  2MAX  LH12MAX' 

[18]  LS3  :-+(0  = ppAM12+COMMON3  A Ml  2X  , AM  1 27  , Ml  2MIN  ,M1  2MAX ) 4 0 

[19]  -( ( X2  3 MIN<X 2 3 ) *X2  3MAX>X2 3- ( X 02 3 x A 2 3 ) tM02 3 ) /0 

[20]  -(  (Y23MIN<Y_2  3 )*7  2 3M.4X>72  3-(  7(72  3 x A 2 3 ) *M02  3 ) / 0 

[21]  -(  0 >MUV-+(  (K023*A23*2)iM023  )-+/(X23,723)*2)/0 

[22]  -(  (K23MINHK23  )*X23MAX>K23+-MUV*  . 5)/0 

[23]  AMM3X*-MM3MIN  ,MM3MAX 

[24]  -(LS6,LS5,QLC’[l]  + l)[2  + xXC>43] 

[ 25  ] AM43X-C  1 + 2 x ( M4  3 MIN*  74  3 MAX)  > L)  4 M4  3 MJ.V  , ( 0X04  3x44  3 ) 4X4  3M4 
X 

[25]  AM43X-AM43X, ±84  ( 8*(MM2MAX*XM2MIN)  >L)$' MM3MAX  Li XM 3 MIN' 

[27]  — LS5 

[28]  LSS  : AM43X-(  1 +2* (MM3MIN* XM3MIN ) < L) i MM3MIN , ( L+-XOM3*  AM3  ) i XM 

3 MIN 

[29]  AM4  3 X—AM4  3 X , ± 8 4 ( 8*  ( MM3MAX*  XM3  MAX  ) < L)  <$> ' MM3  H AX  HXM3MAX' 

[30]  LB 5 :AMM3Y-MM3MIN ,MM3MAX 

[31]  — (LS8 ,LS7 ,DLC[l]  + l )[2  + xy.04  3 ] 

[ 32  ] A/V4  37-(~l+2x  (MM3MIN*YM3MAX)  > L)  i MM  3 MIN  , (L-+YOM3*AM3  )4743M4 
X 

[ 33  ] AM437-AM437  , ±84  ( 8 x ( MM  3 MAX*  YM3MIN ) > L ) 4>  • MM  3 MAX  HY  M3  MIN  ' 

[34]  — LS7 

[35]  L88  : AM4  37—  ( 1 +2* ( MM3MIN* YM3MIN ) < L) i MM3  MIN , (07043x443  )*74 

3 MIN 

[36]  AM437-AM437  ,±84  ( 8 * ( MM3  M AX*  YM3M  AX ) < L)  <$> ' MM3MAX  LHM3MAX' 

[ 37  ] LB  7 :-(  0 =p  p AMM3-+C  OMMON2  AMM3X  , AMM3Y  , MM  3 MIN  ,MM  3 MAX ) 4 0 
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[38]  K^3MINS*-K^3MIN*20  KH3MAXS+K^3MAX*2 

[40]  CK+KO^  3 x j44  3 * 2 

[41]  ->-(  O^CXY-'-C  04  3*2  )+(Y043*2  ) ) * AH  3*  2 ) + LB9 

[42]  -*■(  ( {M'\3MINxKK3MINS)>CK)'tCK>M'\3MAX*K^3MAXS)1t3 

[43  ] A/tf4  3K+-  (-l  + 2x(A/4  3MIN*K,+  3MAXS)  >CK)  *M13MIN  ,CK*K*3MAXS 

[44]  AM4  3K*-AAf4  3K,ilO  + ( 10* (MH3MAX* KH3MIN S)  >CK)  <t>'  M^3MAX  CK*K 

43 MINS' 

[45]  +LB 10 

[46]  LB9  :-*■(  ( KH3MAXS<K2P ) , ( ( KH3MINSSK  2P)  *K2PZKH3MAXS)  ,K^3MINS>K 

2P*-(CK*2)m*CXY)  / (LB9+1)  , LSI 1,0 

[47]  A.V4  3JT«-(  ( ( CK-MUV ) ,CK+MUV->-(  (C  K*  2) ->A  *CXY*  L)  * . 5 ) * 2 x £ , L+KH3MI 
NS ,KK3MAXS)ll  243] 

[48]  L+-COMMON 2 AM43£[3  4 ] . A-V4 3 

[49]  -*-(  ( Q = ppL)*0=ppMUV<-COMM0N2  AM43£[l  2],A.V43)tO 

[50]  AM43^(  ( (0=ppL)40  0 * p p L ) , ( 0 = p p Af  JJV  ) + 0 0*  p pMUV  ) / L ,MUV 

[51]  -*(Z*-l)/0 

[52]  LSI  1 : AM43£*-(  ( CK-MUV ) ,CK+MUV+(  (C K*2) *CXY* KH3MIN S) * . 5)*2* 

K^3MINS 

[53]  LB10  ;-*■(  0 = pp  LM^3+-C  0MM0N2  AM4  3 K , AM4  3 ) /O 

[54]  Z«-l 
7 


7 BBOIff  00; LO 

[1]  n AMI.V  PROGRAM 

[2]  +00+00=1 

[3]  GFSTUP 

[4]  OFOL/1 

[5]  ' PLEASE  SELECT  MECHANISM :' 

[6]  ’ 1)  FOUR-BAR ' 

[7]  ’2)  SLIDER-CRANK1 

[8]  * 3)  (vMTr  II* 

[9]  ’ 4)  STEPHENSON  III1 

[10]  -*■(  ( 0 = p pM_TY)  a a /M_TY  €.  1 1234'  )+("l  + » 1234’  i 1*M_TY+-B)  ©LB1  , LB2 

, LB3  , LB4 

[11]  CM  BBO 07  7 

[12]  -□LC[l]-2 

[13]  LBhFBSPEC 

[ 14]  ’ PLEASE  WAIT . . . ' 

[15]  DASHES 

[16]  LBS:ERR+0 

[17]  FBDVA  1 ~1  [ S.TYPE1 

[18]  -+LB5  * i ERR-  0 

[19]  GFCLA 

[20]  'AN  ERROR  OCCURRED  DURING  THE  DISPLACEMENT  ANALYSIS 1 

[21]  'OF  THE  H-BAR.  FROM  THE  LIST  BELOW,  SELECT  ALL  ITEMS' 

[22]  'YOU  WANT  TO  CHANGE  ( USE  SPACES  TO  SEPARATE 11  EM ) OR' 

[23]  ’ ENTER  4 FOR  QUITTING.' 

[24]  ’ ’ 

[25]  ’ 1)  U 1 VI  U 4 K4  = ’,(t01),’  ’,(T71),’  ’,(tO 4),’  » ,t74 

[26]  ’ 2)  LENGTHS  OF  12,  23,  AND  43  = ’,(tM12),’  ’,(t423),'  ’ 

, T /44  3 

[27]  ’ 3)  TYPE  OF  DYAD  234  IS  ’,(2  ~2 [ iTYPEl ) * ’ +1 -1 ’ 
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[28] 

[29] 

[30] 

[31] 

[32] 

[33] 

[34] , 

[35] 

[36] 

[37] 

[38] 

[39] 

[40] 

[41] 

[42] 
[43  ] 

[44] 

[45] 
.46] 
[47] 


' 4)  QUIT 1 

QARBOUT  SI 
FBTXT 

UARBOUT  lOplO 

-+(DLO[  l]  + 4 )X!  A/(M=4  )v(¥P7«-[]0  e i 3 
UARBOUT  7 
’ TRY  AGAIN 1 
-□LO[  l]-3 

->(  (□SC[1]  + 1 ) ,0)[ l+A/¥07e4] 

FBOFO  ¥07 
-+IB6 

LB5 :AFLAG+12  2p'AT' 

OFOL/1 
’ ENTER:1 

1 1)  FOi?  GRAPHIC  DISPLAY  OF  DATA  FROM 1 

’ OIO.  404LYSIS;’ 

’ 2)  IF  A/O  OISPI4Y  IS  DESIRED 

* 3)  FOB  QUITTING 1 

-+(  ( 0=p p¥07) a a /¥07e  ’ 123*  )+ ( 1 + ’ 123*1  l + ¥07-’-Q  ) ©IB7  , IB8 
□4BBOOI  7 


0 


[48]  ’IBY  40410’ 

[49]  +DL0[l]-3 

[50]  LB1 :FBKOUT 

[51]  GF  ACT  1 

[52]  OF, ¥07  0 7 70 

[53]  D4BBOOI  13 

[54]  -+LB5  +2 

[55]  IB8  :F12+-F2  3«-F4  3«-2  lpO 

[56]  F12+-F2  3-«-F4  3-<-(  2 ,0+1  ) pO 

[57]  112^-123^143^(0+1  )pO 

[58]  L0*-I04PS 

[59]  FB  ST  ATIC 

[60]  -(10=0 ) + DLO[ 1 ] + 3 

[61]  CASE+'L' 

[62]  STADYNOUT 

[ 63  ] ¥1  2— ¥2  3 -+¥4  3+-Z1  2-Y2  3-14  3-Y1 2-Y2  3-Y43-B1 2-02  3-043-0 

[64]  IB9  :¥PIO 

[65]  FBDYNAID 

[66]  FBDYN 

[67]  FBTOTAL 

[68]  CASE+'T' 

[69]  STADYNOUT 

[70]  * WOULD  YOU  LIKE  TO  RETROCEDE  AND1 

[71]  1 MODIFY  THE  CURRENT  MASS  PARAMETERS? 1 

[72]  i’-’  , 5+ ("5  + 5x’ Y0*  i 1 + 0)  4>f  IS9  □I0+3DI0+1' 

[73]  UARBOUT  7 

[74]  -C|L0[l]-2 

[75]  MP LIMITS 

[76]  — OLP[  1 ] + 2 

[77]  LBlliMPLIMMOD 

[78]  FBMPOLIM 

[79]  'INDICATE  ALTERNATIVE  YOU  WANT  TO  ADDRESS:1 

[80]  ’ 1)  SHAKING  MOMENT  SYNTHESIS 1 

[81]  ' 2)  SHAKING  FORCE  AND  INPUT  TORQUE  SYNTHESIS 1 
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[82]  -*•!  5t  ( 5 + 5 x ♦ 21 ' i USYNCASE+U)<$>'  LB10  D££+3ri££+l  ' 

[83]  UARBOUT  7 

[84]  -[]££[l]-2 

[85]  SMRP 

[86]  £512  :+LBll*\Q=ITSMSYN 

[87]  SMT+CURVE 

[88]  £51 3 : WFWIN 

[89]  EXPLN_ 

[90]  -►(£B11,£B12,£B13,0)[  tFBOPTl 

[91]  £51 0 : SMRP 

[92]  £514  :-*LBll*\  0=ITSMSYN 

[93]  ITT+CURVE 

[94]  -*-£51 1 x i 0 -SFSYN 

[95]  SFT+CURVE 

[96]  £51 5 : WFWIN 

[97]  EXPIN_ 

[98]  -*(£511, £514, £515, 0)[  iFBOPT  ] 

[99]  £52: 

[ 100 ] £53 : 

[ 1 0 1 ]  £54 : 

V 


7 CALC  L;UTX 

[1]  n FUNCTION  CALLED  BY  'CALCULATOR' . PERFORMS 

[2]  n OPERATION  INDICATED  BY  L. 

[3]  UT X*-'  ERR1+1O+0  ' 

[4]  -»-£16xlL=16 

[5]  TT+TO  ZZ+ZO  YY+YO  XX+XO  VVSTK+VSTKO  +sl'L',t>L 
[11]  £1  : X+X*  . 50  -*-££’1 

[13]  £2  : X+1OX0  +LC1 

[15]  £3  :X+2oXO  -*-££1 
[17]  £4  :X+X*20  -*-££1 
[19]  £5  :**-”lo*0  +££1 
[ 21  ] £6  :X+  2OX0  -*-££1 
[23]  £7  :X+*XO  +LC1 
[ 25  ] £8  : X+10*XO  -►££! 

[27]  L9:X+Y*XO  +LC2 
[29]  £1 0 : X+&XO  -*-££1 
[31]  £11:  X+l  0&XO  ->££1 
[33]  £12: X+*XO  +LC 1 
[ 35  ] £1  3 : X+0OX0  +LC1 
[37]  £14:  X-*-  3OX0  +LC1 

[39]  ££1  :VSTKL*  ;>3+(  ' R'  ,tR+R+1  ) , ' '0+0 

[41]  £1  5 : T-^ZO  Z+YO  Y+X 

[44]  VSTK+VSTKL2  3 4 ; ] , [ 1 ] VSTKl  4 ; ]0  -*-0 

[46]  £16: +ER* \ (11=0 ) a J=  1 

[47]  MUV+X 0 X+XXO  XX+MUV 
[50]  MUV+YO  Y+YYO  YY+MUV 
[53]  MUV+ZO  Z+ZZO  ZZ+MUV 
[56]  MUV+TO  T+TTO  TT+MUV 

[59]  MUV+VSTKO  VSTK+VVSTKO  VVSTK+MUV 
[62]  +(l+0W[l])x(^mTi[4;l])«'ff5' 
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[63]  ±MUV  ,MUV  - 1*0  +0 

[65]  L17:T*-ZO  Z+YO  Y+X 0 X+IANG 

[69]  VSTK+VSTKl  2 3 4 ;],[1],I4  '0  ->0 

[71]  L18-.GFM0V  Q,H+H-10 

[72]  UARBOUT  13  27  59 

[73]  T+-ZO  Z<-YO  Y+-XO 

[77]  VSTK+VSTKl 2 3 4 ; ] , [ 1 ] 3 + ( 1 N ' , fN+N+ 1 ) , » * 

[78]  UARBOUT  27  560  -*-0 
[80]  £19 : X+Y*XO  +LC2 

[ 82  ] L20:X+-Y*XO  -+LC2 
[ 84 ] L2l:X+Y-X0  +LC2 

[86]  L22-.X+Y+X 

[87]  LC2-.Y+-Z0  Z+-T 

[89]  VSTK+VSTKllil  , [ 1 ] KSM  1 2 ;]  ,[l]3  + ( 'R'  ,tJ?^J?+1  ) , » '0+0 
[91]  ER-.ERR2+- 10  -*-0 
V 


V W+CALCULATOR  IN PANG ; II ; I ; N i KIM ; LV ; ACC_I A ; UN IT ; RV ; MUV ; I AN 
GiX;RiY;Z;Ti VSTK ; H ; CHPOS ; COL ; ROW ; ERR  1 ; ERR2 

[1]  n CALCULATOR  - CALCULATES  VALUE  OF  A FUNCTION  FOR  EACH 
ELEMENT 

[2]  A OF  IN  PANG.  LINKED  ' SEGMENTS'  OF  DIFFEREN 

T FUNCTIONS 

[3]  A CAN  ALSO  BE  CONSIDERED . 

[4]  A INPANG  - VEC . OF  INP.  ANG . VALUES  IN  INTERVAL  0-2 

PI. 

[5]  INPANGLpINPANGl+6 . 28 

[6]  25  GFREG  150  200  850  7000  25  GFSCL  150  200  850  700 
[8]  GFACT  2 50  TRC*-'T' 

[10]  KIM+( 16  0 1 7 , [ 1 ] 5 3 p 1 8 ) , 19  20  0 21  0 22 

[11]  KIM-((6  3p(i5),(3p0),7  8 9 ,(  3p0),10  11  12), 13  14  0,3pl5) 
, KIM 

[12]  LV+~  10  W<-ACC_IA+iOO  UNIT+'U' 

[15]  LB2 : GFCLAO  UARBOUT  10  27  59 

[17]  * YOU  WANT  TO  CONSIDER  THE  INPUT  ANGLE' 

[18]  'FROM  ' , (vLV+LV< 0),'  TO  WHAT  VALUE?' 

[19]  + ( 2 + l\ULC)  x x~'  U'  eUNIT 

[20]  '(YOU' 'RE  GOING  TO  BE  ASKED  ABOUT  THE  UNIT  SUBSEQUENTLY ) 

f 

! 2 1 ] fllMHO  UARBOUT  10O  -*-LAl*\~'U' eUNIT 
.24]  'NOW  SPECIFY  THE  UNIT  BY  PRESSING :' 

L25]  » R FOR  RADIANS' 

[26]  ' D FOR  DEGREE' 

[27]  LA2  :+£41*i (UNIT+WU)  e ' RD' 

[28]  UARBOUT  70  -*-LA 2 

[ 30  ] LA  1 :-*■(  3 + 1 + QLC)  x ( 3 60 , 02  ) lMUV+1  + ' R ' eUNITl  ) vRVZLV+LV<0 

[31]  UARBOUT  70  +£B3 

[33]  IANG+(  ( (LV*MUV)<INPANG)MNPANG<RV*MUV+(  ( (Ol)*l80  ) , 1 )IMUV 
] )/ INPANG 

[34]  UARBOUT  10 
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[35]  ’ I s the  function  constant  for  this  interval  of  the  inp. 
ANGLE?' 

[36]  £B4  :-+LBS*\  1 = ( + /MUVe  ' N ' ) + + / (MUV+U)  e'  Y ' 

[37]  UARBOUT  70  -£B4 

[39]  LB5  :-*LB6*i  ' N ' eMUVO  UARBOUT  10 

[41]  ' WHAT  IS  THAT  CONSTANT? ' O X*-{pIANG) pD 

[43]  £B9 :+LA3*\ {pX  ,l)*pIANG ,10  GFCLG 

[45]  25  PLOTiACC  _IA  ,IANG)  ,[  . 5lW ,X0  GFACT  1 

[47]  ( 250  1 00  , [ 1 . 5 ]3  50  100  ) GFRECT  320  150  ,[l.5]710  150 

[48]  UARBOUT  27  56 

[49]  275  120  GFTXT  'OK'  A 5p  ' '),'UNDO  LAST  PART  OF  GRAPH ' 

[50]  UARBOUT  27  59 

[51]  £44  : -*-£B7  a/(250  100  <CH  P0  S)  , (,C  HPO  S*-2  + GFGCUR)  S3  20  1 50 

[52]  -+LB3  x\A/(350  1 0 0 <CHP0S)  , CHPOSS  7 1 0 150 

[ 53  ] UARBOUT  70  ->-£44  . 

[55]  LB1 :W+W ,X 

[56]  £ ( 8 *RV>  {360  ,&.2S)[_l  + 'R'eUNITl  )+  ' GFCLAO+O' 

[57]  LV*-RVO  ACC _I A+-ACC _I A , I AN  GO  -*LB3 

[60]  £B6  : II*-R*-N-*-0O  T*-Z*-Y*-X*-I  ANGO  VSTK*-H  3p  'IA  ' 

[63]  H*-7  7 00  I*- 0O  GFCLGO  DRWCAL 

[67]  (340  100,[l.5]590  100) GFRECT  530  150,[l.5]710  150 

[68]  UARBOUT  27  56 

[69]  365  120  GFTXT  ' THAT' 1 S IT', (9p*  '), 'RESET' 

[70]  £B2 : GFM OV  200,701-110x11 

[71]  UARBOUT  31  116  8 10  122  8 10  121  8 8 10  93  120 

[72]  LBl-.MUV  GFRECT  (MUV*-(  2 4p(5p0  ),22  44  66  ) + (64  2 p ( 2 2 0 + Jx  5 0 ) , 

630-110xJJ) ) + 64  2 p 5 0 22 
[73  ] GFMOV  5 5 + ,MUVl ;4] 

[74]  UARBOUT  3 1 , , ( " 1 +UA V i VSTK) , 4 4p8  8 8 100  I-I+l 

[76]  £45  :->-LB3*\  a/(  340  1 00  <CH  PO  S)  ,{C  HPO  S+-2  + GFGCUR)  S530  1 50 

[77]  ->-( -LV=~1  ) + LB3*\ a/(  590  1 0 0 ZCHPOS)  ,CHPOS<  7 1 0 150 

[78]  COL+-2  2 0 300  380  460  540  630  640  8 3 0 <CHPOS[  1 ] 

[79]  ^£43xi~i€£,O£-t-BO£A290  370  450  530  630  640  820  900>CHP0S[l 

] 

[80]  ROW*-  430  370  320  310  260  2 50  >CHPOSl  2 ] 

[81]  ->-£4  3xi.  ~leROW*-ROW  *3  8 0 320  310  2 60  2 50  2 0 0 <CHPOSrL  2 ] 

[82]  ERR1*-ERR2*-Q  0 CALC  KIMl  ( ROW  i 1 ) ; COL  \ 1 ] 

[84]  -+LB8*\ERR1=QC>  UARBOUT  10  GFMOV  0.B-B-10 

[87]  UARBOUT  13  27  590  UEMl 1 ; ]0  UARBOUT  27  56 
[90]  ->-£450£43  '.UARBOUT  70  ->-£45 

[93]  £B8  : ->-£4  3 x \ ERR2  =10  11^11*1-15 

[95]  I-(I,0)[l+I=15]O  ->-(£Bl  ,£B2  )[l+J=0] 

V 


V Z+CALCOF 1 iMUVjACTL 

[1]  a CALCULATES  OBJ  .FNC . FOR  CASE  INVOLVING 

[2]  a S. MOMENT  SYNTHESIS 

[3]  AC  TL*-  \ Z*-0 

[4]  FBRIO  F1T_*-F1L+F1I0  F2T+F2  £+ F2JO  F3T+F3  L+  F3IO  F4T-«-F4£+F4r 

[9]  -*-([]£ B[l]  + 2 )xiVF[l]=0 

[10]  +(WORSTlll<~l  + ACTL+ACTL , f/ ( +/F1T*2 )*  . 5 )/0 
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[11]  -*-(DtC[l]  + 2)x1f/F[2]  = 0 

[12]  +(W0RSTl2  ~\<~l  + ACTL+ACTL  , \ / ( + -/F2T*2  )*  . 5 )/0 

[13]  -*-(dLC[l]  + 2)xii/F[3]  = 0 

[14]  ^(V07?5r[3]<~l  + ACTL+ACTL  , T / ( +/F3T*2  )*  . 5 )/0 

[15]  -(□£C[l]  + 2)xlS/F[4]  = 0 

[15]  -*-(i/19i?S2,[4]<  l\  ACTL+ACTL  t\  / (+-fF'\T*2)*  .3) /O 

[17]  FBFTIO  SFT+-SFL  + SFI_ 

[19]  -‘■(□£C[l]  + 2)xii/F[5]  = 0 

[20]  -*(W0RST\_S']<~1*  ACTL+ACTL  t\  / (+JSFT*2)*  . 5)/0 

[21]  ITT+ITL  + ITl^  ~ 

[22]  -•»(□££[  l]  + 2 )xiJ/F[6]  = 0 

[23  ] -+(W0RSTl6]<~l*  ACTL+ACTL  , (I / ITT) ~L / ITT)  /0 

[24]  -(□LC[ l]+2 )xiVF[7]=0 

[25]  -(WORSTZ 7]<  HACTl-ACTL , ( (+/ITT*2 )*pITT)* . 5 )/0 

[26]  -*-(□££[  l]  + 2 )x  iVF[8  ] = 0 

[ 27  ] -►(V0FS.7,[  8 ]<~1\  ACTL+ACTL  ,\  / | ITT)  /O 

[28]  Z^10>  MUV+-WF*  0 

[30]  0F+-+  / (MUV/WF)*(1- ACTL*MUV/ WORST)  *MUV/ EXP 

V 


V Z+CALCOFQ.  ;MUV;  ACTL 

[1]  n CALCULATES  OBJ .FNC . FOR  CASE  INVOLVING 

[2]  n I. TORQUE  AND  S. FORCE  SYNTHESES. 

[3]  r CALCULATES  OBJ.  FUNCTION  FORO  r CASE  OF  S . FO  RC  E+ 1 .TO  RQ 
UE  SYNTHESIS 

[5]  r OF  4 -BAR  LINKAGE ^ ACTL+  \Z«-0O  FBRIO  F1T+F1 L+  FIFO  F2'T+-F2 
L+F2I_ 

[10]  F3T+-F3  L+  F3I0  F4T^F4F+F4I 

[12]  -►(□tC[l]  + 2)xlVF[l]  = 0 

[13]  ->-(WORSTll]<~l  + ACTL+-ACTL  , \ / ( +/F1T*2  )*  . 5 )/0 

[14]  -*•  ( □FF[  l]  + 2)xxyF[2]  = 0 

[15]  +(W0RST12  1<~1*ACTL+ACTL,\’/(+tiF2T_*2)*.  5)/0 

[16]  -►(□££[  l]  + 2 )x1VF[3]  = 0 

[17]  -(J0RSTl3l<~l*  ACTL+ACTL  , T / ( +-/F3T*2  )*  . 5)/0 

[18]  -*■(□££[  l]  + 2 )x  iVF[4]=0 

[19]  -(i/0i?Sr[4]<  1 + AC  TL+-  AC  TL  , T / ( +/F4F*  2 )*  . 5 ) / O 

[20]  F 3 FT  10  FBSMI_ 

[22]  SMT+SM  L+  SMI_ 

[23]  -*>(□££[  l]  + 2 )xiI/F[5]  = 0 

[24]  -*■  ( VORSTL  5 ] < 1 + ACTL+-ACTL,(,{ /SMT)-l/SMT) /O 

[25]  -►(□EC[l]  + 2)xii/F[6]  = 0 

[ 26  ] ->(i/0F5F[6]<-l  + ACTL*- ACTL  , ( ( + /SMT*2  ) * pSMT)  * . 5 ) / O 

[27]  ->-([]FC,[  l]  + 2 )xi!/F[7]  = 0 

[28]  ■+{W0RST\_1'\<~1\ ACTL+ACTL  ,\  / \SMT)/0O  Z«-l 
[30]  MUV+WF* 0 

[ 31  ] 0F<-+  / (MUV/JF)*  ( l--4C’2,tT^i/7/V0F5r)*Wi/7/FJP 


V 
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V Z+COMMON 2 INT 

[1]  n FINDS  INTERSECTION  OF  INTERVALS  INT Cl  2] 

[2]  n AND  INTl 3 4]. 

[3]  n Z = 0 IF  rtfBBB  IS  NO  INTERSECTION . 

[4]  n IF  INTERSECTION  IS  A POINT,  Z[l]=Z[23. 

[5]  n INTlll=INTl2]  OR  INTl 3 ]=INTZ 4 ] IS  OK. 

[6]  -*■(  ( </INTlA-  l])v</INTl2  33)  + Z«-0 

[7]  Z*-INTl{bINT)l2  3]] 

V 


V Z+COMMON 3 INT; MUV 

[1]  n FINDS  INTERSECTION  OF  INTERVALS  INT [1  2], 

[2]  n INTl  3 4]  4/VB  I/1/:T[5  6]. 

[3]  n SBB  COMMENTS  IN  C0MM0N2 

[4]  + {0  = ppMUV+-COMMON2  4 + IA/T)  + Z«-0 

[5]  +{0=ppMUV+COMMON2  MUV  ,~2+INT)  + 0 
[ 6 ] Z«-MB7 

V 


V DASHES; ABC; MUV 

[1]  n CREATES  BUFFERS  FOR  DASHED  LINES 

[2]  n Z/SBB  BY  PB021  4//B  PL0T3 . 

[ 3 3 TRC-*-  * 4 1 

[43  GVRDT+-1  2 3 4,4  12pl  ~12  + G7i?£r 

[53  GF4/7T  20  4B£>2  50  DASHTYPE+3 

[83  LI :{320 , MUV) GFD ALINE  9 8 0 , MUV+-1 5 0 - ABC 

[93  -*-£1  x i 300*  ABC-*-  ABC+2  5 

[10  3 GFACT  30  4BO-2  5 

[123  £2  : (/ Vi/7, 400  ) GFD ALINE {MUV+21 0+ ABC)  ,40 
[133  +L2 xi 300 *4B£>4BC+ 2 5 

[14  3 GFACT  40  4B£«-4  4 

[163  L3  : ( MUV  ,585  ) GFD  ALINE  {MUV-*-^  50  +ABC ) ,195 

[173  ->L3x;  52  8*4B£>4BC’+44 

[18  3 B2+GFRBUF  2 

[193  B3+GFRBUF  3 

[203  B^*-GFRBUF  4 

[213  GFPBUF  2 3 40  TRC+' AT * 

V 


V DRWCAL;  A;B;H; VI ;V2 \W3 ; & ;HD ; LL; UR ; VD; MUV1 ;MUV2 
[13  n BB4V5  4 CALCULATOR  IN  THE  ACTIVE  REGION 
[23  n TO  BE  USED  BY  'CALCULATOR1 . 

[3  3 A*- 2 2 00  B*- 2 0 00  5«-5  0O  V1^7  0 

[73  V2-*-9  0O  W3-*-2*W20  A«-10 

[103  HD+-W 1+A 

[113  ££«-($4  2 p 4 , B ) + 7 £-*-2  4p(5pO),(B+A)x;3 

[123  UR-*-  (64  2p  (A  + Wl ) ,B  + H)  + VD 

[133  ££-^££,(6 4 2 p ( A + HD ) ,B)  + VD 

[143  UR+UR,  (6)4  2 p ( A + HD+W 1)  ,B  + H)  + VD 
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[15]  11- LL  , (04  2p  (,A  + 2*HD)  ,B)  + VD 

[16]  UR-URAS  4 2p (A+W1+2*HD)  ,B+H)  + VD 

[17]  LL-LL,  <.MUV1-A  + 3*HD)  ,B 

[18]  UR* UR , ( 4+Wl+3*HD) ,B+A+2xff 

[19]  LL-LL  A^  2 2f>MUVl-MUVltB  + 2*a+L)+MUV2-2  2p(3pO),S+A 

[20]  UR-URAS  2 2pMUVl-MUVl+Wl,H)+MUV2 

[21]  LL-LL, MUVl-MUVl+lO  UR-UR ,MUV1+W2 ,H 

[23]  LL-LL, MUV1-MUV1+(W2+A)  ,0^  UR-UR .MUV1-MUV1+W3 ,H 

[25]  Ll-LLtHUV2-(A+'l*HD)  ,B 

[26]  UR-UR ,MUV2-MUV2+(W2+A+W3  ) ,H  + 2*H+A 

[27]  CZ>LL,(04  2 p/tfSKl  «-(  AfffPl  [ 1 ] + A ) , S ) + KB 

[28]  Stf«-Stf,(04  2pMZ/Kl+Vl,tf)  + KB 

[29]  LL-LL, MUVl-(A  ,B)-2*A 

[30]  UR-UR ,MUV2-NUV2+(V1 ,H) +3*1 

[31]  LL-LL, MUV 1-10  UR-UR ,MUV 2+1 
[33]  LL  GFRECT  UR S>  lMBBSSI  2 7 5 6 

[35]  ( 25  2 0 + /4  , B ) GFTXT  ' LN' , (3p»  ' ) ,' LOG'  A 3p’  • ) , » 1 / I’  , ( 1 4p  ’ 
'),' KEYBOARD' ,( 7p»  '),'Y+X' 

[36]  (A/Z/Kl-20  1 0 + ^ , B ) GFTX  I ( 1 8 p ' ' ) ,' R' 

[37]  (VSKl^O  20+MSK1  )GFTXT(  18p  * M.’B’ 

[38]  (MSKl«-0  20+MUV1  )GFTXT(  18p  ’ ' ) ,'T' 

[39]  (MUV 1-0  2 0 +MUV 1 ) GFTXT ( 1 8 p ’ '),'N' 

[40]  (MUV1-0  2 0 +MUV 1 ) GFTXT ( 1 8 p ’ '),'EA 

[41]  (MUVl-25  1 8 + /1  , B+  KB«-S+ A ) GFTXT  { 5 p ' ,)>,10f,(4p'  M.T'.Cll 

p'  '),  'NUMBER  THRU  THE'  ,{ 4p  • '),'Y~X' 

[42]  0111111111  SFIMBS(610  2p-VSKl)  + 2 lOplO  6 0 ~4  ~4 

0 6 10  10  “4  4 0 0 4 12  16  16  12  9 9 

[43]  (MUV2-2  lH+MUVl)GFTXT(lp'  ' ) A X'  A 5p  ' '),'X',(Ap'  • ) ,» X' 

[44]  (MUV  l«-2  4 m + A,B  + 2xVD)  GFTXT  'X',(Hp'  '),' SIN ' ,(2p ' ')ACO 
5',(3p'  1 ) , 1 I4B  ' , ( 3p  ' ' ) ACLICK  HERE  TO  INPUT' A 2p ’ *),' 
Y*X' 

[45]  (0  15+MUV1  )GFTXT  ' 2',(5p’  'J.’-l'.Op'  ,),,-l,,(4p’  •), 

. _!» 

[46]  (MUV2-MUV1+0 ,VD) GFTXT  ' X'  A 4p T ’ ) , ’ SIB ' , ( 2p » M.'SSS’.O 

p*  • ) , 'I4B‘  , (3p  • * ) , ’ UNDO ' 

[47]  (.VSK1+0  , KB)  SFIII(43p  ' ' ) A Y*X' 

[48]  {NUV1-MUV2-*  5 ) SFIXI ( 2 9p ' ' ) , » ( IN  RADIANS) ' 

[49]  (MUV 1+0  20 )£FIXI(29p ’ ') A I A- IN P . ANGLE' 

[50]  0 111  SFK4BS(64  2pMUVl)+2  4p  1 5 1 0 1 25  20  0 30  30 

[51]  UARBOUT  27  59 
V 


V DRWFB ; S2  1 ; S3  2 ; S4  3 ; II 1 ; II 2 ; MUV  ; 12 1 ; 12  2 

[1]  n DRAWS  A FOUR- BAR  IN  THE  ACTIVE  REGION. 

[2]  ( Sl-*-7  5 0 10  0 ) GFC  IRC  8 40 

[3]  (S2^800  190) GFC IRC  8 40 

[4]  (B3«-9  6 0 2 3 0 ) GFC  IRC  8 40 

[5]  (S4*-970  10  0 ) GFC  IRC  8 40 

[6]  S21^3*S21*(+/(S21^S2-S1 )*2)*. 5 

[7]  (C1+C21) GFLIN E C2-C21 

[ 8 ] S3  2 -'-8  xS32r(  + /(S3  2^S3  -S2  )*  2 )*  . 5 
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[9]  (C2+C32 )GFLINE  C3 -C32 

[10]  C43-«-8xC,43-S-(+  / ( 04  3-*-(?4  -S3  )*2)*.  5 

[11]  (S3+S43 )GFLINE  S4-S43 

[12]  71 1«-L  . 5 + Cl+8xrilf  ( +/(Tll«-73  5 8 0 -Cl  ) * 2 ) * . 5 

[13]  T12+1 . 5+S1+8 xri2r( + /(T12«-765  80 -SI )* 2 )* . 5 

[14]  0 111  GFVABSS 4 2p711,735  80  765  80,712 

[15]  725  80  GFLINE  775  80 

[16]  MUV*- 725  70  730  80  740  70  745  80  755  70  760  80  770  70  775 

80 

[17]  01010101  GFVABStS)  8 2pM[W 

[18]  721H  . 5 + C,4  + 8x2’21t(  +/(721^955  80-S4)*2)*.5 

[19]  722H  . 5+S4  + 8xT22*( +/(722«-985  8 0 -S4  ) * 2 ) * . 5 

[20]  0 111  GFVABSS) 4 2p721,955  80  985  80,722 

[21]  945  80  GFLINE  995  80 

[22]  MUV*-  945  70  950  80  960  70  965  80  975  70  980  80  990  70  995 

80 

[23]  01010101  GFVABS^B  2p MUV 

[24]  750  100  GFARC  40  "25  75  60 

[25]  101  GFVINCS) 3 2p2  10  '2  10  9 3 

[26]  Oil  GFVABSS 3 2p700  140  700  50  785  50 

[27]  MUV*- 696  132  700  140  704  132  777  54  785  50  777  46 

[28]  011011  GFVABSS 6 2p MUV 

[29]  MUV*- 720  40  725  50  735  40  740  50  750  40  755  50  755  40  770 

50 

[30]  01010101  GFVABSS 8 2p MUV 

[31]  650  0 GFRECT  1020  280 
V 


V DRWFB2 

[1]  a DRAWS  ON  SCREEN  THE  4- BAR 

[2]  r STORED  IN  BUFFER  SI. 

[3]  UARBOUT  SI 

[4]  GF ACT  10  FB AXES'}  FBTXT 
[7]  SFMOF  0 7790  CMFBSST  13 

V 


V DRW STAR  X ; STAR 

[1]  fl  SS4VS  STARS  FOR  PRECISION  PO IN TS ( S YN THE SI S) 

[2]  -»-(  0 = pJT) /0 

[3]  +(l  = pp^)/DLC[l]  + 20 

[5]  STAR*- 8 2 p 0 ~3  0 3 3030  "3  ~3  333  3 33 

[6]  ( 8 p 0 1 )GFVABS*{  8 2p2iX)+STAR 

[7]  -([]SS[  1 ]-l  )xO*pX^2 


V 
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V EXPIN_\  KK  ;NJ  ;I 

[1J  A ALLOWS  USER  TO  SPECIFY  EXPONENT  OF  SCORE  CURVE 

[2]  GFCLA 

[3]  f THE  SCORE  SIGNIFYING  PERFORMANCE  OF  THE  LINKAGE  RELATIV 
E TO  ' 

[4]  'A  DESIRABLE  CONDITION  ( FOR  EX:  MAX  MAG  OF  4 BEARING  REA 
CTION) ' 

[5]  ' WILL  BE  ASSIGNED  THROUGH  THE  FORMULA :» 

[6]  i 00  6 0 p '~'0  i 0 

[9]  ’ SCORE  = (1  -(.ACTUAL  V ALUEiWORST  VALUE) )*EXP' 

[10]  tOO  6 0 p ' ~ ' 

[12]  ' WHERE  MEANS  EXPONENTIATION . ' 

[13]  'NEXT  YOU' 'LL  BE  PROVIDED  WITH  A MEANS  FOR  CHOOSING  AN' 

[14]  'ADEQUATE  VALUE  OF  ''EXP''  FOR  EACH  DESIRABLE  CONDITI 

ON  . ' 

[15]  lOO  'PRESS  "RETURN''  TO  CONTINUE .'  0 KK+-\ 3 

[18]  EXP<-\0r>  NJ+ltpJOINTSQ  I+-1 

[21]  EXP-*- EXP , SC  0 REPLOT  3 lOp  'MAX  MAG  OFFORCE  AT  JOINT  '.(?!) 

I f 

* 

[22]  -+tULCl  1]-1  )*\NJ>I+I+1 

[23]  -*LB1*\SYNCASE='  2 ' 

[24]  EXP-*- EXP , SCO  REPLOT  3 lOp'Afyi*  MAG  OFSHAKING  FORCE  ' 

[25]  EXP*- EXP , SCO  REPLOT  3 9p 'PTP  VALUEOF  INPUT  TORQUE  ' 

[26]  EXP+EXP .SCORE PLOT  3 9p 'RMS  VALUEOF  INPUT  TORQUE  ' 

[27]  EXP-*- EXP  , SC  0 REP  LOT  3 10p' MAX  ABS  VALUE  OF  INP  TORQUE' 

[28]  ->0 

[29]  LB1  :EXP-*-EXP  ,SC  OREPLOT  3 lOp  'PTP  VALUE  OF  SH  AKINGMOMENT 

» 

[30]  EXP+EXP , SCO  REPLOT  3 lOp 'RMS  VALUE  OF  SH AKINGMOMENT  ' 

[31]  EXP+EXP .SCORE PLOT  3 Zp'MAX  ABS  VALUE  OFS .MOMENT' 

N 


V FBAXES;  A ',  I ; FL  AG ; MUV ; MUV2  ‘,X  AXIS',  Y AXIS;  JOINT',  S 

[1]  n DRAWS  LOCAL  AXES  OF  4 -BAR 

[2]  TRC+'T'O  4 -*-3  2 p ( C2 -Cl  ) , (C3 -C2  ) ,<73 -C4 
[4]  I«-lO  LINKS-*- 3 2p’ 122343' 

[6]  LI  iMUV2'*-30*MUV-*-MUV*  ( + / {MUV*-  ,AlI  i ] )*2)*.  5 

[7]  XAXIS+MUV  2,4>(150  ROTVEC  S)  , ( -MUV ) . [ 1 . 1 1MUV-~1  50  ROTVEC  S 
+-6*MUV 

[8]  YAXIS-90  ROTVEC  XAXIS 

[9]  GFMOV  JOINT-*-*'  C ' , LINKS  L I ',  FLAG*-l  + l=  t'  FLAG'  , LINKSlI ',  ]] 

[10]  110  1 GFVINCl ( ( pXAXIS) p . 5 ) + X AX I S* 1 llFLAGl 

[11]  GFMOV  JOINT 

[12]  1101  GFVINCl ( ( pYAXIS) p . 5 )+YAXIS*l  ~1 iFLAGl 

[13]  -*£1*1  4>I«-H-1 

[14]  LIN  KS*-  LINKS  , ( FL  AG1  2 , FL  AG2  3 , FL  4(74  3 )<pLINKS 

[15]  T3Z>( (FLAG12*1 80 ,412 ) , [ 1 ] ( FLAG2 3 xi 80 , 423 ) , [ . 5 1FLAGA 3 * l 3 0 
,443  ) , 3p0 

[16]  GFMOV  0 5000  UARBOUT  13 


V 
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V FBDVA  TYPE;  DEN;  P24  \D\  COSA  ; A ; P43  ; EL\ 02  3 ;043  ;.¥23  ;P43  ;UTX 
Cl]  n DISP./VEL./ACC.  ANALYSIS  OF  4 - BAR 

[2]  UTX*-'ERR*- 10*0* 

[3]  EL*- pPH  1 12*-(  0,  \ N ) * (.  02  ) * N 

[4]  U2*-U1  + , 1 04i?12"*-/412x(  20PHI12  ) , [ . 5]loPffJ12 

[5]  R2H*-{  (ELpUH)  , [.  5]P0p74)-P12+$(P£,2)pPl,71 

[6]  COSL*-(D+(  -443*2  )+423*2  )f  2 x 4 2 3 x ( p<-+ /P2  4*  2 )*  . 5 

[7]  H*-COS&  ARG  TYPE*  ( l-POSA*2  )*  . 5 

[8]  PHI23*-PHI23-(02  )x(02)ZPHI23+tL+R2Hll;)ARG  P24[2;] 

[9]  P4  3-*-(  423x  ( 20PHI23  ) , [ . 5 ] 1 OPHI2  3 ) - R2  4 

[10]  PPJ4  3-*-P4  3 [ 1 ; ] ARG  P43[2;]0  72«-71  + P1  2 [ 2 ; ] 

[12]  U3*-U1  + ,~1  0\V3*-R12  + R2  4 +P4  30  V3*-V1  + ,1  0 + 73 
[14]  U2D*-~  A 12*W  I*  lOP  HI  12^  V2D*-A12xWI*2oPHI12 

[16]  U2DD*-  - A1 2X(VJ*2)X2  OPHI 120  V2  DD*- -A  1 2*  ( VI*  2 ) x loPffll  2 

[18]  PHI  2 3 D*-WI  x 02  3-«-(  4 1 2 xIOPHI  1 2 - PHIH3  ) + 4 2 3 x DEN*-loPHIH3  - PHI  2 3 

[19]  PHIH3D*-WIxG<A3*-(  A12xloPH 1 12 -PH  123)  iAH3x DEN 

[20]  U3D*-~  AH3  xpHIH3Dx  lOPff  J430  73  P«-44  3 xp#I432?x  20PPI4  3 

[22]  H23*-{  ( -443xG43*2)+(412x20PPI43-PflI12)+423x(G23*2)x20PPI4 
3 -PHI 23  ) ■*•  .42  3 xJ)EN 

[23]  P4  3-<-(  (.423xG2  3*2)  + (412x20PPI23-PffI12)-443x(G43*2)x20PffI2  3 
-PHIH3 ) t 44  3 xDEN 

[24]  PHI23DD*-H23xWI*2C>  PHIH3DD*-HH3  xWI*2 

[26]  U3DD*-~  AH  3 x ( PH IH3  DDx  10PHIH3  ) + ( Pff  J4  3 D*2)x  20PHIH3 

[ 27  ] 73  0P-«-44  3 x ( PH I H3  DDx 20PH I H3  ) - ( PHIH3D*  2 ) x lop# J43 

V 


V FBDYN ; MXY 23 G; VEC 

[1]  n DYNAMIC  4.V4L  . OP  4-54P  POP  GIVEN  MASS  PARAMETERS 

[2]  MXY23  G*-PH  123  R0TVEC3<A  ( ( 1 ) , 2 ) p ( *02  3 , 702  3 ) x 4 2 3 

[3]  7PP«-(  ( - J MXY  23  Gx  V2  DD , [ . 5 ] U2DD ) + PH  123  DDx  K02  3x423*2)  , [ . 5]Pff 
IH3  DDx K0H3 x 44  3*  2 

[4]  F3I*-  ( (7PO[l;]xP2  2)-7PO[2;]x012)  * DETC 

[5]  F3I*-F3I , [ . 5 ] ( (011x7  PO  [2  ; ] )-P21x7PP[  1 ; ] )*  DETC 

[6]  P4r-^P3I+(  (J04  3xJ.V44  3 ) -704  3 xDANIH3  ) ,[  . 5 ] ( *04  3 x DAN  14  3 ) +704 
3 x IN 44  3 

[7]  F2I+F3I- (M23xU2DD ,1 . 5 ] 72 DD)+( (X 023 x IN A23 ) - Y023 x DAN 123 ) ,[ 
. 5](X023xDANI23)+Y023xINA23 

[ 8 ] PlI«-(  ( (X012xINA12)-Y0l2xDANI12  ) , [ . 5 ] ( *01  2 x DAN II 2 ) + Y01  2 x I 
NA12)-F2I 

[9]  iri^412x-/(  ( 10PPJ12  ) , [ . 5]20PPI12  )xP2I0>  SFI+-F1I+FHI 

V 


V FBFTI 

[1]  r CALCULATES  S . FO RC E+ IN P . TO R . DUE  TO  OPTIMUM  M.P. 
[2  ] ITT*- A 1 2 x - / ( ( loPHI  12)  , [ . 5 120PPI  1 2 ) x P2 1 

[3]  SFI_*-~  F1I+FHI_ 

V 
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[1] 

[2] 

[3] 

[4] 

[5] 
[8] 
[12] 

[13] 

[14] 

[15] 
[IS] 
[19] 
[22] 
[23] 
[ 24] 


[27] 

[28] 
[29] 

[32] 

[33] 

[34] 

[35] 


FBGEO  Z\MUV 

A ASKS  FOR  4 -BAR  GEOMETRY 
-►LSI  x i ~le  Z 

1 ENTER  COORDINATES  OF  PIVOTS  1 AND  4 (U 1 VI  U 4 74 ) : » 
"►(□LC[l]  + 4)xl5  = pi,  MUV+U 
UARBOUT  7 0 ' TRY  AGAIN '0  +ULCZ  l]-3 
Ul+MUVZl]$  V1+MUVL210  S4-.VSF[3]0  74<-MS7[4] 

LSI  : -*-LB2  * i ~2e  Z 

• ENTER  LENGTHS  OF  CRANK,  COUPLER,  AND  LINK  43* 

’ (LI 2 L2  3 L43 ) : ’ 

"►(nLL’[  1 ] + 4 ) x i 4 = pl  ,MUV-*~U 

UARBOUT  7 0 ' TRY  AG  AIN  ' 0 -*£|LC[l]-3 

A12+MUV\_  1 ] 0 A2S-MUVI  2]0  A43«-Aftf7[  3 ] 

LB 2 :->(~3<£Z)  + 0 
MUV+-0 

LI: 'WHAT  IS  THE  TYPE  OF  DYAD  234?’ 


(9x,VS7)4’  3)  HELP » 

-*■(  □LL’[  1 ] + 3 )x  x ( 0=ppIIPF)  aa/(7YPF«-D)  €(  -MUV)*  * 12  3* 
□4PS0SI  7 0 ->DLC[l]-2  0 -{~MUV+TYPE='  3’  ) + 0 
* the  rips  is  +1  if  when  you  stand  at  joint  2 and 1 

’LOOK  AT  JOINT  4,  JOINT  3 IS  AT  YOUR  LEFT.'1 
’ OTHERWISE  THE  TYPE  IS  - 1 . ’ 

’ ’ 0 -L10L2  '.SPEC 


V F BRI_\  MXY 23  G\  VEC 

[1]  r CALCULATES  REACTIONS  DUE  TO  M.P.  FROM  OPTIMIZATION 

[2]  MXY2  3 G-*~PHI2  3 ROTV  ECZ<*  ( ( N + 1 ) , 2 ) p U02  3 , 102  3 ) x A2  3 

[3]  7£C«-(  ( --/MXY23G*V2DD , [ . 5 ] U2DD ) +PHI2  3 DD*  K02  3 *A  2 3*  2 ) , [ . 5]PS 
J43SSxXS43  x443*2 

[4]  F 3I«-(  ( 7FS[  1 ; ]xp22  )-7F<?[  2 ; ]xS12  )*DETC 

[5]  F3I+-F3I , [ . 5](  (Sllx7SS[2  ; ] )-C21*VECl  1 ;]  )t  DETC 

[6]  F4I«-F3I+(  (X04  3xJJ|M4  3 ) - Y04  3 *DANI'*3  ) , [ . 5 ] ( XS4  3 x£  4//I4  3 ) + Y0  4 
3 x IN  A 4 3 

[ 7 ] F2I_+-F3I_-{M023*U2DD , [ . 5]72SL>)  + ( U02  3 *IN  A 2 3 ) - Y02  3 *DANI2  3 ) , 

[ . 5 ] (X02  3xDANI23  ) +Y  02  3*  IN  A 2 3 

[ 8 ] FlI«-(  ( (X012*IN  A12  ) - Y012*D  AN  112  ) , [ . 5 ] ( X01  2 *DANI1  2 ) +Y01  2 xj 

V412 )-F2I 

V 


V FB SMI 

[11  r CALC.  S. MOMENT  DUE  TO  GIVEN  M.P.'S 

[2]  SMI+-(  (4>$F4r)-  . x(S4,t/4  )-RPOINT)  + { (4>$F1I)-  . x ( Si , 71  ) - RPOIN 
T) + IT I 

V 
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7 FBSML 

[1]  n CALC.  S. MOMENT  DUE  TO  EXT.  LOADS 

[2]  SM  L*-  - ( (^PtJ3  )-  . x(04,74  )-RPOINT)  + i (4><$FB_  12)-  . x ( 01 , 71  ) - /? 
POINT) + ITL 

7 


7 FBSMI 

[1]  n CALC.  S. MOMENT  DUE  TO  M.P.  FROM  OPTIMIZATION 
[ 2 ] SMI*-  - ( ( 4>$F4I)  - . x ( Z/4  , 74  )-RPOINT ) + ( ( <t*J}FlI)  - . x ( 01 , 71  )-RPOIN 

T)  + 1 TI_ 

7 


7 FBSPEC 

[1]  a OSFB  FOB  5PF0.  OF  4-B4B 

[2]  GFCLG5  GFACT  10  DRWFB 

[5]  Bl+GFRBUF  1 0 GFPBUF  10  FBTXT 
[8]  FBGEO  1 2 30  SPFO 

[10]  FB4O12fFB4O23fF£4O43f0 

[11]  LBl'.GFCLGU  UARBOUT  BlO  FBTXT 
[14]  FBAXESC>  TRC*-' AT'O  UARBOUT  13 

[17]  ' 4BF  FOB  £0040  COORD.  SYSTEMS  SHOWN  OK ?' 

[18]  i , 5t  (“5  + 5x»tfY*  x 1 + [5)(|>»LB2  OOGOODBC+l' 

[19]  CMBBOOT  70  -*-l1£0[1]-2 

[21]  LB2  : ’ WHAT  COORD.  SYSTEMS  MUST  BE  CHANGED ?' 

[22]  » 1)  LINK  12* 

[23]  ’ 2)  LINK  23' 

[24]  * 3)  LINK  43  » 

[25]  * FOR  EXAMPLE,  YOU  MAY  ENTER:  2 31 

[26]  -*-£B3  x t (a  / MUV  e 1 2 3)a3>p  ,MUV-U 

[27]  UARBOUT  70  ' TRY  AGAIN1 0 f[]L0[l]-3 

[30]  £B3  : F£  401  2-<- i ( (leA/B7)  + ,~'  ) , 'FB4012’ 

[31]  FLAG2  3-*-i  ( (2etf07)+  ) , 'FLAG2Z' 

[32]  F£4043«-i(  ( 3eMUV)+  'FB4043  ’O  -*-£Bl 

7 


7 FBSTATIC ;F120;F230;F430;7FO 

[1]  a 4-B4B  STATIC  ANALYSIS 

[2]  F120-*-(PffI12xl80*Ol  )R0TVEC2l  14+  ( 14x1=  1 + p F 1 2 ) ( ’ F 1 2 ’ , 1 1 p ’ 

»),’$(  (0  + 1 ) , 2 ) p FI 2 ’ 

[3]  F2  3 G*-  ( PH  12  3x180tO1)  ROTV  F02fl4+(14xl=  1 + p F2  3 ) 4>(  ’ F2  3 ’ , 1 1 p f 

• ) , ’$(  OV+1 ) ,2)pF23»  _ 

[4]  F4  3 G*~  ( PBJ4  3 xl  80  vOl  ) ROTV  EC2&  14  + ( 14x1=  1 + p F4  3 ) 0 ( ’ F4  3 , 1 lp 

•),'$(  U+l  ) ,2)pF43» 

[5]  VEC+-  ( ( - / F2  3 B[  2 1;]xF23)-T23),[.5](-/F430[2  1 ; ] x F43  ) - T43 

[6]  DETC*-{  (01 1-*-- 423xloPBI23  ) x 02  2 + -44  3 x2oPBI4  3 ) - ( 01  2f  4 2 3 x 2 oPtf 
12  3 )x021f 443xloPBI43 

[7]  FL 23  f ( (7F0[ 1 ; ]x022 )-7F0[ 2 ; 1x012 )*DETC 

[ 8 ] F2L*-FL2  2_*-FL2  3_ , [ . 5 ] ( ( 01 1 x7F0[  2 ; ] ) - 02  1 x 7F0[  1 ; ] )iDETC 
[9]  F4£f F£  43-*-F£2  3_-F43 
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[10]  F2L—FL1  2_—FL2  3_+F2  3 

[11]  F1L+FL_12+-FL12_+F12 

[ 12 ] ITL-(A 12 x( -FL1 2_[ 2 ; ] *2oPH 1 12 ) + { F LI 2_[ 1 ; ]xioPtfJ12 ) ) + ( -/£1 
2 6[  2 l;]xF12)-T12 

[13]  SF£«~FL_12+F&_43 

[14]  J0INTS-LINKSll;3  4 ] , [ 1 ] LINKSl  ; 3 4] 

V 


V FBTOTAL 

[1]  n CALC.  TOTAL  REACTIONS,  S. FORCE  AND  I. TOR 

[2]  F1T-F1L  + Fllf>  F2T-F2L+F2IO  F3T-F3L+F3I 

[5]  F4r-«-F4£+F4rO  SFT-SFL+SFI 0 ITT-ITL+ITI 

V 


V FBTXTi C2 1 ;<?3  2 ;C43 ;T11 ; T1 2 \MUV\T2A \T 22 

[1]  r CFF/ITFS  TEXT  IN  4 -BAR  PICTURE 

[2]  740  112  ffFrxr  » 1* 

[3]  730  195  GFTXT  ’2' 

[4]  945  240  GFTXT  ’3* 

[5]  985  110  GFTXT  '4’ 

[6  ] 795  134  GFTXT  ' . 1 

[7]  795  120  GFTXT  ' 4>  =CONST.' 

[8]  803  115  GFTXT  '12* 

[9]  785  30  GFTXT  'U' 

[10]  690  38  GFTXT  'O' 

[11]  688  135  GFTXT  'V' 

[12]  GFMOV  0 7790  UARBOUT  13 

V 


V FBDYN A ID \MUV\ S2 3 ; C23 ; 543  ;F4 3 

[1]  r GENERATES  VARS.  FOR  FBDYN 
[2  ] INA12-{  MUV*--A12*VI*2  )*2oPHI12 

[3]  DANI12+MUV*loPHI12 

[4]  IN A23--A23 x(PHI 23 DD*S23-loPH 1 23  ) + ( M UV-PHI2  3 D*  2 ) * C2  3-2QPH 
123 

[ 5 ] DAN  123- A2 3 *(PHI23DD*C2 3 ) -52  3 *MUV 

[ 6 ] T/M  4 3 - 4 4 3 x ( PH  1^3  DD*  S4  3+-1  OFF  J4  3 ) + ( ^UV-PHI^3  D*2  ) * C4  3-20PH 

1 43 

[7]  DANIli3—Al\  3 x ( PH  1^3  DD*  C4  3 ) - 543  *MUV 

[8]  X012-M12*X12iA120  Y01  2-MI  2 *Y  1 2 4 A.1 2 

[10]  X023-M23*X23iA230  Y023-M23*Y23 i A23 

[ 12  ] K02  3— M2  3 x ( {X2  3*  2 ) + ( Y2  3*  2 ) +K2  3*  2 ) 4 A2  3*  2 

[13]  X04  3^43x14  3 4/14  30  Y04  3*-.V4  3 xJ43  4 44  3 

[15]  F<?4  3—M'4  3 x ( (Y4  3*2)+(Y43*2)+F4  3*2)4. 44  3*2 

V 


V Z*-FBITCOEF\  272  23  ;D323  ;£>4  4 3 \CO\IN 

[1]  r CALC.  COEFS.  FOR  EQUALITY  CONSTRAINTS 

[2]  r FROM  I. TORQUE  SYNTHESIS 

[3  ] PRC PTSL+PH 1121  IN*-l+l  . 5 +PRCPTSLZ  1 ; ]*3  60* if]  , [ . 5 ] PBC PTS £ C 2 ; 

] 

[4]  PRCPTSLD*-(PRCPTSLZ  1;]*180tO1)  ,[.  5 1PRCPTSLZ  2 ; ] 

[5]  D 223*-((  (U2D*U3DD)  + V2DxV3DD)  + (U3D*U2DD)  + V3Dx  V2  DD)*WI 

[6]  £>3  2 3-*-(  ( ( U3D*V2DD)-V3D*U2DD)-(U2D*V3DD)-V2D*U3DD)iWI 

[7]  DW3-*-{{U3D*U3DD)  + V3D*V3DD)*WI 

[8]  C0<-(PRCPTSLl2;l-ITLLlN])®D223lINl  ,D3  2 3 ZIN_  ] , [ 1 . 5 ]£443  [Iff] 

[9]  NEC  .CONDITIONS  ARE  GOING  TO  BE  CHECKED 

[10]  Z^O 

[11]  -*■(  ( (^02  3WI/7-A:02  3«/1Z)<C,C>[1])^£<7[1]^Z(92  3W/5J-^02  3A/I/V)/0 

[12]  -*■(  (Y023MINSC0Z2])*C0Z2]ZY023MAX) /0 

[ 13  ] -►(  ( (K023MIN+K0^3MIN)<C0l3  ] ) *COl  3 ]ZK02  3MAX+KO^  3MA  X ) /0 

[14]  KITC 0*~IT_C O-y  ITCO*-CO 
[16]  KCURVE+CURVE 

[ 17  ] CURVE*-ITT*-ITL  + ITI*-(  ITCOi  1 ]*  £2  2 3 )+(ITCOZ  2 ]*£>3  2 3 ) + ITCOl  3 ]* 
£>4  4 3 

[18]  KPRCPTS+PRCPTSLD'U  KFIRST+FIRSTCr  FIRST*-!*- 1 
[21]  KRMS+- RMS'D  KPTP*-PTP*>  KMIN+MINC  KM  AX*- MAX 

V 


[1] 

[2] 

[3] 

[6] 

[8] 


FBKO UT ; MUV ; S ; SEL 
R USED  FOR  PLOTTING 


KINEMATIC  PARAMETERS 


R CALCULATED  BY  FBDVA. 
GFCLGy  UARBOUT  BlO  FBAXES 
TRC*-'  AT'  0 FBTXT 
PLEASE  ENTER  SELECTION 


[9] 

' 1) 

POS. 

OF  JOINT  2' 

[10] 

' 2 ) 

VEL  . 

OF  JOINT  2' 

[11  ] 

* 3 ) 

ACC . 

OF  JOINT  2’ 

[12] 

1 4) 

POS. 

OF  JOINT  3' 

[13] 

’ 5) 

VEL. 

OF  JOINT  3’ 

[14] 

’ 6) 

ACC . 

OF  JOINT  3* 

[15] 

» 7 ) 

ANG . 

POS.  OF  LINK 

’ ,MUV-{ 1=FLAG23)$'  23’ 

[16] 

' 8 ) 

ANG . 

VEL.  OF  LINK 

’ ,MUV 

[17] 

* 9) 

ANG. 

ACC.  OF  LINK 

' ,MUV 

[18] 

1 10  ) 

ANG  . 

POS.  OF  LINK 

' ,MUV*-(1  = FLAG^3  )4>f43» 

[19] 

* 11  ) 

ANG. 

VEL.  OF  LINK 

’ ,MUV 

[20] 

’ 12  ) 

ANG. 

ACC.  OF  LINK 

' ,MUV 

[21] 

-*•(□££[  l]  + 4 ) 

*i  ((0  = ppS)AA/(S«-[])€il2) 

[22] 

QARBOU 

T 7 0 

'TRY  AGAIN ' 0 

-+ULCI  1 ]-3 

[25] 

SEL*-  ( » 

U 2.  [ . 

l]72’,(4pt  • 

)),[!]( 'U2D,l. l]V2D' , ( 

] ’ U2DD 

, [ . 1 ] V2DD' 

[26]  SEL+SEL, [l](  » U3 , [ . 1 ] K3  * , ( 4p  * 
')).[. l]f U3DD  ,Z . 11V3 DD' 

[27]  SEL+SEL, lll(  'PHI23'  , (8p»  ')) 

PHI23DD' ,6p'  ' 

[28]  SEL*-SEL  ,Lll(  'PHIH3'  , ( 8 p * •)) 

PHIH3DD'  , 6p  » ' 


* )).[!]( ' U3D ,Z 
[ 1 ] ( ’ PH 123  D' 
[1](  ' PH IH3  D' 


' ')),[. 

1]73£' , (2p ' 


( 7 p * 
( 7 p » 


' )).[.1]» 
')),[.  1]» 
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[29]  GFCLO  MUV+- ' [ 1 , ( ( 6 <5)  / ' . ' ) , ' 1 ] ' 

[31]  SEL+'PHI  12,  * ,MUV  , , SELLS  \ ] 

[32]  FNC+{-3<S)±'PL0T3'O  TRC+-  tAFLAG\.S*  ] 

[34]  £’  ( 10x5)  » ,F/VC,  ' ' ,SELO  AFL  AGl  S ; ]■*-'  T ' 

V 

V FBMPOLIM;MUV 

[1]  n 04L0.  LIMITS  FOR  LUMPED  MASS  PAR  AMS. 

[2]  X012MAX*-(MUV+M12MIN  ,M12MAX)[1+X12MAX>0]*X12MAXtA12 

[ 3 ] X012MIN<-MUVll+X12MIN<0  ]*X12MIN*A12 

[4]  Y012MAX-MUVll+Xl2MAX>0  ]x712¥47*412 

[5]  Y012MIN+MUVI 1+7 12MIN<0 ]*I12MIN*A12 

[6]  K012MAX*-M12M4X*(.  ( ( ( I XI  2MIN  ) [ \Xl  2 MAX  ) * 2 ) + ( ( ( | Y12MIN)  T I 712 
M47  )*2)  + 712M47*2)-s-412*2 

[7]  KOI  2MIN*-M1  2MIN*  ( ( ( ( | 712.^1,7)1  | 71  2MAX ) * 2 ) + [ ( ( |712MI7)l  I 712 
jV/47)*2)  + 712WJ.V*2)t^12*2 

[8]  702  3MAX*-(MUV*-M2  3 MIN  ,M23MAX)L  1+X2  3 A/47  >0  ]x72  3 MAX*  423 

[9]  702  3MIN*-MUVL  1 +72  3MIN<0  ] x72  3 A/J7-5-  42  3 

[ 10  ] Y023MAX+MUVI 1 +72 3A/47>0  ]x72 3A/47*42  3 
[ 11  ] 702  3MIN+MUV[.  1+723  A/T7<0  ] x 72  3 MIN*  A2  3 

[12]  K023MAX*-M23MAX*{  ( ( ( |72  3A/I7)  T | 72  3 A/47 ) * 2 ) + ( ( ( | Y23MIN)  T | 72  3 

MAX)*  2 )+K23MAX*2  )*423*2 

[ 13  ] 702  3 MIN-*- M2  3 MIR*  ( ( ( ( |723¥I7)L  | 72  3A/47  ) * 2 ) + ( ( ( | Y23MIN)l  | 72  3 
MAX)*2)+X23MIN*2)iA23*2 

[ 14]  704  3 MAX*-  ( MUV-*-M^3  MIN  , A/4  3 A/47  ) [ 1 +74  3 A/47  >0  ] x74  3 MAX*  44  3 

[15]  X0^3MIN*-MUV\.  1+743 MIN  <3  ] * 74 3 MIN*  443 
[ 16  ] YOA3MAX+-MUV\_  1 +7 4 3 A/47  >0  ] x 74  3 A/47  * 44  3 

[17]  Y0^3MIN-*-MUV[.  1 +74  3A/T7<0  ]x743A/T7*443 

[18]  704  3A/47-*-A/4  3A/47x  ( ( ( ( |74  3 MIN)  l I 74  3 MAX ) * 2 ) + ( ( ( | 74 3 A/17)  T | 74  3 
A/47  )*2)  + 743M47*2)t443*2 

[19]  704  3A/I7+-A/4  3A/I7x(  ( ( ( |74  3A/I7)L  | 74  3 MAX ) * 2 ) + ( ( ( |743A/I7)L  I 743 
MAX ) * 2 ) + 74  3 A/J7*  2 ) * 44  3 * 2 

V 

V Z+FBOPT ; PAIR \I\J\ MUV ; PQ 

[1]  n PERFORMS  DYNAMIC  OPTIMIZATION 

[ 2 ] L IM- M 2 3 MIN , M2  3 MAX ,702  3 MIN ,702 3MAX 

[3]  -*•(  SYNC  ASE=  ' 2 * )/DLO[  1 ] + 5 

[4]  7012 + SMC  0 [ 1 ] O 702  3 *- SMC  0 [ 4 ] O 704  3+-SMC01  6 ] 

[7]  -+LS4 

[8]  702  3+-ITC  0 [ 2 ]<>  701  2-+5FOO[  3 ] + 702  3 

[10]  LBA : GFCLA 

[11]  'THE  FREE  MASS  PARAMETERS  FOR  THE  OPTIMIZATION ' 

[12]  ' PROCESS  THAT  IS  ABOUT  TO  BEGIN  ARE:''M23''  AND  ,,7023,, 

t 

[13]  'WHERE,1 

[14]  ( 7 p ’ ' ) , (vLIMl 1] ) , '<M23<'  , (vLIMl2l ) , ’ AND  ' , ( f LIMl 3 ] ) , » 5 

702  3 S • ,t£IA/[4] 

[15]  i0 

[IS]  ' PLEASE  ENTER  A PAIR  OF  NUMBERS  INDICATING' 

[17]  'HOW  MANY  EQUALLY  SPACED  VALUES  IN  THE  ABOVE ' 

[18]  'INTERVALS  YOU  WANT  TO  CONSIDER,  INITIALLY,  FOR' 

[19]  ' ' ' M23' ' AND  ' ' X023' ' , IN  THAT  ORDER:' 

[20]  -+(□£<?[  1 ] + 3 )x  i 3 = p 1 .PAIR+U 

[21]  UARBOUT  70  -*•□£<?[  l]-2 

[23]  AA/23«-( -/LIMl  2 1 ] ) + P4I7  [ 1 ] -1 

[24]  A7023-+C  - / LrA/[  4 3 ] ) * PAIRl  2 ] - 1 

[25]  I-J+-0F+-0 
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[26] 

[27] 

[28] 

[29] 

[30] 
[31  ] 
[35] 
[38] 

; 39; 

[40] 

[43] 

[45] 

[46] 

[47] 
[50] 

[53] 

[54] 

[55] 
[56  ] 

[57] 

[58] 

[59] 

[ 62  ] 

[ 63  ] 

[64] 

[65] 

[ 66  ] 

[67] 

[68] 

[70] 

[71] 

[72] 
[76] 

[ 77  ] 

[ 81  ] 
[82] 

[83] 

[84] 
[86] 

[ 87  ] 
[91] 
[92  ] 

[94] 

[95] 
[96  ] 
[97  ] 
[98  ] 
[99] 

[ 100  ] 
V 


75  1 : MlJ2  3-*-L  IMZ  1 ] + IxaW23 
LB2  :X022*-LIM\_ 3 ]+Jxa£523 
-753x; 0=1  ' FINDMP_'  , SYNC  ASB 
-753 0 = i' CA LCOF ' ,5755455 
— 753  x t 55S  55 

OF+-OFO  M23+-M0230  £523-£5230  £523-£523 
£512-£5120  £54  3 -£54  30  £543-£543 


753  :— 752  x \ 54 75 [ 2 ] >7—7+1 
— 751  xiPAIRl  1 ] >7— 7+1 +7-0 
□455557  7 0 M5K-D57  10  □455557 


7 


555740  -755 xi 05*0 

'SORRY,  NO  SOLUTION  COULD  BE  FOUND.' 

-+LB& 

LBS  :M02  3-A7230  £02 3-£02 3^  £523-5023 
£512-55120  £54 3— £54  30  £543-5543 
*£5.  (75  DESIGN  PVINTS  CONSIVERED=' ,yx /PAIR+1 
'RANGES:  ' , (t7IA#[1])  , ' <M23<'  , (t77M[2])  , ’ 47/5  ' , (t77A/[  3 ] ) , 
’ <£023<’  ,t>LIMZ  4 ] 

’555:7  SOLUTION :' 

( 1 5 p ’ '),  'SCORED'  ,t  55 

( 1 5 p ’ '),' (IDEAL  SOLUTION  WOULD  SCORE  100)’ 

( 1 5 p ’ ' ) ,' M23  = ' ,VM023 

( 1 5 p ’ ’ ) , ’£523  = ’ ,t£5230  i OO  *0 

LIM+M02  3 + ( ( - LM2  3 ) , A,V  2 3 )*M023*M23MIN ,M23MAX 

LIM+-LIM  , £02  3 + ( ( -A5523  ) , LX  023  ) * X023XX 023MIN ,X 023M AX 

’575455  MAKE  YOUR  SELECTION :’ 

' 1)  PERFORM  ANOTHER  SEARCH  AROUND  CURRENT  OPTIMUM ' 

' 2)  555V  PLOT  FOR  DYNAMIC  PROPERTIES' 

-(  ( 0 = pp,V57)  AA/.V57e  ’ 12’  )+ ( l + ’12’i  1 + 7757-0  ) 4>754 ,757 
□455557  70  -D75[ 1 ]-2 
757 : 5455- ' 7 ’ 

STADYNOUT 

GFCLG/>  DRWFB 20  '755  OPTIMUM  MASS  5454.7757555  455: ’ 0 lO 
(tAM12[1]),»  <,M12<  ' , T A/V 1 2 [ 2 ] 

’ £512  = ’ ,t£5120  ' 751  2 = ' , t>751  20  t0O  MUV+\/LINKSl  2 ; 1 2]*7  75 
£5  [ 2 ; 3 4 ] 

' 77’  , (55-M574>’  2 3 ’ ) , ’ = ' ,VMO 2 3 
' £’  ,55,  ’ = ’ , T ( £2  3,423  - £2  3 ) [ MUV + 1 ] 

’7'  ,55,’=’  , T (72  3 , - Y23~J[MUV+1  ] 

'£’ ,55, ’ =’  , t£2 30  i 0 

( T A.V4  3 [ 1 ] ) , • <M43<  ' , (tAV43[2]  ) , ( (pV57)x4=pAV43  ) + .V57-' 

55  * , (T  (“2+AM43  )[  1]  ) , ’ SA74  3 < ’ , T>  ( ~2  + AM4  3 ) [ 2 ] 

' £543=’ ,t£5430  ’ 754 3 = ’ , t754 30  ’ £54 3 = ’ , t£54 30  iO 
'NOTE:  £12  CAN  ASSUME' 

' ANY  ARBITRARY  VALUE. '0  iO' 

756: 'PLEASE,  CHOOSE  ONE  ALTERNATIVE: ' 

' 1)  RETROCEDE  TO  MODIFY  MASS  PARAMETERS' ' LIMITS' 

' 2)  RETROCEDE  TO  REDO  SYNTHESIS' 

' 3)  RETROCEDE  TO  MODIFY  OBJECTIVE  FUNCTION' 

’ 4 ) END ' 

-(  ( 0 = ppZ ) a a/ (Z-0) e ’ 1234 ' ) + 0 
□455557  70  -D75[ 1 ]-2 
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V Z+FBSFCOF.F\U\ V\CO\IN 

[1]  n CALC.  COEFFS.  FOR  EQUALITY  CONSTRAINTS 

[2]  n IN  S. FORCE  SYNTHESIS. 

[3]  IN_-*-~2  + 1 AIN_ 

[4]  C0+{  ,SFLl  ;IN1-PRCPTSL)%(U  ,[l]7)  , ( - V+-V2  DDL  IN  ] , [ 1 . 5 ] 73  DDL  I 

in)  ,LllU+U2DDZlin,ll.SlU3DDlIjn 

[5]  r\NOW  NEC.COND.  ARE  GOING  TO  BE  CHECKED 

[6]  Z+0 

[7]  -*■(  ( (X01  2MIN+M2  3MIN-X02  3 MAX  )zC0ll'})*C0Lll£X0l2MAX+M23MAX- 
X023MIN) /O 

C 8 ] -*•(  ( {Y012MIN-Y023MAX)<iC0Z  3 ] )*COl  3 1<Y01  2MAX- Y02  3 MIN)  /O 

[9]  -*■(  ( (X0^3MIN+X023MIN)<C0Z2l)*C0l2~]<X0>A3MAX+X023MAX) /O 

[10]  -*-(  ( ( Y0^3MIN+Y02  3MIN)<C0ZH]  )*COL  4 ] <Y54  3MAX+Y02  3 MAX ) /O 

[11]  KSFCO+SFCO 0 SFCO+CO 

[13]  KCURVE+ CURVE 

[14]  SFI«--(  2,ff+l)p<  (5,[l]7)  , (-7^7255, [1.  5] 7355)  , [ 1 ] U+U2DD  , [ 1 . 
5 ]5355) + . *SFCO 

[15]  CURVE*-SFT*-SFL  + SFI_ 

[16]  KPRCPTS+-PRCPTSL 

[17]  KI A IN+-IAIN 

[18]  KFI1T5T+FTRST 

[19]  KR-R 

[20]  FIRST+Z+1 

V 


V Z-FBSMCOEFi  U 1 ;71  ;52  ;72  ;U3  ;73  ;54  ;74  ; IN_ ; CO ; D3  1 2 ; 54 1 2 ; £2  2 3 ; 
D323  ; 54  4 3 ;5343  ;5243 

[1]  n CALC.  COEFFS.  FOR  EQUALITY  CONSTRAINTS 

[2]  n IN  S. MOMENT  SYNTHESIS 

[3  ] PRC  PTS  L-*~PHI1 2 ZlN-*-l+l  . 5+PRCPTSLl  1;]*360*F],[.5  1 PRCPTSLl  2 : 

] 

[4]  PRCPTSLD<-(PRCPTSLl  l;]xi80fOl),[.5  hPRCPTSLl  2 ; ] 

[5]  Ul+Ul-RPOINTlllQ  V1+V1-RP0INTZ2] 

[7]  U2+-U2  -RPOINTZ  1 ]0  V2+-V2 -RPOINTZ  2 ] 

[9]  U3+U3-RPOINTZHQ  V3+V3 -RPOINTZ  2 ] 

[11]  U^U*- RPOINTZ  HO  J£4-**74  -RPOINTZ  2 ] 

[13]  D312*-(U1*U2DD)  + V1xV2DD 

[14]  D^12*-(U2*V2DD)  -V2*U2DD 

[15]  52  2 3-«-(  (U2xV3DD)-V2*U3DD)  + [1l3xV2DD)-V3xU2DD 

[16]  D 3 2 3«-(  (U2xU3DD)  + V2xV3DD)-(U3xU2DD)  + V3xV2DD 

[17]  5443+-(  53x7355)  - 73x5355 

[18]  5343«-(  U^xU3DD)  + VkxV3DD 

[19]  5243-<-(  54x7355) -74x5355 

[20]  C0*-(SMLZIN1-PRCPTSLZ  2 ; ] )@53  12[IV]  ,54  12 [IT/]  ,52  2 3[_T£]  ,5323 
[15]  ,5443  [J5]  ,5343[I£]  , [1  . 5]5243[I5] 

[21]  *NOW  NEC .CONDITIONS  ARE  GOING  TO  BE  CHECKED 

[22]  Z+-0 

[ 23  ] -*-(  (r012MI5<55[  1 ] )*55[  1 ]ZY01  2MAX ) /O 

[24]  -*•(  ( ( (X012MIN+M23MIN  + K023MIN)-2xX023MAX)<C0Z2l  )*55[2]<J51 

2MAX+M2 3MAX+K02 3 MAX- 2xX02 3 MIN ) / 0 

-*■(  ( (*02  3AfI5-£52  3M/};OsC,0[3]  )*55[  3 ]<Z52  3MAX-K02  3MIN) /O 
-*(  (752  3AfJ/i/<55[4  ] )*55[  4 ]ZY02  3MAX ) / 0 


[25] 

[26] 
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[27]  +( ( (K023MIN+K0H3MIN)ZC0lSl)*C0l5lZK023MAX+K0H3MAX) /0 

[ 28  ] ->(  \Y0^3MIN<,C0l  6 ] )*C0  [ 6 ] <7£4  3 MAX  ) /0 

[29]  + ( ( UOH3MIN-KOH3MAX)ZCOm)*COl7lzXO*3MAX-KO±3MIN) /O 

[30]  KSMCO*-SMCO 0 SMC 0+ CO 
[32]  KCURVE*-CURVE 

[ 33  ] SMI*-{  SMC  0[ljx£>312  ) + ( 5AfC0[  2 ] x Z74 1 2 ) + ( SMCOl  3 ] x D2  2 3 ) + 5WC’(7[  4 ] 
*D323  — 

[ 34  ] CURVE*-SMT*-SML+SMI*--SMI+(  SMCOl  5]x£)4  43)  + ( SMCOl  6 ]x  D3  4 3 ) + SMC 

0[7]x£>243  

[35]  KPRCPTS+PRCPTSLD-y  KFIRST+FIRST FIRST*-Z*- 1 
[38]  KRMS*rRMS/>  KPTP+PTPC  KMIN*-MIN'>  KM  AX*- MAX 
V 


V Z*-F.INDMP1 

[1]  A SOLVES  EQUALITY  CONSTRAINTS  OBTAINED 

[2]  A FP0A7  S. MOMENT  SYNTHESIS 

[3]  Z*-0 

[4]  -( (K023MIN<K023 )*K02  3M  AX>  K02  3*-X02  3 -SMC Ol  3 ] )/0 

[5]  -(  (ZOl2A/JiV<JOl2  )*ZOl2M4J>JO12«-SMC0[2]+(  2x^02  3 )-M023+K023 
)/0 

[6]  -►(  (,K0^3MIN^K0^3  )*K0^3MAX>KQi+3*-SMC0l5l-K023  )/0 

[ 7 ] -*•(  kX0^3MIN<X0K3  )*X04  3MA  X>  XOK  3*-SMC0l  71+K0H3  ) /O 

[8]  Z+BASICMP 

7 


V Z*~FIN  DMP2 

[1]  A 5/7£K5’5  EQUALITY  CONSTRAINTS  OBTAINED 

[2]  A FffOM  I .TOR+S .FOR  SYNTHESIS . 

[3]  Z+0 

[4]  -*•(  (,X012MIN<X012)*X012MAX>X012*-SFC0111  + X023-M023  ) /O 

[5]  -*■(  (X0^3MIN<X0H3  )*X0H  3 MA  X>  JQ4  3*-SFC0  [ 2 ] -X02  3 ) /O 

[6]  ->-(  ( Y 0l\3MIN<:Y  0^3  )*YOH  3MAX>Y0H  3*-SFC0  [ 4 ] - Y02  3 )/0 

[ 7 ] -*(  (K023MIN<K023  )*K0  2 3MAX>K02  3*-X02  3 -ITCOj  1 ] ) /O 

[8]  -*(  (A:^4  3/^I//<A:04  3')wjf04  3Af^Z'>A:C>4  3-<-IfC,(7['3T-A:(92  3 )/0 

[9]  Z*-B  ASICMP  

V 


V Z*-ITSMSYN ; FLAG;  MUV ;MINPTS;  KEEPSM ; KC URV E ; PRC PT S ; KEEP ; DIST 
; .Z7VZ? ; PRCPTSL  ; T;C  ASE;  FIRST ; PP;  RMS;  KRMS',  PTP',  KPTP;  MAX ; KM  AX  • 
tf.ZW ; KM  IN 

[1]  PERFORMS  I .TOR  OR  S. MOMENT  SYNTHESIS 

[2]  TRC*-'T' 

[3]  P.4SP-*-14  ~1  2[  iSIWC’/JSP]  + 1 SHAKING  MOMENTINPUT  TORQUE ' 

[4]  FIRST*-Z*-FLAG*-  0 

[5]  MIN  PTS*-(  7 6 12  12[iM  m,  3 3 6 6[*tf  T7]  ) [ iS7tf<74SP] 

[6]  ~ 

[7]  1 TO  PERFORM  '.CASE,'  SYNTHESIS' 

[8]  'YOU  MUST  SELECT  AT  LEAST  ' , (tMINPTS) , ' PRECISION  POINTS 

f 
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[9]  ' I F MORE  TRAN  1 , (yMINPTS) , 1 ARE  SELECTED , THEY  WILL • 

[10]  ’ Bff  APPROXIMATED  IN  THE  LEAST  SQUARES  SENSE.1 

[11]  \0O  » PRESS  11  RETURN11  TO  CONTINUE  . * 0 MUV+- ffl 

[14]  0 

[15]  DP*-  2 2[iSy/VC\4S£]  + 'SMIT1 

[16]  S.DP , 1 C 0-*-  KC  URV  E 1 

[17]  CURVE*-i(2  21  iSYNCASE]*  1 SMIT1  ) , »T' 

[18]  PRCPTSLD*-PRCPTS*-2  Opl 

[19]  LB2iGFCLG 

[20]  444  PLOT  PH  1 12 , l . 5lC URV E 

[21]  PRCPTS*-GFLTP  PRCPTSLDC>  GFACT  10  DRWSTAR  PRCPTS 

[24]  5 60  GFTXT  1 FOR  THIS  CURVE : RMS=' ,(y9  3iRMJT7'  PTP=',v9 

~3vPTP 

[25]  +(□£<?[  l]  + 3)xlo  = pKPtfP7P 

[26]  5 40  ffFTTr  ' POP  B4S27  PPPI7P:  RMS=',( T>  9 -3t  KRMS),1  PTP=',V 

9 3 T KPTP 

[27]  5 2 5 GF2Wl6p  » ’ ) , 'MIN-  1 , (t9  '3  iKMIN),1  MAX-1  ,t9  “3t  KM  A 
X 

[28]  GFMOV  0 7790  □4PBOP7’  13 

[30]  -IBS*\FLAG=0 

[31]  FLAG*- OO  ’ P0B5  TBIB  CURVE  SATISFY  YOU?1 

[33]  i' -*■□££+'  , 1 + ( ~ 1 + » Y47  * i l + CZDct)’  341* 

[34]  UARBOUT  70  -ULCl l]-20  +~Z^1 

[37]  BB3 : 1 TO  SELECT  A POINT  PRESS  ANY  KEY  OTHER1 

[38]  ' THAN  E,T,S,I  OR  Q,  THEN  PRESS  RETURN . 1 O i0 

[40]  'USE:  ''E11  TO  ELIMINATE  A SELECTED  POINT 1 

[41]  ' ''T11  TO  THROW  AWAY  ALL  CURRENTLY 1 

[42]  * SELECTED  POINTS 1 

[43]  ' ''S11  TO  SYNTHESIZE  '.CASE 

[44]  ’ "I"  TO  INTERCHANGE  THIS  CURVE  WITH1 

L45]  ? THE  PREVIOUS  ONE  ',(19x0 = p KC URV E)  i 1 ( C AN  1 1 T BE 

USED  NOW)1 

[46]  » ''Q11  TO  QUIT  SYNTHESIS  AND  RETROCEDE 1 

[47]  ' TO  MODIFY  I N E Q . C ON STR AINTS . 1 0 lO 

[49]  'NOTE:  2 BEEPS  MEANS  NECESSARY  CONDITIONS' 

[50]  1 CANNOT  BE  SATISFIED . 1 

[51]  LB  4 3+ ("3+3x102  106  1 17  114  1 1 6 i ~1  + KEEP+-GFGCUR ) <$> 1 LB5 

BB9BB0000BB5  LB1 1 

[52]  BB5  :-*■(  ULCl  1 ] + 3 )x  x 0*  ltp  PRCPTS 

[53]  UARBOUT  7 0 -*LB^ 

[ 55  ] IND*-l+lDIST<-+-/  (PRCPTS -§(bp  PRCPTS)  p-l  \KEEP)  * 2 
[56]  -*■(□££[  1 ] + 3 )x  x HZDISTlINDl 

[ 57  ] UARBOUT  70  -*-£B4 

[59]  ( .PRCPTSliIND] )GFCIRC  4 16 

[60]  PRCPT^RI  WINDbPRCPTS O +BB4 

[ 52  ] LB1  :->(DBP[  1 ] + 3 )x  l ( 4 SWzKEEPl  1 ] ) A ( KEEPl  1 ]£ 97  8 ) *~KEEPl  1 ] ePRC 

PTS [ 1 ; ] 

[ 63  ] UARBOUT  70  -*-£B4 

[55]  DRWSTAR  2+ KEEP 

[66]  PRC  PTS*-PRCPTS  , 2*KEEP/>  -*-BB4 

[68]  BB6  :-*■(  ULCL  1 ] + 3 )x  \MINPTS^  1 + p PRCPTS 

[69]  UARBOUT  70  ->-£B4 

[71]  GFACT  44  4 ^ PRC PT SL+GFPTL  PRCPTS A GFACT  1 
[74]  -*LB3  x x 1 = i ( 4 2 p ' FBSCWA  ST  1 ) [ zM  TY ; ] , DP  , 1 COEF 1 


182 


[75]  UARBOUT  7 0 T*UDL  10  □ ARBOUT  7 0 + 554 

[ 79  ] 558  :F540^10  + 552 

[31]  LBO  :PRCPTS_LD*2  OplO  + 552 

[83]  559  : + ( []5C[  1 ] + 3 )"i  0*pKC  URV  E 

[84]  UARBOUT  70  + 554 

[86]  MUV*CURVEO  CURVE* KCURVEO  KCURVE+MUV 
[89]  MUV-FIRST 0 FIRST+KFIRSTO  KFIRST*MUV 
[92]  MUV*RMS/S>  RMS+KRMS*>  KRMS*MUV 
[95]  MUV*PTP'>  PTP*KPTP 0 KPTP+MUV 
[98]  MUV-MIN^  MIN*KMIN?>  KMIN+MUV 
[101]  MUV*MAXO  MAX*KM  AX~>  KMAX-*-MUV 

[104]  MUV*tDP,'C0'O  SlDP,'C0'  ,'*K'  ,DP,'CO'^  i' K' , D P , ' C 0*MUV ' 

[107]  GFACT  444 

[108]  MUV+GFPTL  PRCPTSO  PRCPTSLD*KPRCPTS 0 KPRCPTS*MUV 0 FLAG*FI 
RSTy  + 552 

V 


V W-*- KEYBOARD  IN  PANG IN  PA  ; LI ; LV ; RI  i RV ; MUV 

[1]  a PSFO  FOP  IPPOF  OF  LOADS  BY  TYPING  IN 

[2]  A F40P  INDIVIDUAL  045PF. 

[3]  IN P A* IN P AUG* 180*01 

[4]  LBS-.GFCLA 

[5]  V+i57+0 

[6]  LI* A 

[7]  LBAi'YOU  WANT  TO  CONSIDER  THE  INP  ANGLE ' 

[8]  ’FPOM  ’ ,(T5tO,'  TO  WHAT  VALUE?  {IN  DEGREE )' 

[9]  + ( 3+[]5P[l]  )xt  (0  = ppff7)A(  a / 50<P7 ) a a / ( PlM] ) < 3 60 

[10]  CMPSOPr  7 

[11]  +[]5P[ l]-2 

[12]  Py+IPP,4[PI+lU  IIPP4-PK] 

[13]  ’00F5  FPF  PARAMETER  TO  BE  INPUT ' 

[14]  ’P47F  CONST.  VALUE  FOR  THIS  INTERVAL ' 

[15]  ’ OF  FPF  IPP.  ANGLE? ' 

[16]  4 ’+’  , 5+  ("5  + 5x»  YP’  ; l + a)<})»LSl  582  G50+1' 

[17]  Q<4PS0PF  7 

[18]  +Q50[ l]-2 

[19]  LBli'WHAT  IS  THAT  VALUE?' 

[20]  t/+P,  ( 1 + RI-LI)  pl  + G 

[21]  LI+-RI+1 

[22]  +553 

[23]  552 PLEASE  ENTER  VALUE  OF  PARAMETER ' 

[24]  ’FOP  IPP  4P05F  = * ,TlPP4[5I] 

[25]  W*W, 1 + G 

[26]  +(Q50[  l]-2  )*\RIZLI*LI+1 

[27]  553  : LV*INPA [ LI-MUV*( pIN  PA) < LI] 

[28]  *LBA*\MUV=0 

[29]  *LBS*\ A/f/  = y[l] 

[30]  TRC*'T' 

[31]  GFCLG 

[32]  7.89  PLOT  IN PANG ,[ . 1]P 

[33]  □ .IPS  OOF  13 

[34]  'SELECT:' 
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[35]  * 1)  TO  CONTINUE ' 

[36]  ' 2)  TO  RE-ENTER  DATA ' 

[37]  ■+(  ( 0 = ppty£/lO  AA/MZ/Kel  2 ) + i 5+ ("5  + 5x1  2 i MUV*-U  ) <t> T ULC+  3 LB6  ' 

[38]  UARBOUT  7 

[39]  -+ULCI  l]-2 

[40]  GFDREG  7.89 

[41 ] LB  5 : GFACT  1 

[42]  GFCLA 

[43]  GFMOV  0 779 

[44]  UARBOUT  13 
7 


V Z*-LO  ADS ; SH  ; I •,  MUV ; LIN  ; II GFLAG  ;WHILE 

[1]  n HANDLES  INPUT  OF  EXTERNAL  LOADS 

[2]  GFCLA 

[3]  i'DRW'A  4 2 p ' FBSCWAST ' ) [ iM_TY ; ] , * 2 * 

[4]  Stf-*-l  t p LINKS 

[5]  ' INDICATE  ALL  LINKS  UNDER  EXTERNAL  LOADS :' 

[6]  I<-1 

[7]  ' '.(•Pi),’)  OlffF  1 , LINKSlI ; 3 4] 

[8]  ■+■(  DLC[  1 ] -1  )x  x SH^I*-I+ 1 

[9]  ? ' .(Tl+Sff).*)  THERE'' S NO  LINK  UNDER  LOAD' 

[10]  -*■(□£<?[  l]  + 3)xx  {A/NUV=l+SH)v*/{MUV+-U)€\Sa 
111]  GdflBOW  70  +ULCI l]-2 

1 13  ] LIN  KSUL*-(.  ( \SH)efiUV)/\SH 

1 14]  -+(  (SHrn:)eMUV)iZ*-0 

[15]  2 3 p ' * ' 0 ' 1)  FOFOF  OFFI' 

[17]  » 2)  MOMENT  ONLY' 

[18]  ' 3)  FORCE  AND  MOMENT' 0 23p'*' 

[20]  'B/4SFD  ON  OPTIONS  ABOVE,  INDICATE' 

[21]  ' KIND  OF  EXTERNAL  LOAD ' 

[22]  LOKND*-\ 00  Z«-J«-l 

[24]  LB1:'F0R  LINK  ' , ( LIN  KSL  LIN  KSULl1 ] ; 3 4]),':’ 

[25]  ■+([]£C’[l]  + 3)xl(0  = p pMUV  ) a a / (MUV+-Q  ) g x 3 

[26]  UARBOUT  70  +ULCI l]-20  LOKND-LOKND ,MUV 
[29]  -*-LBl  x x ( pLIN  KSUL  ) S J-+J+1Q  Ji-1 

[31]  LB1  :LIN*-LINKSULlIl O -+BB2  x t 2 = LOKNDl  I ] 

[33]  'FOB  BIBF  ',  (BIBBS [£IAH-£I/VFSO£[I ] ;3  4]),',  Z"S  IFF  POIBI 
OF  APPLICATION ' 

[34]  ' OF  FIT  FORCE  CONSTANT?' 

[35]  i’-+*  ,5+("5  + 5x'y^*xl  + n)^»ifl3  £B4  □iO+1* 

[36]  UARBOUT  70  -*-GBC[  1 ] -20F83  : i 0 

[39]  'ENTER  THE  X-Y  COORDS  OF  THAT  POINT:' 

[40]  -»-(nOO[l]  + 3 )xi  3 = pl  ,MUV+-U§  UARBOUT  70  1 ] -2 

[43]  sl'E'  ,LINKSLLIN‘,1  2],'«-2  lp  ',  vTBLl  LIN  ; 2 3 ] + IBF[  BIB;  1 ]BOIr/ 
FO  M UV 

[44]  ->-BB50BB4  :II+-lO  WHILE*-  ( 2 , (B  + 1 ) ) p 0 0 iO 

[48]  ’70V  00  JOB  <7.401  TO  ENTER  THE' 

[49]  ( 'XY' [II] ) , ' -COORDINATES  OF  THE  POINTS  OF  APPLICATION ' 

[50]  'OF  FBI  FORCE  ON  LIN  K'  , LIN  KSl  LIN 3 4],»?» 

[51]  ' 1)  CALCULATOR ' 

[52]  ' 2)  KEYBOARD ' 
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[53]  -(  3+[]L<7[l]  )*i  (0  = ppMUV)**/(MUV<-U)  el  2 

[54]  □ ARBOUT  70  -+ULCZ  l]-2 

[56]  WHILELII; >j>( 10+ ( 10*MUV-1 )§' C ALCULATORKEYBO ARD  ' ),’  Pffl 
' ,£IVKS[1 ;1  2] 

[57]  -*-(LB4  + 2)*i2  = JI^JJ+l 

[58]  WHILE+($((N+l) ,2)pTBLlLIN;2  3D  + TBLI  LIN  ;1]R0TVEC  WHILE 

[59]  ,LI/V£S[LIP;  1 2 ] , » £ ’ 

[60]  LS5  :II+10  WHILE*-(  2,  (N+l  ) )p0 

[62]  'indicate  coppp  system  for  the ' 

[63]  FORCES  {LINK' ,LINKSlLIN\3  4],*)  YOU"  RE  GOING  TO  EN 
TER : ' 

[64]  » 1)  GLOBAL' 

[65]  * 2)  LOCAL' 

[66]  -*■(□££[  l]  + 3)xi  {0  = ppGFLAG)**/{GFLAG-U)el  2 

[67]  CMPflOJ/r  70  +QLCZ l]-2 

[69]  LB3  :'  HOW  DO  YOU  WANT  TO  ENTER  THE' 

[70]  ((2  2p 'UVXY' )lGFLAG;IIl ) , » -COMPONENT  OF  THE  EXT  FORCE  ( L 
INK' ,LINKSlLIN;3  4],')?* 

[71]  » 1)  CALCULATOR' 

[72]  1 2)  KEYBOARD' 

[73]  -(3+QLC[l])xi  (0  = ppMUV)**/ (MUV+n)  el  2 

[74]  □ ARBOUT  70  -*■□££[  l]-2 

[76]  WHILELII;  ]«-j>  ( 10+  ( IQxMUV-l  )<t>' CALCULATORKEYBOARD  '),'  PHI 
', LINKSl 1;1  2] 

[77]  -LB3*\ 2=11^11+10  +(D£C[l]+2  )xx  (GFL4G=1  ) 

[79]  WHILE-i (180*oi)x4»pffi* , LI0KS[ LIP; 1 2 ] ) R0TVEC2  TBLl LIN ; 1 ] 
ROTVEC  WHILE 

[80]  i'F' ,£IPKS[LIfl;l  2l,'+WHILE' 

[81]  LS2  i-*-LB&x\  l=-L0KNDlI1 

[82]  '#<?(/  DO  YOU  WANT  TO  ENTER  THE' 

[83]  'EXT  MOMENT  ON  LINK ' , LINKSL LIN ; 3 4],’?i 

[84]  ’ 1)  CALCULATOR ' 

[85]  » 2)  KEYBOARD' 

[36]  + (3+[]LC,[l])xi  (0  = ppMUV)*K/ (MUV+U)  €l  2 

[87]  UARBOUT  70  -□££[  1 ] -2 

[ 89  ] i * Z"  , LINKSiLIN ; 1 2],,-*-',Ti(10+(  lOxtfpy-i  )4>'CM£Cm4T0P££YB 
0y?F0  » ),’  PHI' , LINKSl  1 ; 1 2] 

[90]  £56 ;i( llx( pLINKSUL) >1+1+1 H 'DRW' , (4  2 p 1 FBSC WA ST ' ) [ iM_TY ; ] 
, » 20*  LB1  ’ 

V 


V MPINjMUV; VAL 

[1]  n HANDLES  INPUT  OF  MASS  PARAMETERS 

[2]  GFCLGO  t'  DRW'  , (4  2 p ' FBSCWAST ' ) [ ±MJTY;  ] , ' 2 • 

[4]  I«-lO  • VP/4T  MASS  PARAMETER' 

[6]  'DO  YOU  WANT  TO  CHANGE?' 

[7]  MUV*f(t' M' ,MUV) ,{i'X' , MUV ) , (l'Y' ,MUV ) , * ' K ' , MUV+LIN KSlI i 1 

2] 

[8]  ' M'.MUV,'  X' ,MUV ,'  Y ' ,MUV , ' K' , (MUV+LIN KSlI ;3  4]),'=’, 
MUV 

♦(□CC[l]-2)xi ( 1+p LINKS)*I*I+1 


[9] 
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[10]  * ( FOR  EX:  ENTER  ' ' M'  , ( LINKSl 2 ; 3 4]),,T'  TO  ADRESS  THE  MA 
SS  OF  LINK  LINKSL 2 ; 3 4] 

[11]  'OR  ANSWER  " NONE " TO  CONTINUE )' 

[12]  £B2  : + (l  = A/’/V0/y£’=4tAfi/tML)*O 

[13]  -*■(□££[  l]  + 2 )xl~A/(‘2tMneTi9 

[14]  -t-tBlx!  (4  = p ’ 1*  ,MUV)*{  ( jt  2AMUV)  e£ , LINKSl  ;3  4 ] , ( 1 + p C INKS)  p ’ 

' )a( 1+VtftO e'MXYK' 

[15]  UARBOUT  70  +DLC[l]-4 

[17]  LSI: 'WHAT  IS  THE  NEW  VALUE  OF  '.MUV,'?' 

[18]  -+{ULCl  1 ] + 3 )x  i 0 = ppIMZ>[] 

[19]  UARBOUT  70  +ULCI l]-2 

[21]  i ( ItAfZW)  ,LI/VX5[  ( a/LI^X5[  ;3  4 ]=  ( ( 1 + pLINKS)  , 2 ) p 2+MUV)\l;l 

2] , ,^K4L* 

[22]  'WHAT  M.P.  DO  YOU  WANT  TO  CHANGE  NEXT?'0  +LB2 
V 


V MPLIMMOD \MUV \V AL 

[1]  n ALLOWS  MODIFICATION  OF  MASS  PARAMETERS 

[2]  GFCLGD  £ ' DRW  * , ( 4 2 p ’ FB SOW  AST ' ) [ t M_TY; ] , 1 2 ’ 

[4]  I-*-l 

[5]  LIMIT  DO  YOU  WANT  TO  MODIFY?'')’  lO 

[7]  MUV+l{£'M'  ,MUV)  ,{t'X'  ,MUV ) ,{t'Y'  ,Aftf7)  ,1’X'  , MUV+LINKSlI ; 1 

2 ] , ’ MIN ' 

[8]  ' M',MUV,'  X'.MUV,'  Y ' ,M  UV  , ' K'  ,{M  UV*- LINKS  [ I ; 3 4 1,'MI 
N')  ,'='  ,MUV 

[9]  +([]£C[l]-2)xi  il*pLINKS)H-I+lO  lOO  1+ 1 

[12]  MUV+y(£'M'  ,MUV ) , {£' X'  .AflW)  , (±' Y'  ,M)  , i’K’  tMUV-*-LINKSlI  i 1 

2] , » AMY* 

[13]  * M'.MUV,'  X'.MUV,'  Y'.MUV,'  K'  , {MUV+-LIN  KSlI  \ 3 4 ] , ’ AM 
X'  ) ,'='  ,MUV 

[14]  — ( □£C’[  1 ]-2  )x!  ( lipLINKS)>I+I+l 

[15]  lOO  '{FOR  EX:  ENTER  ''M', {LINKSl 2;3  4 ] ) , * MIN ' ' * , * TO  MOD 
IFY  IT' 

[17]  'OR  ANSWER  ''NONE''  TO  CONTINUE )' 

[18]  LB2  :-*(  1 = */' NONE'  =4  + MIAMI])  + 0 

[19]  -*•(□££[  1 ] + 2 ) * i~(  A/ ( 2+ l4>AWy)€»  1234567  89*  )A7=p’l'  ,MUV 

[20]  +LBl*\U.±2+lpMUV)ei,LINKSh3  4 ] , ( 1 + p LINKS ) p ' ' )a((1  \MUV) 

e ’ MXYK ' ) a ( a / ’ MIN ' =~3  + MUV) v * / ' MAX' = 3 A MUV 

[21]  UARBOUT  70  -*ULCl l]-4 

[23]  LB1 : ' WHAT  IS  THE  NEW  VALUE  OF  ' , M UV , ' ? ' 

[24]  -*-(  ULCl  1 ] + 3 ) x t 0 =p  pV  AL+-U 

[25]  UARBOUT  70  -+ULCI  l]-2 

[ 27 ] i( 1 +MUV) .LINKSl (a /LINKSl ;3  4]= ( ( 1+p LINKS) , 2 ) p2+ IfyMUV ) i 1 ; 
1 2], ("3  + MUV)  ,'*-VAL' 

[28]  'WHAT  LIMIT  DO  YOU  WANT  TO  ALTER  NEXT? * 0 +LB2 
7 
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V REG  PLOT  AO  ; KEEP;  EQ',M\MUV 

[1]  n PLOTS  SCALAR  FUNCTIONS  OF  INP. ANGLE 

[2]  REG  GFREG  450  195  978  585 

[3]  REG  GFSCL((*M)x( 0, (-M) , 0 , ( 8 x0 -MIN ) +MIN ) * EQ+MIN=MAX ) + ( 0 
,MIN+l/A0l2il) ,360 , MAX+X / AO [ 2 ; ] 

[4]  GFACT  REG 

[5]  KEEP+GVARGl 10  11  12  13] 

[6]  GVARGLlO  11  12  13>1  10  0 

[7]  450  195  GFRECT  978  585 

[8]  01101  GFVABSS5  2p450  0 450  779  445  769  450  779  455  7 
59 

[9]  01101  GFF/1BS65  2p450  195  1010  195  1000  200  1010  195 
1000  190 

[10]  +L1*\  (*>MUV)  v{MUV‘-lKEEPlHl  + . 5)  >7  6 5 

[11]  ^ ( 2 3 p ( 1 95  )A(AfZ/K<585  ) )/»  {A30 ,MUV) GFLINE  978, MUV' 

[12]  (450 ,MUV)GFLINE  AA6 ,MUV 

[13]  ( 450  , (Afi/7-4  ) ,2  )GFTXT  '_0' 

[14]  LI : 4 50  531  11  GFTXT  9 ~3 T ( ( EQ* 1 = *MAX ) v~F Q ) + MAX+ EQ*MAX* *MA 

X 

[15]  450  191  12  GFTXT  9 3v ( (EQ*~ l = *MIN)v~EQ) * MIN- EQ*MINx*MIN 

[16]  1000  170  GFTXT  'IA' 

[17]  UARBOUT  B4 

[18]  5 82  1 75  1 GFTXT  ' 90' S 7 19  175  2 GFTXT  '180* 

[20]  851  175  2 GFTXT  ,270,O  983  175  2 GFTXT  '360’ 

[22]  GVARGLlO  11  12  13 1-KEBP 

[23]  (0, (~l+"l  + p/40) pi )GFV  ABSiAOL 1 ; ]xl 8 0*01 ) , [.l]40[2;] 

[24]  PTP+MAX-MIN 

[25]  RMS*-(  (+/40[2;]*2)*“l+p4tf)*.  5 

V 


V REG  PLOT3  AUV ; COL ; MAG ; KEEP; A\ B ; MIN 

[1]  r PLOTS  VECTORIAL  FUNCTIONS  OF  INP. ANGLE 

[2]  REG  GFREG  800  220  980  400 

[3]  REG  GFSCL  0 0 ,R  ,R<-[  / MAG+( +■/ AUV12  3;]*2)*.5 

[4]  GFACT  REGS  KEEP+GVARG [ 1 0 11  12  13] 

[6]  GVARGLlO  11  12  13>1  10  0 

[7]  01101  GFVABSSS  2p615  220  1000  220  990  225  1000  220  9 
90  21  5 

[8]  01101  SFK4BS65  2p800  40  800  420  795  410  800  420  805 
410 

[9]  795  410  1 GFTXT  'V'O  990  200  GFTXT  'U' 

[11]  800  220  GFCIRC  1800  GVARGLlO  11  12  131+KEEP 

[13]  (0, Cl+~l*pAUV) pl)GFVABS  AUVZ 23;] 

[14]  GVARGLlO  11  12  13>1  10  0 

[15]  B+(  (-/AUVl2;A])ARG-/AUVl3;A<--\(~lJtpAUV)  ,1]) + 0(180  + 30  ~3  0 
)*180 

[16]  10  1 GFVINCL  . 5 + ( A , ( -A*-(  2oB[  1 ] ) , lofl[  1 ] ) , [ 1 . l](  2oB[  2]  ) , lo 
B[2]  )xio 

[17]  + E,B1*\R=MIN+-L/NAG 

[18]  800  20  12  GFTXT  'MIN.  MAGNITUDE  = ’,8  3 vMIN 

[19]  -+2  +[]LC,[  1 ] 
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[20]  LB  1:800  1 5 12  GFTXT  1 MAGNITUDE-CONST ’,8  3 VR 

[21]  GVARGllO  11  12  13h-*KEEP 

[22]  (REG+l)GFREG  800  450  980  750 

[23]  (REG+DGFSCL  0 ~360,F,0 

[24]  GF ACT  REG* 10  KEEP+GVARGZ 10  11  12  13] 

[26]  GVARGllO  11  12  13>1  100 

[27]  620  450  GFRECT  930  750 

[28]  01101  GFVABSty 5 2p800  750  800  430  805  440  800  430  795 

440 

[29]  0 110  1 GFV ABS§ 5 2p980  750  1000  750  990  745  1000  750  9 
90  755 

[30]  990  730  GFTXT  'U'O  795  430  2 GFTXT  'IA ' 

[32]  ( 6 p 0 1)GFVABS<*6  2p665  750  665  753  710  750  710  753  755  75 
0 755  753 

[33]  ( 6 p 0 l)GFtMBSS?6  2p  845  750  845  753  890  750  890  753  935  75 
0 935  753 

[34]  980  760  4 GFTXT  8 ~3tFO  797  760  GFTXT  *0’ 

[36]  620  760  4 GFTXT  9 -3T-i?0  UARBOUT  B_2 

[38]  610  74  5 1 GFTXT  ' 0 » --  6 1 0 6 70  2 GFTXT  '90' 

[40]  610  595  3 GFTXT  '180'0  610  520  3 GFTXT  *270' 

[42]  510  445  3 GFTXT  '360*0  GVARGllO  11  12  1 31+-KEEP 

[44]  (0,("l  +COL-*  l + pAUV)pl)GFVABS  AUV12-.1  ,1  .ll-AUVll;Ji<-AUVll; 
] x 1 80  tOI 

[45]  (REG+2)GFREG  270  220  570  400 

[46]  ( REG+2 ) GFSCL  0 0 360 ,RO  GFACT  REG+2 

[48]  KEEP-*  GVARGllO  11  12  13] 

[49]  GVARGllO  11  12  13>1  10  0 

[50]  270  40  GFRECT  570  400 

[51]  01101  GFVABSS5  2p270  220  590  220  580  215  590  220  580 

225 

[52]  0 110  1 £F174B5$5  2p270  400  270  420  275  410  270  420  265 

410 

[53]  265  410  1 GFTXT  'F'O  580  200  GFTXT  '14' 

[55]  ( 6 p 0 DGFVABSS6  2p270  265  257  265  270  310  267  310  270  35 
5 267  355 

[56]  ( 6 p 0 1 ) GFV ABS§S  2p270  175  267  175  270  130  267  130  270  85 

267  85 

[57]  260  397  3 GFTXT  8_~3t RO  260  217  1 GFTXT  'O' 

[59]  260  37  9 GFTXT  9 3 y-RO  UARBOUT  B3 

[61]  266  20  GFTXT  'O'O  345  20  1 GFTXT  '90' 

[63]  425  20  2 GFTXT  'laO'O  500  20  2 GFTXT  '270' 

[65]  575  20  2 GFTXT  '360' 

[66]  GVARGllO  11  12  131+-KEEP 

[67]  (0,(  1+COL)  pDGFVABS  AUVll  3;] 


V 
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Cl] 

[2] 

[3] 

[4] 

[5] 
C6] 


V Z+ANG  ROTVEC  VEC\C\S\A 

n ROTVEC  - ROTATES  VECTORS  IN  ' VEC ' BY  ANGLE  IN  ' ANG' . 
pi  ANG  - POS.iCCW)  OR  NEGACW)  ANGLE  OF  ROTATION  IN  D 

EGREE . 

A VEC  - 2*N  MATRIX  OF  N VECTORS.  IF  N= 1.  EITHER  p VEC 
= 2 1 

n OR  pVEC=2  IS  OK. 

n Z - ROTATED  VECTORS  (pZ=pVEC). 

Z-*-  ( 2 2p<7  , ( -5)  , (S+-104  ) ,C’-<-20  4<-4tfGx0Tl  30  )+  . x7BC 

V 


V Z«-<4  R0TVEC2  V\D\C  \ S\T\R\I 

[ 1 ] r ROTVEC 2 - ROTATES  EACH  VECTOR  IN  'V'  BY  CORRESPONDING 


[23 

[33 

[43 


fl 

n -4 
p4>2 

fl  7 


.4//GLB  I/V  *4’  . 

ROW  VECTOR  WITH  ROTATION  ANGLES  IN  DEGREE ; 
MATRIX  WITH  VECTORS  TO  BE  ROTATED-,  p7=2,p4 


[53  r Z - ROTATED  VECTORS ; pZ  = p7  . 

[63  A+-  ( ( 0-*-p  4 ) , 1 1 ) p/4xOvl  80 

[73  T+{C , -S)  , [ 2 3 ( (S-lOA)  ,0204) 

[83  Z^  + /2,x7,[237^(Z7,1  2)p$7 
V 


[13 

[23 

[33 


E 


4 

5 


V 


V 


Z*-A  ROTVEC 2 V D C ; S T ; R\  I 
n SAME  AS  R0TVEC2  EXCEPT  ARGUMENT 
n MUST  BE  IN  RADIANS  HERE. 


A<-{  (Op  4)  ,1 

T-(C ,-S) ,[2 
Z^-H/rx  7,  [2 


1 ) p4 

( ( 5-*- 10/4  ) .£<-204  ) 
7«-(2?,l  2)pfi?7 


'A  ' 


WHICH 


V WW^ SC 0 REPLOT  X I MUV LIN E 
[1  3 444  SFffBG  600  200  1000  600 

[23  444  GFSCB  0011 

[33  GF  ACT  1 

[43  LBl-.GFCLG 

[53  600  620  GFLINE  600  200 

[63  600  200  GFLINE  1020  200 

[73  600  600  GFDALINE  1000  600 

[83  1000  600  GFDALINE  1000  200 

[93  1000  200  GFDALINE  600  600 

[103  1000  183  2 GFTXT  'WORST' 

[113  600  183  1 GFTXT  'O’ 

[ 12  3 800  183  3 GFTXT  ' +■ ACTUAL 

[133  600  595  2 GFTXT  * 1 ' 

[ 143  ( 600  4 5 5 , /V  UV+-1  + 1 + p X ) GFTXT  'SCORE:' 

[ 15  3 ( 600  435  ,MUV)GFTXT  XL  1 ; 3 
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[16]  (600  415 ,MUV)GFTXT  Xl 2 ; ] 

[17]  ( 600  395  ,MUV)GFTXT  XL  3 ; ] 

[18]  GFMOV  0 7 7 90  UARBOUT  130  INLINE*- 0 

[21]  LB2:' INPUT  A TRIAL  VALUE  FOR  ''EXP''' 

[22]  'TO  SEE  CORRESPONDING  SCORE  CURVE:' 

[23]  -(□£C[l]+4)*l ( 2=pl ,WW) a0<1+V&MD 

[24]  LINE+LINE+1 

[25]  UARBOUT  70  -+ULCI l]-3 

[27]  GFACT  444 

[28]  ( 0, 199pl)GFVABS  ABSS , [ . 5 ] ( MSS* • 0 0 5 x t 2 0 0 ) * VW 

[29]  GFACT  10  GFMOV  0 .779 

[31]  UARBOUT  13,  (LINE-LINE+5  )plO 

[32]  'IS  THIS  VALUE  OF  ''EXP''  OK?' 

[33]  -*-±5t  (~5+5*'YN'  i 1*MUV*U)  4>  ’ 0 0 00  0[]Z;C,+  4[]LG,+1  ' 

[34]  LINE* LINE* 1 

[35]  UARBOUT  70  ->DL(7[l]-3 

[37]  GFMOV  0 7790  UARBOUT  1 3 , ( LINE+L INE+ 3 ) p 1 0 
[39]  -*-(LB2  ,LB1  )[  l+4=J-‘-J+l] 

V 


V SMRP ; MUV 

[1]  i0 

[2]  MUV*-  9 ~12ltSYNCASE~)i  ' SYN THE  SI SOPTIMIZ ATI  ON ' 

[3]  • PLEASE  ENTER  THE  SHAKING  MOMENT  REFERENCE  POINT' 

[4]  'TO  BE  USED  IN  THE  ' ,MUV ,'  PROCESS :' 

[5]  LB  1 :-»-([]£C,[l]  + 3)xi3  = pl,  RPO  INT*~U 

[6]  UARBOUT  70  -»-DLC[l]-2 

[8]  -*•(  ( ' 1 * =SYJ[C  ASE)*  ( a / RPO INT-U1  ,V1  ) ) + LB2 

[9]  'SORRY,  C AN'  ' T CHOOSE  JOINT  1 TO  BE' 

[10]  » THE  REFERENCE  POINT.'O  +LB1 
[12]  LB2 : FBSMIO  FBSML 0 SMT+SML+SMI 

V 


V SPEC 

[1]  'CHOOSE  UNIT  FOR  ANGULAR  VELOCITY ' 

[2]  ’ 1)  RADIANS  PER  SEC.' 

[ 3 ] ’ 2 ) RPM ' 

[4]  -(□tC[l]  + 3 ) x x tO=ppUNIT)**/(UNIT*U) e'  12' 

[5]  UARBOUT  7 

[6]  -*■□££’[  1 ] - 2 

[7]  'NOW  ENTER  INPUT  ANGULAR  VELOCITY' 

[8]  ' ( +)  IF  CCW  ; (-)  IF  CW:' 

[ 9 ] WI-U 

[10]  'WHAT  IS  THE  INPUT  ANGLE  INCREMENT  IN  DEGREE  ?' 

[11]  '(A  SUBMULTIPLE  OF  360  IS  RECOMMENDED)' 

[12]  /V**L  . 5 + 3 60 


V 
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V ST  ADYN  OUT \ I \MUV \S\LL\ FLAG\MUV2 

[1]  TRC-'T' 

[2]  MUV-W  ^JOINTS 

[3]  LBliGFCLG 

[4]  i'DRW'A  4 2p  'FB SCW AST'  )ZiM_TY;  ] , • 2* 

[5]  'WHAT  VARIABLE  DO  YOU ' 

[6]  TO  SEE  PLOTTED ?' 

[7]  I- 1 

[8]  ’ • ,(t  I),')  X!?  JOI/l/r  ' ,(Tl),’  ( ACTING  UPON  LINK  ' 

, ] , ' )* 

[9]  + (□£<?[  1 ] - 1 )*\MUV>I+I+1 

[10]  ’ * , (vMUV+1 ) , ' ) SHAKING  FORCE ' 

[11]  » » , (yMUV+2  ) , ’ ) IWPtfT  TORQUE ' 

[12]  ' ' , (t/VJ/K+3  ) , » ) SHAKING  MOMENT' 

[13]  ' ' , (yMUV+H ) , » ) CONTINUE  WITHOUT  PLOTTING' 

[14]  -(□£<?[  l]  + 3)*i  (0  = ppS)AA/(S«-D)  e\MUV+H 

[15]  UARBOUT  7 

[16]  ->QZ.C[  1 ] -2 

[17]  -{S-MU  17+4M0 

[18]  -(S*MUV+3ULB2 

[19]  'PLEASE  ENTER  THE  U~V  COORDS' 

[20]  'OF  REF.  POINT  FOR  SHAKING  MOMENT ' 

[21]  -(  3=pl .RPOINT-U) *ULCl  1]*3 

[22]  UARBOUT  7 

[23]  -ULCZ  l]-2 

[24]  je  ( 4 2p  'FBSCWAST'  )[f  M_TY;  ] , ' SML' 

[25]  -LB5*\  'L'eCASE 

[26]  ±(4  2p  'FBSCWAST'  )ZiM_TYi  ] , 'SM'  ,MUV2-'II'  Ll  + 'T'  eCASEl 

[27]  i'  SM'  .CASE  , '-'  ,-9SML+i'  SM'  ,MUV  2 

[28]  LBliGFCLG 

[29]  444  PLOTii' PHI' .LINKSZl ;1  2]), [.5 li'SM'.CASE 

[30]  GFACT  1 

[31]  500  620  GFTXT  'RMS=',(t>  9 3t  RMS),'  PTP=',V9  3y  FTP 

[32]  -LB 4 

[33]  LB2i-{S*MUV+2)*LB3 

[34]  GFCLG 

[35]  444  PLOTii  'PHI'  .LINKSZl  ;1  2 ] ) , [ . 5 ] j>  ' IT ' , CASE 

[36]  GFACT  1 

[37]  500  620  GFTXT  ' RMS- ' , (t  9 ~3yRMS),'  PTP=', T9  _3t PTP 

[38]  -LB^ 

[39]  LB  3 :-*■(  £B4  -2  )xi~(a/a/(i'5F'  ,C455’)=0  )aS=MP7+1 

[40]  GFCL/I 

[41]  'SHAKING  FORCE  = Z£i?<7’ 

[42]  +LB4 

[43]  GFCLS 

[44]  444  PLOT3(l'PHI'  .LINKSZlil  2 ] ) , [ 1 ] s.  3 + ( 3 *S=MUV+  1 ) <J> » F ' , ( tS 
) ,C4SJ?,  'SF*  ,CMS£ 

[45]  £54 : GFACT  1 

[46]  GFAfOF  0 779 

[47]  CMPBOi/T  13  10  10 

[48]  'PRESS  "RETURN"  TO  CONTINUE' 

[49]  LL-{ 30  -LB  1 


V 
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Cl] 

[2] 


[5] 


7 


WFWAID 

-(  QLCl  1 ] + 3 )*  i ( 0 <-l  + MUV ) a ( 3 = p 1 , MUV ) a ( 1 0 SI  \MUV ) aO  ^1  \MUV*-U 


UARBOUT  7 
"■□LC  C 1 ] -2 
WF+WF  ,MUV\_  1 ] 
WORST+WORST ,MUVl 2 ] 


V 


V WFWIN ; NJ ; I ; MUV ; Z;WF_ 

[1]  ffFCL/l 

[2]  WF+WORST+\ 0 

[3]  NJ+lipJOINTS 

[4]  1-1 

[5]  'INPUT  WEIGHTING  FACTOR  (0<VF<10)  AND  WORST  VALUE  (IN  TH 
AT  ORDER ) FOR:' 

[6]  LB1:'MAX  MAGNITUDE  OF  FORCE  AT  JOINT  ',Vl 

[7]  WFWAID 

[8]  -LBlxiNJil+I+l 

[9]  -+LB2*xSYNCASE='  2' 

[10]  'MAX  MAGNITUDE  OF  SHAKING  FORCE' 

[11]  WFWAID 

[12]  ' PEAK-TO-PEAK  VALUE  OF  INPUT  TORQUE' 

[13]  WFWAID 

[14]  'RMS  VALUE  OF  INPUT  TORQUE' 

[15]  WFWAID 

[16]  'MAX  ABS  VALUE  OF  INPUT  TORQUE' 

[17]  WFWAID 

[18]  -LB 3 

[19]  LB2 :' PEAK-TO-PEAK  VALUE  OF  SHAKING  MOMENT' 

[20]  WFWA ID 

[21]  'RMS  VALUE  OF  SHAKING  MOMENT ' 

[22]  WFWAID 

[23]  'MAX  ABS  VALUE  OF  SHAKING  MOMENT' 

[24]  WFWAID 

[25]  LB 3 :WF+-(Z+WF*0  )/WF 

[26]  WF*-Z\100*WFi  + /WF 

V 


V MPLIMITS',I\J\PAIR\VAR 

[1]  GFCLG 

[2]  t'DRW'A'*  2p  'FB  SCW  AST'  )[  tM_TY;  ] , » 2* 

[3]  J+I+l 

[4]  'PLEASE  ENTER  AN  ORDERED  PAIR  OF  NUMBERS' 

[5]  'FOR  ''MIN"  AND  ''MAX''  VALUES  OF:' 

[6]  LB1  : (VAR*-'MXYK'  [«H  ) ,LINKSlI;3  4] 

[7]  -*-  ( □ Zi  C [ l]  + 3 ) x i 3 = p 1 , P AIR+-U 

[8]  UARBOUT  7 

[9]  +ULCI l]-2 

[10]  iVAR  .LINKSLI  ;1  2 ] , ' MIN+- ' , vPAIRl  1 ] 

[11]  iVAR  ,LINKSlI;  1 2 ] , » MAX*-'  , VPA  IRl  2 ] 

[12]  ■+LBlx\A->J-*-J+l 

[13]  -+LB1*\  ( ltpLINKS)zl+-I+J+l 

V 
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V Z*SFSYN ;MUV ; FLAG;  FI RST  •,  MIN PTS ; KCURVE;  PRCPTSi  PRCPTSL  ; LIVE 
; IAINj  KIAINx  KFIRST  ;R-,KR 
[1]  GFCLAL>  TRC*'T'0  FIRST*Z*FL AG* 0 

[4]  'DO  YOU  WANT  COMPLETE  SHAKING  FORCE  BALANCE ?' 

[5]  -*■(  (□££[  l]  + 3 ) ,LB1  ,ULC\_  1]  + 1 )[  'YN'  i 1 + Q] 

[6]  UARBOUT  7 0 *ULC\_  l]-2 

[8]  -*•(  ( (X012MIN+M23MIN-X023MAX)  >0  )v0  >X01  2MAX+M2  3MAX-X02  3MIN) 

/LB 

[ 9 ] -*•(  ( (Y012MIN-Y02  3 MAX)  >0  ) vO  > Y01  2MAX- Y02  3 MIN)  / LB 

[10]  -*-(  ( (XOH3MIN+X023MIN)>0  ) vO  >X0*\  3MA  X+ X02  3MAX ) / LB 

[ 11  ] -*(  ( (Y0H3MIN+Y023MIN)  >0  )vO  >Y0H3MAX+Y023MAX)/LB 

[12]  SFC0*0  0 0 00  -*-(Z«-l)/0 

[14]  LB:' SORRY  THE  INEQ.  CONSTRAINTS  ON  THE' 

[15]  ’ MASS  PARAMETERS  CANNOT  BE  SATISFIED .' 

[16]  'PLEASE,  SELECT  AN  ALTERNATIVE: ' 

[17]  » 1)  RETROCEDE  TO  MODIFY  INEQ.  CONSTRAINTS' 

[18]  ’ 2)  CHANGE  TO  PARTIAL  SYNTHESIS' 

[19]  ->■(  ( 0 = ppMUV)  AA/MUVe  ' 12  f ) + ( 1+’ 12 ' i 1+MUV*U)$Q  ,LB  1 , □&<?[  1 ] + l 


[20] 

[22] 

[23] 

[24] 

[25] 

[26] 

[27] 

[28] 

[29] 

[30] 

[31] 

[33] 

[34] 

[37] 

[38] 

[40] 

[41] 

[42] 

[43] 

[44] 

[45] 
[47] 

[50] 

[51] 
[54] 
[56] 

[58] 

[59] 

[63] 

[64] 
[65  ] 
[66] 
[ 67  ] 
[ 68  ] 


UARBOUT  70  -*-□££’[  1 ]-2 
LB  1 :MINPTS*2  2 3 3 liMJTY] 

GFCLA 

' TO  PERFORM  SHAKING  FORCE  SYNTHESIS  USE' 

' CROSS  HAIRS  TO  FIND  INPUT  ANGLE  (IN  DEGREE)' 

' CORRESPONDING  TO  THE  HODOGRAPH' ' S POINT' 

' YOU  WANT  TO  SHIFT.  PRESS  ''A'',  TYPE  IN' 

' THE  VALUE  FOR  THAT  ANGLE,  THEN  USE' 

' CROSS  HAIRS  INTERSECT  TO  INDICATE  THE' 

' NEW  DESIRABLE  LOCATION  OF  POINT.' 

' REPEAT  THIS  PROCESS  FOR  A SECOND  POINT . '0  iO 
'NOTE:  TWO  BEEPS  MEANS  NECESSARY  CONDITIONS' 

' CANNOT  BE  SATISFIED . ' O lOO  MUV*U 

KCURVE*SFCO*  2 Opl 
CURV  E*SFTO  IAIN * 1 
PRC  PTS  L*PRC  PTS  *2  Opl 
LB 2 : GFCLG 

444  PLOT 3 PHI12  ,\_1~\CURVE 
GF ACT  444 

PRCPTS*GFLTP  PRCPTSL 
GFACT  10  DRW STAR  PRCPTS 
GFMOV  0 7790  UARBOUT  130  LINE*- 0 
-*-([]£0[l]  + 4 )*Xa/2  0=p KCURVE 

' FOR  LAST  CURVE : MAX.  MAG.  = ',(t>  8 3y  KR)0  i OO  LINE*2 
*LB3*\FLAG=0C>  FLAG* 0 

'DOES  THIS  CURVE  SATISFY  YOU?  * O'  LINE* LIN E+ 1 
i' *ULCll]+'  ,l^•(_l  + ,Y/V’^l  + □)4>,451, 

UARBOUT  70  LINE*LINE+ 10  ^□L(7[l]-30  *~Z*1 
LB3 :' USE : ''A' '-TO  TYPE  IN  ANGLE' 

' ' ' T ' ' -TO  THROW  AWAY  POINTS' 

' ''S' '-TO  SYNTHESIZE ' 

' ''I' '-TO  INTERCHANGE' 

' " Q" -TO  QUIT  SYNTHESIS' 

' AND  RETROCEDE' 
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[69]  » OTHERS-TO  SELECT  NEW  LOCATION » 

[70]  ’ FOR  A HODOGRAPH'  ' S P0INT'^>  LINE+LINE+  3 

[72]  £54  :i’ , 3+ (“3+3x98  106  117  114  116i  1 + KEEP+-GFGCUR)  4>’  LBS  I 

B9 LBO 000  LB6 LB7 ' 

[73]  LB  5 : £ 1 -+□  £ C [ 1 ] + ’ , »31» [l  + ( ( o 1 . IAIN  ) *2  +MUV ) v2  =MUV*-~1  t p PRCPTS 

]0  UARBOUT  70  -+£B4 

[76]  GFMOV  0 7790  CMFBOtfT  13,£I/VFplO 

[ 78  ] -+(  (KEEPt'lMArN)  *{.1<>KREP)  A(tf+1  )>KEEP*-l+l  . 5 + ( 1 + □ ) * 3 60  iN)  / 

□ £C[  1 ] + 4 

[79]  UARBOUT  70  LINE+LINE+20  -♦■□££[  1] -3 
[82]  IAIN^IAIN  , KEEPO  £I/VF<-£I/VF+20  +-£54 

[85]  LB  9 : — ( ( 0 ^ 1 + pFCBBFF)  a ( p 1 , IAIN_ ) = 2 +MUV+-  1+ pPRC PTS) /ULCL l1+3 


[ 86 ] UARBOUT  70  -+£54 

[88]  MUV+CURVEO  CURV E+ KC URV EO  KCURVE+MUV 
[91]  MUV-FIRST 0 FIRST+KFIRSTO  KFIRST+MUV 
[94]  MUV+RO  R+-KRO  KR+-MUV 


[97]  MUV+-SFC0 0 SFCO+KSFCOO  KSFCO+MUV 

[100]  GFACT  444 

[101]  MUV+GFPTL  PRCPTSO  PRC PT SL+KPRC PT SO  KPRCPTS+MUV 


[104  ] 
[ 108  ] 
[111] 
[112] 


[118] 


I4I^+KI4I770  KIAIN+MUV/>  -+£52 
LBO  : PRC  PT  SL-*-2  Opl^  IAIN*-  10  -+£B2 

LB1  :-+(  ( (MUV+3  ) = pl, IAIN_)  av/o  1 =#57- “it  p PRCPTS ) /ULCl 
UARBOUT  7 0 -+£B4 

DRW  STAR  2+ KEEPO  PRC  PTS+-PRC  PT  S , 2 + KE  EPO  -+£54 
£56  MINPTS = 1 + pPFOFTS;)  /□£rTTJ  + 3 
□ 4i?S£>i/r  7 0 -+£54 


11  + 3 


[120]  GF.4CT  4440  PRC PTSL+- GFPTL  PRCPTSO  GFACT  1 

[123]  -+£B3  x i 1 =i  ( 4 2p  1 FBSCWAST'  )LiM_TY;  ' , 'SFCOEF' 

[124]  UARBOUT  7 0 T+-[]5£  10  □4FSOI/27  7 0 -+£B4 

[ 128] £B8  :F£/1G+-1 

[129]  +-£52 
V 
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